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HAUPTAUFSÄTZE 


Zur Theorie der kompressiblen Potentialströmungen Ill. 


Charakteristiken, Bicharakteristiken und Eikonal der linearisierten 
Grundgleichung. 


Von M. Pinl (Messerschmitt-AG.) in Augsburg. 


Es handelt sich um die Integration der Differentialgleichungen der Cha- 
rakteristiken und Bicharakteristiken, die zur Grundgleichung für das Ge- 
schwindigkeitspotential einer ebenen adiabatisch-kompressiblen Strömung 
gehören. Diese Grundgleichung wird in der mit Hilfe einer verallge- 
meinerten Minkowskischen Stützfunktion &(a,ß,y) linearisierten Form 
(Gleichung (1,4) vorausgesetzt. Die Integration der Charakteristiken und 
Bicharakteristiken kann nach einer Reihe von Transformationen der un- 
abhängigen Veränderlichen a,ß,y auf Quadraturen zurückgeführt werden. 


$ 1. Einleitung. 


Der Ortsvektor der Geschwindigkeitspotentialflächen einer ebenen adiabatisch-kompres- 
siblen Strömung läßt sich nach einem früheren Ergebnis’) als Gradient eines Skalars 
(2(a,ß,y) darstellen, wodurch mit Beschränkung auf homogene Lösungen erster Ordnung 
die quasilineare partielle Differentialgleichung zweiter Ordnung 


A+lu—D) DL. + u PD] Dex —2 Di D, Day + +au DDP, +uP] Po, =0 .: - (1 


für das Geschwindigkeitspotential ? (x, y) in die lineare Differentialgleichung zweiter Ordnung 


Iua+ 233 + = - — ee . ie 5 
k: (u I) + pP) +iy 


übergeführt wird, wenn man partielle Ableitungen nach x,y bzw. a,ß,y durch Indizes zum 
Ausdruck bringt. Dabei bedeutet 4 das Quadrat der Schallgeschwindigkeit im ruhenden 
Br ‚ . ne 

mit x als adiabatischem Exponenten. 


Medium und « die Konstante —, 


1) H. Behrbohm und M. Pinl: Zur Theorie der kompressiblen Potentialströmungen ], Neue Linearisierung 
der Grundgleiehung der ebenen adiabatisch-kompressiblen Potentialströmung: Z. angew. Math. Mech. Bd. 21 (1941), 
S. 193 bis 203 
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KEMEENE N 
Betrachten wir statt P(x,y) die Funktion z (x, ») = yi P(x,9y), so geht (1,1) mit den 


üblichen Bezeichnungen 


_02 A ' 0’ 2 
Br 0x Er y' i we 7. h dxody’ o® y 
in die Differentialgleichung 
+ —Dp+tueflr-2pqas+li+lu—Df+uplt=0.. .. . (13) 


und (1,2) in die Differentialgleichung 





= = -y—=0 De IR 
(u 1)(@® +) Hy?" 


(r Ouat @83Tt 











über, die wir nunmehr als Grundgleichung für die in Rede stehenden Probleme ansehen. 
Der Grundgleichung (1,4) ordnen wir die folgenden für ihre Integration wichtigen Hilfs- 
differentialgleichungen (bzw. Hilfsformen) zu: 
a) die Charakteristikenbedingung für einen Ausnahmestreifen erster Ordnung €, °): 


dp\’ 


| =0 Eu te - 


n,=(y2) +6; 


+ an al 
Ja dp (u — 1) (a? + B°)+ y’\0y 

b) die Charakteristikenbedingung für die Grundmannigfaltigkeit €, von C, (Hamilton- 
sche Differentialgleichung des Bicharakteristikenkomplexes’)): 


oy\’  (Oy\ ze 
H, Irre De+tm+r RE... 


e) die ternäre Riemannsche Metrik®): 


(u — 1)(a? +?) 
. ( 


ds®’=da +dp+' 1 ae ERBE ; 2, 


‘ 


d) die Mongesche Differentialgleichung des Bicharakteristikenkomplexes?): 


(u —VD)l(a’+Pp?) ,. dp dy 
Be y® ' "de ’ da (18); 
e) die Differentialgleichüng für den geodätischen Abstand oder das Eikonal 7°): 
oI? /dT? (ua—1)(a+P°)+y?/oI\ 
(—) +45) +" SEITE | -1=0 (1,9). 
da \0Pp y” \0y 


/ 


Das Ziel der vorliegenden Untersuchung ist die Integration der Hilfsgleichungen (1,5), 
(1,6), (1,8) und (1,9). Diese Integration gelingt (im Falle der Gl. (1,9) und (1,5) bis auf (Qua- 
draturen) nach einer Reihe von Transformationen der unabhängigen Veränderlichen, an deren 
Ende wir die transformierten Gleichungen (Grundgleichungen wie Hilfsgleichungen) mit einem 
Querstrich charakterisieren: G=0, H,=0, H,—0, M=0, E=0. Dabei ist der konform- 


2) Courant-Hilbert: Methoden der mathematischen Physik; Bd. TI, Kap. VI, $ 1: S. 346 bis 356; Berlin 1997. 

3) Vgl. 2), insbesondere S. 353 bis 354. 

4) Vgl. ?2) Kap. VI, $ 2, S. 356 bis 367; $ 9, S. 430 bis 433; den niehteuklidischen Charakter der Metrik (1,7) er 
kennt man am besten durch Berechnung des Riemann-Christoffelschen Krümmungstensors R,.3,g9 in Zylinder- 
koordinaten %,, %a, 43 nach GI]. (2,1). Man erhält 

ER, w+u—l)u ( u”) u 

Rızı2 = Rızıs= Rıaazs = Rıaza = 0, Rısıa= +, 3, Raszs= ("-—1) 
"3 u. + le — 1) u 
3 1 

Daher gilt im Falle «= 1 für die Klassenzahl k der quadratischen Differentialform (1,7) k=1. Da das System 

Rıaı2 = QRıı Qa2 — Rıe Rı=t, Rısıs = Rıı Rıs — Rıs Rai = 

Rı22s = Rı2 Rs — Qıs Rat, Rısıs = Qıe ( 

Rısıs = Rıı Ass — Qıs RsızEV, Rosas = Rn ( 

auf Widersprüche führt, gilt nach den Kriterien von K. H. Weise (Beiträge zum Klassenproblem der quadratischen 

Differentialformen, Math. Ann. Bd. 1104 (1994), S. 522 bis 570) k l. Indessen ist die Einbettung bereits in einem 

fünfdimensionalen euklidischen Raum mit den kartesischen Koordinaten X;, X2, X3, X4, X möglich, wie die von Herrn 
K. H. Weise angegebenen Gleichungen 


233 — Rıs Ra =) 
I38 


Qıs Ar=FV 


x =u 608 (Vu — ı1g us)ecosu, Xu sin (Va — ı lg us) cos un, 
X, Hg 
x = Ur eos (Vu — ılg us)sin us, % = u, sin (Vo — ılg us) sin us, 
zeigen. Für die Klassenzahl gilt also k =2. 


5) Vgl. ?). 
6) Vgl. 2) Kap. Il, $ 9, S. 98 bis 101; Kap. VI, 9, S. 431 bis 433. 
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invariante Charakter des Bicharakteristikenkomplexes (1,8), der ja den isotropen Komplex’) 
eines dreidimensionalen Riemannschen Raumes mit der Maßbestimmung (1,7) darstellt, von 
besonderer Bedeutung. Er sichert die Kovarianz der Transformationsprozesse: 


G=0>6G6—=0 G=-0>G =0 G=0>G =0 
Y Y bzw. Y Y bzw. Y Y 
M=0>M=0 H,=0>H,—0 H,=0>N, —0, 


Das Integrationsresultat der Gl. (1,8), (1,6) und (1,5) ist somit konforminvariant cha- 
rakterisiert und unempfindlich gegenüber Multiplikation der Grundgleichung mit einem nicht- 
verschwindenden Faktor. Demgegenüber ist das Integrationsresultat der Eikonalgleichung 
(1,9) metrisch-invariant charakterisiert und gegenüber Multiplikation der Grundgleichung mit 
einem Faktor im allgemeinen nicht kovariant. Insbesondere erhalten wir, wenn die Grund- 
gleichung nach einer Parametertransformation einen Faktor abspaltet, unter Umständen das 
nichtkovariante Schema: 


G—=0 > G—=0, 
Y Bi E=E E*. 
E=V>E=0; H=0, 


Während also die Reihenfolge von Zuordnung (zur Grundgleichung) und Transformation 
(auf neue Parameter) für die Integration der Differentialgleichungen (1,8), (1,6) und (1,5) be- 
langlos ist, ist dies für die Integration der Eikonalgleichung (1,9) nicht der Fall. Dort führt 
vielmehr der Weg G=0>G=0>E*=0 zu einer einfacheren Gleichung als der Weg 
G=-0>E=0>E=0O (vgl. $ 3). 


Für die Integration der Gl]. (1,8), (1,6) und (1,5) wählen wir den Weg 6 =0>6=0 
>M=0 bzw. G=0>G=0>H,=0 bzw. G6=0>G=0>H,=0 und benutzen außer- 
dem den dualen Charakter der Mongeschen und Hamiltonschen Differentialgleichungen °) 
(vgl. $$ 2 und 4). Dabei handelt es sich in allen Fällen stets um dieselben Transformationen 
der unabhängigen Veränderlichen. 


$ 2. Zyklische Parameter und Bicharakteristiken. 


Die Koeffizienten der Differentialgleichung (1,4) enthalten sämtliche Parameter explizit. 
Gelingt eine Reduktion der Anzahl dieser explizit auftretenden Parameter durch eine Trans- 
formation der Differentialgleichung auf neue Parameter, so handelt es sich mit Bezug auf 
diese in den Koeffizienten der transformierten Differentialgleichung explizit nicht mehr auf- 
tretenden Parameter um sogenannte zyklische Parameter®). Benutzen wir die Zylinder- 
symmetrie der Koeffizienten der Grundgleichung (1,4) in den Parametern a, ß,y, so gewinnen 
wir durch die Transformation 

. + > ) 
a=u,cosu,, P=usmu, Y=uU; uv=Va+P’, w=arcig -, u,=y.. &) 

3 ’ . > a 
aus (1,4) in Zylinderkoordinaten «,,%,., die Differentialgleichung 
1 u; 1 
| 3 22 
7) > — 0. T ua - Olalis T BE. 5 am, 
WETTE u em 


deren Koeffizienten den zyklischen Parameter «, nicht mehr enthalten. Eine weitere Re- 
duktion gelingt durch die Transformation '") 


‚ ptl2 . — nVR. . 1 en u a 9% 
„met, weU, wet; u,=7, % 2 us 7 VE © 
3 
Sie verwandelt die Gl. (2,2) in 
on R . I+{u— Ver ‚1+(u+VDv} 
e-2":l(l + ur) orıv, Er Ovire Tr Orr tr A Vvsr; T v Wwv; rd, 
1 1 l 


spaltet also den Faktor e=-2'» ab. Es verbleibt die reduzierte Differentialgleichung 


”, D. h. die Nullinien der quadratischen Differentialform (1,7). 

*s) Vgl. ?) Kap. II, $ 4 und 5; S. 70 bis 74. 

9, G. Prange: Die allgemeinen Integrationsmethoden der analytischen Mechanik, Eneyklopädie der mathe- 
matischen Wissenschaften IV 1, II, Heft 4; $ 9, S. 571 bis 580; Leipzig 1935. 

10) Die Kenntnis der Transformationen (2,3) und (2,5) verdankt Verf. einer freundlichen Mitteilung des Herrn 
K. H. Weise. 


20% 
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' P En ‚i+@-1be ‚I+@+Dei _0 024 
(i+ uv) Orr, av, Ovıta tt Orvavg 2. Ovzvg 2 dv, Di, = (2,4), 
1 1 


deren Koeffizienten nurmehr den einen Parameter v, enthalten. Schließlich führt die Trans- 
formation '') 


| 
>, w,„—5-gl+umw),n,=m;; w, 1, Wm=rT,, g(il+uvi),w,=v, (25) 
’ BR” ’ 2 Be 


wiederum auf orthogonale Parameter, in welchen die Differentialgleichung (2,4) die Gestalt 





— A 1+(u—1)w; l+(u— 1) w; 
Tr (1 Tu m) Odw,w, 7 9 Ows wa T VDwzwz 
' l+umw DR 
(2,6) 
1+(u+1)w! 1— (u—]1) w; 
F ; Ow, Tr = Ow, =V 
w, 1-+ um; B 











annimmt (,, w, als zyklische Parameter). 

Mit den Transformationen der Grundgleichung (1,4) sind die 
Mongeschen Differentialgleichung (1,5) des Bicharakteristikenkomplexes kovariant verknüpft. 
Man erhält an Stelle von (1,8) kovariant zu (2,2) die Mongesche Gleichung 


Transformationen der 





na Dit ,„ u ' * N.” en 
l1+uWu;+ - u —=0, u, u,;: ne ie 
> u; ’ du, ’ de, 
kovariant zu (2,4) die Mongesche Gleichung 
Bun, hal dev, s dev. 
| +2 v, v’, + (Al La v*) u v v” e 0, v, en = v R n 3 £ J . ’ ? B (2,8). 
2 4 n :: de, dv, 
und schließlich kovariant zu (2,6) die Mongesche Gleichung 
— 1+(u—-D)w? , dw i dw. 
M Ü tl +umwm)w?+wn?=0, w,= ’,w,= u. .., MM 
l+umw 1 . ’ de, de, 








Offenbar spielen auch in diesen Gleichungen die Parameter «,,r,,®,,%,,, die Rolle 
zyklischer Parameter. Das gilt ebenso für die den Mongeschen Gleichungen (1,8), (2,7), 
(2,8), (2,9) dual zugeordneten partiellen Hamiltonschen Differenti: ılgle ichunge :n erster Ord- 
nung. Wegen der Kovarianz dieser Zuordnung können wir unmittelbar an die Grundgleichung 


(2,6) anknüpfen. Ihr entspricht die Hamiltonsche Gleichung 














en m, dw,”  (d )[1 1) w; 
M,=-(+ um: (5 ') +1 +0 Dm]l, .) + Fir+e- Dil. am 
dr, n 
Die Integration dieser Gleichung führt nach Division durch den Faktor 1 + (u — 1) w! =E0 
und Separation der Variablen auf das vollständige Integral 
5 l+(u-1)&® 1+(u-1D)E# |, 
==60 + Ve r + F & 2 
w, n ;\] 0 (i+u®e (+ ® BE. 5. SER 
bzw. auf das allgemeine Integral 
Ri: l+u—- DE 1+a—1)# 
= a: 2 sis “ - £ 7 ») 1» 
wi; om,til] ” At) A+u2)& ua ER. 


wenn wir den Parameter r als willkürliche (analytische) Funktion des Parameters o ansetzen. 


Das Integral 


‚ira—)® ‚I+ 


1 R ’ 
J(w, o-\y- (tue +üirup) ee A 


11) Vgl. 10), 
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läßt sich durch quadratische Irrationalitäten und elementare Transzendente in geschlossener 
Form darstellen‘). Das Ergebnis ermöglicht nunmehr die integrallose Darstellung des 
jicharakteristikenkomplexes mit Hilfe der willkürlichen Funktion 7 (6) und ihrer ersten 
beiden Ableitungen, in die das Integral (2,13) als bekannte Funktion der oberen Grenze , 
eingeht. 


Für die Bicharakteristiken ergibt sich die explizite Parameter- 
darstellung 





| e.—w,, m, =w,(W;; =F iJ, +T, w.=w,(w,; 0)=H "re: ioJ, tor | . (2,14), 





für die Fokalkurven des Komplexes die implizite Parameterdarstellung 





[77 


|: iJ+ow,— w,=t(o), +iJ, +w,=tr' (0), tiJ ss =T" (0) 





(2.18). 
0.J %J BE: d’r 


[7 


„= : A — a. we A - - 
o do 0° 0 0° do do? 


In (2,14) erscheinen die Bicharakteristiken als Schnittkurven benachbarter Integ: “- 
flächen gegeben (Parameter w,), in (2,15) die Fokalkurven als Gratlinien aller aus dem voli- 
ständigen Integral (2,11) entspringenden einparametrigen Integralflächenscharen (Parameter 0). 
Eine explizite Darstellung auch dieser Kurven in der Form 


„ ’ [77 


/ ’ ’ ’ ’ 
m, w.(0,7,7T,7 ), W=W,10,7,7,T ), Wm,=w,(0,7T,7T,T) 


setzt die Möglichkeit der Auflösung der GIn. (2,15) nach w, voraus. 


S 3. Das Eikonal. 


Die der Grundgleichung @=0 zugeordnete Eikonalgleichung E=0 (1,9) verwandelt 
sich durch die Transformation (2,1) in die Differentialgleichung 


oI® 1/9I® u-Dwtu/oöl® 
tt ” (du. BR 7 2 re 


sodann durch die Transformation (2,3) in die Differentialgleichung 


om: I /d NM? Fa oI\? Se PY; 
(3 - u: (3 3 +[1 + (u Dell; u an) e- *a 0 J M r , ö (3.2). 


schließlich durch die Transformation (2,5) in die Differentialgleichung 


— 128 im, 4 Eier 
ü o0w, 1+umwidm, nv; \0 = 
i (3,3) 
1+(u l)w! 01 01]? zw —-—Ig(l+uw;) 
+-[1 u—1)w?]| a mw e e 0 
I+( ) || l+um! dm, 'd | 


Dagegen lautet die der transformierten Grundgleichung @ = 0 (2,6) zugeordnete Eikonal- 
gleichung 





1 2} l+ um! [277 w; +E , 


E* ; +: ; + ; 
l+umi\ow, 1+(u —1) wi \d - 1+ (u — 1) mw} \o u, 


1=0 .. 84). 











Der Parameter w, ist also zyklischer Parameter der Differentialgleichung (3,4), aber 
nicht ein solcher der Differentialgleichung (3,3). Die Vertauschung der Reihenfolge von Zu- 
ordnung (der jeweiligen Gleichung) und Transformation führt zu verschiedenen Eikonal- 
gleichungen (E== E*). Dieser Unterschied wird durch die Transformation (2,3) verursacht, 
nach welcher sich von der Grundgleichung (aber nicht von der Eikonalgleichung) der Faktor 
e"?v» abspaltet und so für die Grundgleichungen (2,4) und (2,6) aber nicht für die Eikonal- 
gleichungen (3,2) und (3,3) einen weiteren zyklischen Parameter v, bzw. w, liefert. Für die 
Integration wird man natürlich der Eikonalgleichung (3,4) den Vorzug geben, da deren Ko- 





effizienten nurmehr von dem einen Parameter w, abhängen. 
12) Zunächst führt die Substitution S?—=n auf die Zerlegung 2 J= J; — "Ja des Integrals J, worin die Integrale 


Jı und Js nurmehr quadratische Irrationalitäten enthalten; die weitere Berechnung vermittelt die Substitution 1 + « i=L, 
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Wir setzen © 2 = 1 betrachten das Syst 
ir setze = : =», ‚, und betrachten das System: 
EBENEN Zn, Di7 dw, nahe Te A nn TER, 
Y 1 3 1l+ um s Dr 4 f 0 
ar Pr sp,” =VU), 
ya I+ um: 1+ (u HarP: 1+(u Intl» 
(3,2). 
X, Ps a,=0, | 
X, P, —a,=0 


Seine Poissonschen Klammern?) verschwinden identisch: 


3 X,oX, 9X,9X, 


2.2. 0, 
ein Zu VO PRdm; Opr = 
k=1 
V.X Se X,0X, 0X,0 X) 
(A,,A;) £, \öpr dwı px dm; - 
k=1 
5 i \"/0 X,0X, 9X,0X, 
(A,,X,) R AwT) N N 
u \O Pr! Wi; OPrO Wr 
1 


k 
Das System (3,5) ist also ein Involutionssystem. Für die Ableitungen p,,P,,P, ergibt sich 


(1-+umj)’ „ww (l+um;) 


_ — 2 — - 9 4,8 Be er 
+ (u—1)m} ’I1+(u—)m! P: 2, Ps a 


B} 


»,=+ l+u m a; 
und damit das totale Differential 


1-+ um?)? .„w(l+-um?) m 
I ı — a dw +a,dwr,+ta,dw, .. (N 


u +(u—1)n} 474 (u—1)w; 





l+umwi—a; | 
und das vollständige Integral 


(1+ um))’ 


wi (1-+ um}) ER 
Era 7. : —  dwtan,+aw, . . (38. 
+ (v ) mw? 


1+(u— 1)w; 


i=6# \ | l+uw! —a a; 
Sind dann h Relationen W,, (a,,«a,,a,)=0 zwischen den Parametern a,,a,,«a, des voll- 
ständigen Integrals (3,8) gegeben (h< 3), so führt die Elimination der Parameter 4, aus 
dem System 
h 
ol ’, 8 W 
3 u, —!.9 y=1,23, h<3 

OB, u oa, 
k1 
auf drei Gleichungen zur Bestimmung der Parameter a,,a,,a, als Funktionen von w,,m,,w, 
und damit in bekannter Weise'*) auf das allgemeine Integral. Damit hat sich ergeben: die 
Integration der Eikonalgleichung (3,4) der linearisierten Grundgleichung 
(2,6) läßt sich auf eine Quadratur zurückführen. 


$ 4. Die charakteristischen Mannigfaltigkeiten. 


Jede einparametrige Lösungsschar y=y (a, ß;c) der Differentialgleichung (1,6) führt in 
der Form g (a,ß,y)=e auf Lösungen der Differentialgleichung (1,5)'°). Um zu einer direkten 
Integration der Differentialgleichung (1,5) zu gelangen, knüpfen wir an die transformierte 
Grundgleichung (2,6) an und ordnen ihr die charakteristische Gleichung 


0g 


1-+-(u— 1)w; (04 
7 | B ...: 


l+(u iu be : 
l1+ um; dm, | 


H,=(A+u wi) | ) + ) + Pen \O m. 


13) G. Hoheisel: Partielle Differentialgleichungen; Kap. IV, $ 2, S. 86 bis 92; Berlin 1928, 
14) Vgl. 13), 
15) Vgl. ?2) Kap. VI, $ 1; S, 346 bis 356, 
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zu. Sodann setzen wir ähnlich wie bei der Integration der Eikonalgleichung in $ 3 ES 
- dw, 
09 og h 
dr 5 w, 9» 5 "Goa 4, und betrachten das Involutionssysten 
1+(u— 1)w! 1+(u— 1)w] 
3 l+umwi)g + ' + n \g? 0, 
ini ui 7 Ya l+um; 1 w; 4 
® ua ar (4,2). 
r, Ga N EEE 7 ED 
Y, 4a b, 0 
Für die Ableitungen g,,4,, 4, ergibt sich 
’ 1+(u—1)w} „1+(u—1)w! 
A b? “ —, — + b5 \ > 2 (4,3) 
’ (d-+um;)’ ’(1+ uw?) w; 
und damit das totale Differential 
1+(u— 1) mw? 1+ (u — 1) w? 
Io si + bi — de, +b,dw,+b,den, . . . 4 
” 5 ” (1+ umw;) ” wi (1+ um;) rannte, (4,4) 
und das vollständige Integral 
1+ (u — 1) w? l+(u— 1)»; 
o=b, Hill 1 — +6 "dw, +b,w,+b,w. . (4,5) 
’ ns | £ 1l+um) £ wi (1+ um?) WERTE  - (4,8). 


Damit ist auch die Integration der charakteristischen Differential- 
gleichung (15) auf eine Quadratur zurückgeführt. 

Weiterhin entsteht die Frage, ob es nicht möglich ist, die hier und früher auftretenden 
(Quadraturen eventuell nach weiteren Transformationen in geschlossener Form zu be- 
rechnen. Im Zusammenhang damit erscheint auch eine nähere Untersuchung der Konform- 
krümmungseigenschaften der quadratischen Differentialform (1,7). von Interesse '®). 


Zusatz während der Korrektur: Die quadratische Differentialform (1,7) ist nicht konform- 
euklidisch, da der ihr entsprechende kubische Tensor 


2 D.w= Di bus» — Du Las; 4, u, v=1,2,3 


nicht identisch verschwindet. Dabei bedeuten: 


1 
L ‚—=Ls Rus E= n R Qu»; R.ı - R; "g B= Burg’; 


Ö | o | o | i ß 
D; Lu: I; ı u Los i In»: a. MM m DB 1.2. = 

0%, lu ıJ d r r E 
Die Berechnung dieser Tensoren vereinfacht sich beträchtlich, wenn man nach einer Bemerkung 
des Herrn G. Herglotz unter x, Parameter versteht, welche die Differentialform (1,7) durch 
die Transformation 


’ 


ame-Aossa,, P=errhsne,, Yy=z, 


erhält. Der Hauptteil der Grundgleichung (1,4) kann demnach niemals auf die Laplacesche 
Form ©, x + @xsxs + ®x3 x, transformiert werden. 118 


16) Auch die in der Einleitung erwähnte notwerdige Beschränkung anf h« mogene Lösungen erster Ordnung der 
Grundgleiehung (1,4) führt auf eine bemerkenswerte Konsequenz. Für die zweiten Ableitungen einer solehen Lösung 
bestehen notwendig die Homogenitätsrelationen 


aa + 303 + 179,0, 
. (CO u 
. 0) ( = — 0 
ew er E d(a,3,Y) 
ao, 0,37 ©, 0, 


Bereehnet man daraus «, ?,y als Funktionen der zweiten Ableitungen und setzt die so hervorgehenden Ausdrücke in 
(1,4) für «, ?,y ein, so ergibt sich die Differentialgleicehung 


2 2 
(vaa MER Bin ;) 
[0 PR W242 + Vu, ==0. 


6) y 
k 2 N t DB. L Be 22 „2 
(Oaa 33 Or ;) 4 (u — 1) [(Wg3 03, Vay®3B) (ga 03, ®ay@3a)?] 


Ihre Koeffizienten sind konstant, der lineare Charakter der Grundgleichung (1,4) ist aber jetzt wiederum zerstört. 
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Die Ausbreitung von Biegungswellen in Stäben. 
Von W. Flügge in Berlin-Adlershof. 


Es wird gezeigt, daß auch in einem auf Biegung stoßartig beanspruchten 
Stab elastische Wellen ron endlicher Laufgeschwindigkeit möglich sind. Ihre 
(Geschmindigkeiten werden bestimmt und ihre Eigenschaften erörtert. 


1. Longitudinalwellen. 


Wir knüpfen zunächst an Wohlbekanntes an und betrachten einen geraden Stab von 
konstantem Querschnitt F, der bei 2=0 beginnt und sich in positiver Richtung ins Unend- 
liche erstreckt. Bis zur Zeit {=0 sei er unbelastet und in Ruhe. Läßt man von t=0 an 
auf seinen Anfangspunkt eine konstante Zug- oder Druckkraft P wirken, so durchläuft pr 
bekanntlich eine Dehnungswelle, die sich mit der Schallgeschwindigkeit e=yE/joe=yEF/] 
fortpflanzt (E = Elastizitätsmodul, o = Massendichte, Fo= u). Vor dem Kopf dieser Welle 
P 
ue 
Die Differentialgleichung dieses Vorganges ist bekanntlich dieselbe wie für die longitudinalen 
Schwingungen des Stabes: 


ist der Stab noch in Ruhe, hinter ihm hat er die konstante axiale Geschwindigkeit u = — 


BE ee ee ee 


wobei Striche Ableitungen nach .r, Punkte Ableitungen nach £ bedeuten. 


Ganz entsprechende Vorgänge wie beim Längsstoß spielen sich ab, wenn man den Stab 
durch ein Torsionsmoment belastet oder wenn man einen gespannten Faden seitlich auslenkt: 
Immer kommt man auf eine Differentialgleichung von demselben Bau mit entsprechend ge- 
änderter Bedeutung der darin vorkommenden Größen. Immer gibt es eine endliche, für alle 
Störungen gleiche Ausbreitungsgeschwindigkeit e, und es ist ce? gleich dem Quotienten aus 
einem Parameter, der die Rückstellkraft beschreibt (Torsionssteifigkeit @ Ir beim Torsions- 
stoß, Spannkraft S beim Faden), und einer Trägheitsgröße (Drehträgheit o /, beim Torsions- 
stoß, «u beim Faden). 


2. Biegewellen nach der einfachen Theorie. 


Ein wesentlich anderes Bild bietet sich, wenn man statt der Axialkraft P eine Biege- 
belastung anbringt, also eine Querkraft Q oder ein Biegemoment M. Dann liefert die Balken- 
biegetheorie für die Durchbiegung » die Differentialgleichung 


cl ae ae 





die nur noch eine Ableitung der höchsten Ordnung enthält. Zu ihr gehören keine laufenden 
Wellen mit einem ausgeprägten Wellenkopf, vor dem der Stab noch in Ruhe ist. Das wird 
schon durch eine Dimensionsbetrachtung nahegelegt, denn EI/« hat nicht die Dimension 
eines Geschwindigkeitsquadrates, und es wird durch die bekannten strengen Lösungen be- 
stätigt. Stützt man z. B. den Anfangspunkt des Stabes («=0; w=(0) und läßt man dort von 
t=( an ein konstantes Biegemoment M (0,1)—=M wirken, so heißt die Lösung der Gl. (2)'): 


en 0 M c (s( N v3) c( na ) 


EI n® Vat 
Zen a BEN (3), 
n” 3° ‚inxz ‚na . nm X n”x\| 
TR s(, u P ut ls u 


wobei 
4 


u 
- — 
4EI 
ist und S und € die Fresnelschen Integrale bezeichnen: 


BR. RR. DR 
S(z2)= \ sin a C(z)= \ cos dl. 


0 0 


1) Die Lösung findet sich bei W. Prager: Ing.-Arch. Bd. #4 (1933), 8. 30, weitere elementare Lösungen bei 
H. Baumann: Ann. Phys. V. F. Bd. 24 (1955), S. 49, 
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Dies Ergebnis, das in Bild 1 dargestellt 

ist, zeigt, daß sofort beim Aufbringen 

der Last der ganze unendlich lange Stab 

in Mitleidenschaft gezogen ist und daß v.Eyr 
nun überjeden seiner Punkte viele Wellen 
hinlaufen; denn die in Bild 1 gezeichnete 
Durehbiegungsverteilung streckt sich pro- ß 2 
portional y f, wobei sich außerdem die ! 
Ordinaten proportional f vergrößern. Dies 

Ergebnis ist eine zwangläufige Folgerung 

aus der Differentialeleichung (2), die M 
natürlich wie alle Elastizitätstheorie ein- 
schränkende Voraussetzungen enthält, 
aber doch nur solche, die sich insbe- 
sondere bei schlanken Stäben überall 
bewährt haben, und von einem unend- Bild 1a und Ib. 

lich langen Stab muß man wohl an- 

nehmen, daß er schlank genug ist, um die Voraussetzungen der Gl. (2) zu erfüllen. Es ist 
also zu erwarten, daß die Lösungen von (2) das Wesentliche der Vorgänge, die sich beim 
Biegestoß abspielen, richtig wiedergeben. Trotzdem besteht gegen diese Lösung ein grund- 
sätzlicher Einwand: Daß sich die Störung, wenn auch. zunächst nur mit sehr kleiner Am- 
plitude, unendlich schnell über den ganzen Stab ausbreitet, paßt nicht recht zu der Vor- 
stellung, daß ein Biegemoment doch aus Zug- und Druckspannungen besteht, von denen man 
erwarten muß, daß sie sich längs der Zug- und Druckzone des Stabes mit der Schall- 
geschwindigkeit Y E/o ausbreiten werden. 





: wa N LI 





3. Querkraftdeformation. 


Zur Aufklärung dieses überraschenden Verhaltens der Biegewelle gehen wir von der 
eben gebrauchten Vorstellung einer Zug- und Druckzone aus. Wie Bild 1 zeigt, hat das Biege- 
moment viele Nullstellen; es gibt also nicht einfach eine Zugzone unten und eine Druckzone 
oben, sondern zwischen Zug und Druck findet auf jeder Seite des Stabes ein häufiger Weclısel 
statt. Und eben dieser Übergang z. B. der Zugkräfte von der einen Randzone in die andere 
und wieder zurück ist das Entscheidende. Er wird bekanntlich durch die Querkraft vermittelt, 
und die unendlich schnelle Ausbreitung der Biegestörung liegt daran, daß die übliche Balken- 
theorie die Querkraftdeformation vernachlässigt. Das bedeutet, daß der Stab gegenüber einer 
speziellen Kraftwirkung unendlich steif gemacht wird, und da die Laufgeschwindigkeit von 
elastischen Wellen proportional der Wurzel aus der jeweils in Frage kommenden Steifigkeit 
ist, so muß damit in der Tat eine unendliche Laufgeschwindigkeit hervorgerufen werden. 
Das schließt natürlich nicht aus, daß die Werte für die Durchbiegung, die die übliche Theorie 
ergibt, mit guter Annäherung richtig sind, da ja die Querkraftdeformation ihrem Betrage nach 
in einem schlanken Stab sicher keine wesentliche Rolle spielt. 

Nimmt man die Querkraftdeformation in die Gleichungen auf, so ist die Formänderung 
des Stabes nicht mehr durch die Durchbiegung w allein zu beschreiben; denn da jetzt die 
Quersechnitte nicht mehr senkrecht auf der verbogenen Stabachse stehen, so ist ihre Drehung y 

nicht mehr gleich d9w/d2= mw’ und wird zu einer zweiten 














Pr selbständigen Verschiebungsgröße (Bild 2). Das Elastizitäts- 
9 gesetz besteht jetzt aus zwei Gleichungen 
7 en M EI ’  T /..t 
BT nz M=—K#Hlvy, Q): (7 F, (w y) e s s (4a, b). 
alf —t, mn x e - . : i 
\ | Be I mmar Die dynamischen Gleichungen heißen wie bisher 
u u ’ .. 7 - 
u Q=uiw, Er. ° 
| 0+0ax Aus 5a) ergibt sie ch; Wlieminaki ee 
en E; Tv, To Aus (4b), (5a) ergibt sich durch Elimination der Querkraft 
— iz e .; Wh . 
v=Ww ee a 
Bild 2. ! GF. a) 


und aus (4a,b), (5b) durch Elimination von M und e) 
a Me ee: ; ME re © 


Das sind zwei partielle Differentialgleichungen für w und y, aus denen wir y eliminieren, 
um auf eine mit (2) vergleichbare Gleichung für w zu kommen. Dazu brauchen wir nur die 
zweite einmal nach x zu differenzieren und dann y aus der ersten einzusetzen: 
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Das ist eine Gleichung, in der zwei verschiedene Ableitungen der höchsten Ordnung 
vorkommen und die deshalb auf Sprungwellen mit endlicher Fortpflanzungsgeschwindigkeit 
führt. Man sieht das am einfachsten, wenn man an der Stelle #=0 zunächst eine periodische 
Störung w (0,#)=asinot# einleitet. Dann heißt die Lösung der Gleichung 


w(at)= —a sin | |, 


und durch Einsetzen ergibt sich: 


Daraus folgt 
1 Eu | u Il u 
ce? 2GEF, x | 2G Fr) BI 


Das ist das reziproke Quadrat der Laufgeschwindigkeit der Sinuswelle. 
Trägt man es über 1/®” auf, so erhält man den in Bild 3 gezeichneten i 
Verlauf. 1/e? hat also ein Minimum, e ein Maximum für 1/0’=0, und 


diese Geschwindigkeit ist 





i 
GF, u [ 
Ce i F 
nax z \ EM 
Schneller kann sich keine laufende Welle irgendwelcher Frequenz aus- 2 a 





breiten, und e,„.„, wird also die Geschwindigkeit sein, mit der eine am 
Anfang des Balkens eingeleitete Sprungwelle diesen durchläuft. Das «| 
heißt aber noch nicht, daß vor der Sprungwelle überhaupt keine Störung Bild 3. 

in dem Stab vorhanden ist. Ein solcher Schluß wäre zwingend, wenn 

es nicht auch noch imaginäre e gäbe. Solche liefert der untere Ast der in Bild 3 gezeichneten 
Parabel. Zu. ihnen gehören reelle Lösungen vom Typus 


VON C 


w(a,t)= et'"* *sinmt. 


- 
Es gibt also neben der Ausbreitung der Störung mit endlicher Geschwindigkeit immer noch 
eine unendlich rasche Ausbreitung. 


4. Rotationsträgheit. 


Wie das Gegenbeispiel der Wärmeleitung zeigt, muß sich zwar durchaus nicht jede 
Störung in Form einer Wellenfront mit endlicher Geschwindigkeit ausbreiten, aber beim Biege- 
stoß gegen einen prismatischen Stab verschwinden die Lösungen, die eine unendlich schnelle 
Ausbreitung darstellen, völlig, wenn man die begonnene Ergänzung der Differentialgleichung 
konsequent zu Ende führt. 

Wir haben bisher als einziges Trägheitsglied das Glied «w in un 
der Gl. (5a) für die zum Stab senkrechten Kräfte. Damit ist aber, 

Ir 
| 





wie aus der Theorie der Biegeschwingungen bekannt, der Trägheits- 
einfluß der nicht auf der Stabachse liegenden Massenteilchen noch 
nicht erschöpfend dargestellt. Diese erfahren ja nicht nur eine lot- 
rechte Verschiebung w, sondern, wie Bild 4 zeigt, auch eine waage- 
rechte u—=zy. Dazu gehören Trägheitskräfte, zu deren Überwindung / | 
ein Moment ol; erforderlich ist, und die Momentengleichung (5b) \ 
heißt also vollständig 42 


Eee een Bild a. 





Führt man damit dieselbe Rechnung durch, so erhält man die Differentialgleichung 


r ” u .. u :.; - 
ET w!\ [EI — t-ollv + -— olu+unm=0 .....2.0. 
"I Au 22 ! 


Die Gl. (8) geht in die vorher diskutierte (7) über, wenn man o/=0 setzt. Wollte man 
anderseits die Querkraftdeformation vernachlässigen, aber die Rotationsträgheit mitnehmen, 
so hätte man G F,=x zu setzen und würde dann eine Gleichung von genau demselben Bau 


erhalten. In beiden Fällen verschwindet das Glied mit © und so, wie Gl. (7) ohne Rotations- 











XUM 


Z. angew. Math. Mech. ’ i , i . - 
Bd. 22 Nr.6 Dez. 1942 Flügge. Die Ausbreitung von Biegungswellen in Stäben 315 





trägheit Wellen mit der Laufgeschwindigkeit e=cy=y @F,/n liefert, muß die Gleichung 
mit dem Rotationsglied, aber ohne das Querkraftglied zu Wellen mit der Geschwindigkeit 
c=Cy=YEIJol=yYE/o führen, die also die normale Schallgeschwindigkeit haben 

Für einen prismatischen Stab aus einheitlichem Material ist der zum Koeffizienten von so” 
hinzugekommene Anteil o/ kleiner als der schon in Gl. (7) vorhanden gewesene; denn es ist 
u/F,>o, E/G >1, aber er ist doch etwa von derselben Größenordnung, so daß es in diesem 
Falle keinen Sinn hat, nur einen der beiden Einflüsse zu berücksichtigen. 

Die vollständige Differentialgleichung muß beide Wellengeschwindigkeiten enthalten. 
Um uns davon zu überzeugen, benutzen wir nicht die Methode der harmonischen Erregung, 
die den Vorzug hat, gleich ein allerdings mühsames Rechenverfahren zu liefern, sondern 
untersuchen einmal die Vorgänge an einer Unstetigkeitsstelle des Spannungszustandes, die 
mit einer dieser Geschwindigkeiten den Stab durchläuft. 

Bild 5 zeigt ein Balken- 























Xy element d.r, das von einem Wellen- 

MN,” m, kopf mit der Geschwindigkeit e in 
/ 4 \ der Zeit dt=dzx/e von links nach 
0, pP B rechts durchlaufen wird. Vor der 

! B, Welle, also rechts in Bild 5, ver- 

Ang p—% — sehen wir alle Größen mit dem 
70, ro, - Index 1, dahinter mit dem Index 2. 


Bild 5. Bild 6. Jede dieser Größen gibt es nur in 

einem bestimmten Bereich der .r - t- 

Ebene (Bild 6). Beide haben eine Grenzgerade der Neigung dt/dr=1/e gemeinsam, die wir 

als Wellenlauflinie bezeichnen. Den endlichen Zuwachs, z. B. M, — M,. den eine Größe er- 

fährt, wenn die Welle über eine Stelle des Stabes hinwegläuft, bezeichnen wir als Sprung- 

größe. Diese Sprunggrößen gibt es nur auf der Wellenlauflinie. Da an dem Wellenkopf der 

Zusammenhang des Balkens gewahrt werden muß, so müssen w und y dort stetig sein: da- 
gegen dürfen sich deren Ableitungen und die Schnittkräfte M, 0 sprunghaft ändern. 

Da die Sprunggrößen nur auf der Grenzgeraden existieren, können sie auch nur längs 
dieser Geraden in der &-t-Ebene differenziert werden. Dazu brauchen wir noch einen Maß- 
stab, durch den wir den Abstand zweier Nachbarpunkte ausdrücken. Dafür können wir u.a. 
eine der beiden Koordinaten x und f wählen. Wir wollen uns für ? entscheiden und be- 
zeichnen “diese besondere partielle Differentiation in der .-f-Ebene längs der Geraden 


d RR 
x —=et+ const durch ein besonderes Symbol dr Dann ist?) 


! 
d Ö 0 
— & 24) 
dt 0° dx 1“ 
Wendet man diese Operation. auf die Gleichung w, ww 0 an, so folgt: 
| sm, ' 
d , ; ‚ 
v7: In EIER 200 te 5 
und entsprechend gilt 
vr, —-Yv,tew;—y)=0 .... 534 


Wir schreiben nun für das in Bild 5 gezeichnete Element den Impulssatz auf. Er heißt für 
Impulse in Richtung der Durchbiegung 


. dx 
uw, —w)de=(Q, — 9.) dt= (9, (),) 
Br 
und für Drehimpulse a | 5 
d.r 


ol (y,—y,)de=(M,— M)dt=(M, — M,) 
Mit Benutzung des Elastizitätsgesetzes (4) auf der rechten Seite wird daraus 
uc(wn,—- w)=GF,(w — wm, +y,—y)=GF, (m — w)), \ 
ß % - € B TE 
oTe (ip) = Ey; — y;) 


und mit (10a, b): 
uc(w — w)=GF,(w, —w)), 


ol (pi y)=EI (pw). 


2, Vgl. die entsprechenden Überlegungen bei W. Kucharski, Z. angew. Math. Mech. Bd. 21 (1041), S, 152, 
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Diese Gleichungen haben zwei nichttriviale Lösungen. Die eine heißt 


’ / 0 


c=@g=y6GF,u, ww, w=UVi - -» - : . . . . (18a) 
und die andere 

c=cy=yYEljel, w&£y, WW ..: .. 0.2.05, .fl18b. 
Ein Sprung in w’ breitet sich also mit der Schubwellengeschwindigkeit ec, aus, und dabei ist 
y" stetig, also auch M, während der Querkraftsprung 


Bi. ea ae 


ist. Anderseits pflanzt sich mit der Biegewellengeschwindigkeit ey ein Sprung in y” fort, 
und dabei sind w’ und @ stetig, während der Momentensprung 


EEE ) 2 m Sul”, > 


ist. Ein Sprung im Biegemoment und einer in der Querkraft können also nicht gemeinsam 
bestehen, sondern laufen verschieden schnell den Balken entlang. 

Für jede der beiden Unstetigkeiten haben wir an der Sprungstelle einen Satz von z. T. 
unstetigen Übergangsbedingungen, die in diesem Punkt die Differentialgleichung ersetzen, für 
den @Querkraftsprung z. B. die Aussagen über die Stetiekeit von w, y und M und die Be- 
ziehung (14a) zwischen den Unstetigkeiten von mw’ und @. 


5. Entwicklung der Sprunghöhe in der laufenden Welle. 

Es bleibt noch die Frage offen, wie ein irgendwo vorhandener Sprung beim Weiter- 
wandern mit der Wellengeschwindigkeit e seine Sprunghöhe ändert. Wir wollen uns diese 
Frage für die Querkraftwelle überlegen. Für die Momentenwelle verläuft die Rechnung ebenso 
und führt zu einem ganz entsprechenden Ergebnis. 

Der Impulssatz ergab nach (11a) für die Querkraftwelle die Beziehung 


(), A ucolmw, ",). 


2 “1 


Differenzieren wir das längs der Wellenlauflinie, so erhalten wir nach (9). 


51 (,-Q)=zuon tom — Wh, —cCoW)). 
Nach (6a) und (i3a) ist in beiden Bereichen i 
iv = ey (mw u") n 
also 
d P) 7 „ 2 ’ ’ .r . 
57, 09, — 9) = ucol[c (m w;)— colyı —y;) + co (wi — w;)]. 


Da y” auf der Querkraftwelle stetig ist, so ist rechts das mittlere Glied null, und die beiden 
anderen lassen sich wieder zu einer Ableitung der Sprunggröße w; — w; längs der Wellen- 
lauflinie zusammenfassen: 





d 


3,9%): (w, — my). 


uct 
 ı 
Anderseits hatten wir nach dem Elastizitätsgesetz die Beziehung (14a) zwischen dem Sprung 


von Q) und von w’ gefunden, die differenziert heißt: 


Y r dD ’ ’ 
-Q)=6R,( wm —m). 


Durch Vergleich dieser beiden Gleichungen folgt 
Li 
JA 


Die Höhe des Sprunges verändert sich also nicht. 


9,.-0Q)=0, 9. 


> 


Q, = const. 


Wir haben damit zwei wichtige Feststellungen gewonnen: 


1. Auch im gebogenen Balken gibt es Kopfwellen, die mit konstanter Geschwindigkeit 
laufen, und zwar deren zwei. Die eine rührt her von der Wechselwirkung zwischen 
Querkraftsteifigkeit und Translationsträgheit, die andere von der Wechselwirkung 
zwischen Biegesteifigkeit und Rotationsträgheit. 

2. Die Höhe eines einmal eingeleiteten Sprunges der Querkraft oder des Biegemoments 
ändert sich beim Weiterlaufen längs des Balkens nicht. 
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Was hinter den Kopfwellen geschieht, ist damit nicht gesagt und läßt sich auch nicht 
so einfach ausrechnen und darstellen wie beim Längsstoß und seinen Analogen. Die mecha- 
nische Anschauung lehrt unmittelbar, daß sogar recht komplizierte Vorgänge eintreten müssen. 
Leitet man z. B. von t=0 an ein konstantes Biegemoment M (0,f)=M am Stabende ein, so 
pflanzt sich die Kopfwelle der Höhe M mit der Geschwindigkeit ey fort und erteilt jedem 
Querschnitt, über den sie hinweggeht, nach dem Impulssatz die Drehgeschwindigkeit 

M 

oO I cy j 

Da zunächst die Querkraft null ist, gehört dazu nach (4b) ein ebenso großes w’, und es muß 
also hinter der Welle eine Querbewegung anlaufen. Um sie in Gang zu bringen, sind aber 
Querkräfte nötig, und wo diese auftreten, muß das Moment sich ändern. Das Ergebnis muß 
schließlich ähnlich aussehen, wie die Lösung (3) der vereinfachten Theorie, natürlich mit 
den Änderungen, die sich aus dem Bestehen der Sprungwellen ergeben. Wie man sieht, 
wird der Vorgang recht verwickelt, und seine Berechnung erfordert letzten Endes die Auf- 
lösung der Diiferentialgleichung (8). Diese Aufgabe soll nicht Gegenstand dieser Untersuchung 
sein, doch wollen wir als wesentliche Voraussetzungen dafür noch einen Blick auf die zu 
Gl. (8) gehörenden Randbedingungen werfen. 


6. Randbedingungen. 
Die Differentialgleichung (8) ist in w» von 4. Ordnung. Wir werden also an jedem Stab- 
ende für w und seine Ableitungen zu allen Zeiten f zwei Bedingungen vorschreiben müssen. 
Am einfachsten ist das, wenn am Anfang des einseitig unendlich langen Stabes ıw und 
M als Funktionen der Zeit vorgeschrieben sind, wenn also z. B. der Punkt gelenkig gestützt 
ist und ein Moment plötzlich aufgebracht wird, oder wenn er einspannungsfrei gehalten und 
plötzlich mit gegebener Geschwindigkeit senkrecht zum Stab ausgelenkt wird. In diesen 
Fällen hat man sofort einen vollständigen Satz von Randbedingungen für die Differential- 
gleichung (8); denn das Moment läßt sich nach (6a) durch Ableitungen von w allein ausdrücken: 


[A 


M=-—EIy'= Ein" +Bl2 #6 .....% 0.3. MB 


GF, 
Hat man es mit Randbedingungen zu tun, die yw selbst enthalten, so braucht man 
eine Beziehung zwischen y und ı, die keine Ortsableitungen von y enthält. Eine solche 
Gleichung läßt sich aus (4a,b), (5a), (5b’) gewinnen. Aus (4b) und (5b’) folgt nämlich 
GF,(w —y)=oIy+ M' 
und durch Einsetzen von M aus (15): 
ol EI 


r 


„rare“ + Gr") ha N ee IE Ss VE 


Ist nun die Endverdrehung y (0, f)=f(t) vorgeschrieben, so ist natürlich auch % (0, f) 
f (t) bekannt, und wir können aus Gl. (16) entnehmen, wie groß deren rechte Seite für 

x2=( und alle t ist. Damit haben wir eine allerdings schon etwas komplizierte Randbedingung 
für w gefunden. 

Ähnlich läßt es sich auch machen, wenn die Querkraft @ (0, t) = g(f) gegeben ist, in- 
dem man mit dem Elastizitätsgesetz (4b) 
3) 
GF,’ 
durch die gegebene Funktion g (f) ausdrückt und links in Gl. (16) einsetzt. Geordnet ergibt 
sich dann 


y (V, t) == mw’ Mm) = 


u, / u RR ” 
Im+gFIM=-EllarntEr)® BETE ra  e re 
und das ist wieder eine Randbedingung für w. 

Haben wir danach w (x, f) bestimmt, so können wir aus (6a) y durch eine einmalige 
Integration nach x berechnen, wobei sich die Integrationskonstante aus der Forderung er- 
gibt, daß in dem von der Störung noch nicht erfaßten Teil y = 0 ist. 


7. Ergänzungen und Einschränkungen. 

Die Differentialgleichung der Balkenbiegung, gleichgültig ob mit oder olıne Querkraft- 
glied, gilt nicht nur für prismatische, homogene Stäbe, sondern für alle schlanken Stäbe, 
mögen sie nun inhomogen sein oder irgendeine Feinstruktur besitzen, wie etwa ein schlanker, 
regelmäßig gebauter Fachwerkträger. Entsprechend lassen sich natürlich auch die dyna- 
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mischen Differentialgleichungen (2), (7), (8) auf alle diese Stäbe anwenden und auch auf solche, 
deren Trägheitsglieder von zusätzlichen Massen herrühren, mögen diese nun kontinuierlich 
über die Stablänge verteilt sein oder in Form vieler kleiner Einzelmassen. 

In Rücksicht auf solche Fälle haben wir in dem Quotienten EI/oI in Gl. (13b) nicht 
(die beiden Trägheitsmomente gegeneinander gekürzt, denn sie können ja etwas ganz ver- 
schiedenes bedeuten und dann auch verschiedene Werte annehmen. Im Zähler steht die 
Biegesteifigkeit. Hier ist also bei der Bildung von I nur über diejenigen Teile des (Quer- 
schnitts zu integrieren, die an der Aufnahme der Biegespannungen beteiligt sind, und wenn 
darunter Teile mit einem anderen Elastizitätsmodul vorkommen, wie z. B. die Eiseneinlagen 
eines Eisenbetonquerschnitts, so sind sie bei der Bildung des Trägheitsmoments ihrem größeren 
oder kleineren Elastizitätsmodul entsprechend höher oder geringer zu bewerten. Im Nenner 
dagegen handelt es sich um ein wirkliches dynamisches Trägheitsmoment, in das alles ein- 
zeht, was Masse hat. Hier sind also nicht nur die tragenden Querschnittsteile in Rechnung 
zu stellen, sondern alles, was mit dem Stab fest verbunden ist und an der Drehung y 
(uerschnitts teilnehmen muß, wie z. B. die schraffierten 
Zusatzmassen in Bild 7, und wenn der Stab aus verschie- 
denen Werkstoffen zusammengesetzt ist, so sind die Ein- SSRASNSNN 
lagen jetzt nicht nach ihrem Elastizitätsmodul, sondern nach l Er 
ihrer Massendichte o zu bewerten. N N N N N N N N 

Im Gegensatz zum einfachen prismatischen Stab, wo un 
das Verhältnis ey:co nur in engen Grenzen veränderlich 
ist, kann man also im allgemeinen Fall beide Geschwindig- 
keiten unabhängig voneinander sehr stark beeinflussen, cy durch Zusatzmassen, wie sie Bild 7 
zeigt, cg durch Schwächung des Schubverbandes, so daß man es in der Hand hat, einen Stab 
zu bauen, in dem die Momentenwelle oder die Querkraftwelle schneller läuft. 

Selbstverständlich sind der Anwendung der hier entwickelten Gleichungen auf die wirk- 
lichen Vorgänge auch Grenzen gesetzt. Für den einfachen prismatischen Stab wird man nicht 
erwarten können, daß Wellen, deren Wellenlänge von der Größenordnung der Stabdicke ist 
oder noch darunter liegt, noch richtig wiedergegeben werden, sondern dann ist natürlich eine 
Behandlung des Stabes als 2- oder 3-dimensionales Kontinuum nötig. Ebenso wird man bei 
Stäben mit Zusatzmassen oder mit irgendeiner Feinstruktur des elastischen Aufbaues nicht 
erwarten dürfen, daß Wellen, deren Länge den Abstand der Strukturelemente erreicht oder 
unterschreitet, oder deren Frequenz gleich der höchsten Eigenfrequenz der mechanischen 
Kette ist (bei der benachbarte Massen paarweis gegeneinander schwingen), von einer Konti- 
nuumstheorie, wie wir sie aufgestellt haben, noch richtig erfaßt werden. Hier hat eine Be- 
handlung des Gebildes als ein mechanischer Kettenleiter einzusetzen. t14 


des 


Bild 7. 


Halbraum mit Halbkugelbelastung'. 
Von Constantin Weber in Dresden. 


Verfasser zeigt, daß die in der Arbeit von Usunoff’) vorkommenden Integrale 
durch Umwandlung des rotationssymmetrischen in ein zweidimensionales 
Problem in geschlossener Form ausgewertet werden können. 


Da die Lösung für den Halbraum mit einer Einzellast bekannt ist, so kann allgemein 
die Lösung des Halbraumes mit beliebiger Belastung in Integralform angegeben werden’). 
Bei rotationssymmetrischer Belastung wird man hierbei erst nach dem Umfang, dann naclı 
dem Belastungshalbmesser integrieren. Schon die erste Integration führt auf nicht elementare 
Integrale. 

Bei der Halbkugelbelastung empfiehlt sich eine von mir vorgeschlagene Methode’). 
Hierbei wird durch den Translationsprozeß T aus dem rotationssymmetrischen Fall ein zwei- 
dimensionaler gebildet, aus diesem durch die Operation D ein neuer zweidimensionaler und 
aus dem letzten durch den Rotationsprozeß A wieder der ursprüngliche rotationssymmetrische 
Fall. Es zeigt sich, daß die Spannungen o, (parallel zur Rotationsachse) und die Spannungen 
Tx„ elementare Funktionen geben. 

1) Nach Drucklegung fand ich den Aufsatz von M. T. Hnber: „Zur Theorie der Berührung fester elastischer 
Körper“, Ann. Phys. Bd. 14 (1004), S. 153, in dem alle Spannungen in geschlossener Form gefunden sind. 


>, N. Usunoff: Über den Spannungszustand im Halbraum bei halbkugelförmiger Druckverteilung. Z. angew. 
Math. Mech. Bd. 22 (1942), S. 262 bis 269. 


‚G. Lundberg: Elastische Berührung zweier Halbräume. Forsch. Ing.-Wes. Bd. 10 (1939), S. 201. 

4,0). Weber: Achsensymmetrische Deformation von Umdrehungskörpern. Z.angew. Math. Mech. Bd. 15 (19:5), 
S. 416 bis 408. C. Weber: Zur Umwandlung von rotations-symmetrischen Problemen in zweidimensionale und 
umgekehrt. Z. angew. Math. Mech. Bd. 20 (1940), S. 117 und 118. 
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Für den oberen Halbraum y=0 sei nach Bild I die Randbelastung 


Oy: } 1 Y” 
Prozeß T, der darin besteht, daß wir die Gesamtbelastung, die in jedem 
auf die Tiefe 2 verteilen, erhalten wir das 


Durch den 
Schnitte senkrecht zur Bildebene vorhanden ist, 
zweidimensionale ebene Verzerrungssvstem mit der Randbelastung (Bild 2) 


{ (1 E» 
Die Potentialfunktion mit diesen Randwerten und den Werten Null für 2 © lautet 
7-1 — ge’) 1 E=3_g, 
. | z’yin 2 +1 r 


Das letzte Glied der eckigen Klammer beseitigt das Unendlichwerden für |z 


I.r 





Bild 1. Halbraum 3ild 3. 
Zweite Belastung der Halbebene. 


mit Halbkugelbelastung. 


Erste Belastung der Halbebene. 


Durch die Operation D bilden wir eine neue Belastung der Halbebene. Wir differen- 


tieren hierzu die Potentialfunktion nach y und ändern das Vorzeichen; wir erhalten 
l 2—1 
y R|Z en, +1): 





Die Werte für y=0 nehmen wir als Randbelastung (Bild 3) 


e | 
wg z.lln(e - D-In@e+d)+1. 


(Der obere Index (2) gibt im weiteren an, daß es sieh um den zwei-dimensionalen Spannungs- 


zustand handelt.) 


Der Spannungszustand dieses Belastungsfalles ist mit 9 (#,y) gegeben: 


” Ö 4 
O4 (x, = q „ 
: I, y 
0g 
7 De y)= Y U 
xy I . O r 
Ö 4 


0,2 z,y)=y+y Iy 
’ ( 


0,7 (a,y)=ulo” +0," )=2up. 


Dieser zweidimensionale Spannungszustand hat die Eigenschaft, daß er durch den 
Rotationsprozeß R in den ursprünglichen rotationssymmetrischen Spannungszustand übergeht. 
Hierbei denken wir nur, daß die senkrecht zur Bildebene unendlich ausgedehnte Belastung 

" R u 
bei der Rotation von «=0 bis a=7 auf der Strecke da den Teil der Belastung auf 
A 
die Oberfläche des Halbraumes abgibt. 


met, Bu 


Hiermit folgt durch nur eine Integration aus den Spannungen ©, 
(die Spannungen des rotationssyimmetrischen Problemes: 





XUM - 
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A 
: al.) 
ah \o,” (rcosa,y) da = \o,” (x, y) 


i er u? [*) 


1 3 2 Fin al r 
Tı«y >= j- \ Tx (tr cos a,y)cosa da= - \ Tx u (2,9) . Br 
m 1) | 1 : 
(=) 
x 
' 7 ' ' a! s 
OR & Lg ı 4 j 
rn \ (0x'” cos’a + 0,” sin’a)da 7 \ (0x a er r 
i i Vı-(-) 
r 
ı() 
4 ( 
1 . , R 2 ) Ad r r 
e7) r- \ N sın’ a +0," cos? a) da = > \ lo, 2 r 06; Re a ö 
0 ; v | | | 
Y 


Für die Rotationsachse gibt die Integration nach a unmittelbar: 
Fr 


a, 3 —d, 0, Re rain, oa: Fa) 


Für die Halbkugelbelastung wollen wir noch die Berechnung von 6, durchführen. Es wird: 


= ! 1 z+1 % @/14 2 1 
u," R| > z In FRE Ya I+yK FA 5 zIn or 
. 1 2—1 l 2 —1 9 E 2 
> | 2 Be a * L 
u | Baia" 7 55 Due € SET Bar FE a Are 
3 iy/ 1 1 
Rig-h; +1 5 |; ira 
Hieraus: 
2 | al“) 
2 1 2 iy/ 1 1 j 
os ala Hit el -itzr 
0 | l | } 
Durch partielle Integration von = - erhält man: 
Vı-() 
*) . = .\2 . 
u ta > | e. y | A FERIEN A | 
Gy r- h 73 In, 1 1—| — | + aresin — | 





Das erste Glied der geschweiften Klammer gibt den Realteil Null. 


XUM 








XUM 
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Für die Integrale erhält man 


d\—)= (: arcsin - 


= z (7) (5 + aresin -). 
Adır | 








A a | a’ 7 l 
N ’ 
Hiernach wird 
) . s ; 2 . ; 
2 „8, #1 ‚fiy I\r _ iy I D 7 na ' 
" BL R\3 2 | +, r)2 ( r "ya iy) | 2 iu 
r u u .-1% r n 22,5 r 
pH rs ı 77 -))ViFiyp (3 N iy) 
” - (> + arcesin “ 4 = aresin )\ 
2r| yvüy-W—r\2 7 TO 1—-iy Vüy+W—r\2 TUT MLHiy)f' 


Die unterstriehenen Glieder heben sich fort. Den Restausdruck schreiben wir nochmal hin, 
wobei wir die Vorzeichen vor allen © wechseln. Man erhält dann dasselbe o,. Bilden wir 
die halbe Summe, so fallen alle Glieder mit aresin fort: wir erhalten hierbei sofort den 
Realteil: 


| rY I 1 Hr) | Y Y 


“v 2 r "r)yäiy- Wr" vüiy+ly— ri] 


iy| r r N 
Ir |\yüiy—- Wr" Vüy+W—r|f' 


Wir bringen die Brüche auf den Generalnenner und zerlegen die Wurzeln der komplexen 
Ausdrücke nach der Formel: 


ı 1 M R I 1 
Vatib | + HP FIR HH V - Ha FR. 


Hiermit wird: 


1 A R 1 - 5. 3 
ai! - ı) +  Ü-y„-rfriyf 


(1—-?—r?) 


1 1 
+y -;(i1-y„—-r) + yil-y—r)+4y 
] Mi - y r’+ ty? 


In ähnlicher Weise findet man 


1 | 
FrH)VY su-F -N+> ua-y- rip 


| 1 
Y sit —-„-r)+zyi-y„—r’ +4y 
Try 9 = 
ı ) (1 y = y})? + 4y 
421 


21 
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Die Spannungsspitzen in gelochten Blechscheiben 
und Streifen. 


Von A. Hütter in Dresden. 

Verfasser sucht zahlenmäßig den Einfluß von Lochungen auf die Spannungen 
in auf Zug beanspruchten Scheiben und Streifen zu erfassen. Im ersten 
Teil wird die unendliche Scheibe mit einer Lochreihe untersucht. Für 
kleine Lochdurchmesser wird die Airysche Spannungsfunktion durch 
Reihenentwicklungen angegeben, während für große Lochdurchmesser zwei 
verwandte Probleme als Näherungen herangezogen werden. Der zweite Teil 
behandelt den gezogenen Streifen mit einem Loch. Die Spannungsfunktion 
wird für kleine Lochdurchmesser wieder mittels Reihenentwicklungen be- 
stimmt. Für große Lochdurchmesser wird die Halbscheibe mit einem Loch 
in der Nähe des Randes als Näherungslösung behandelt. 


KHTTEEITBSTITTITEL DIT TTTITE I I 


ty 
TEE 


il» -. 





Behandelt werden Belastungsfälle nach Bild 1, 2 und 3. 


Die wichtigsten Bezeichnungen: 


2 
iQ 


FTTITITEI TA 


21= Lochabstand 





Derprrennterentnnnn 


: AMEBAAALAABBSLBBBBLLBABLDAL BR ALL AL 
2a = Lochdurchmesser 
PBIEZLILTLLLLITELLLELELLELTLELLLLTDG h 


[A 





2b = Stegbreite zwischen den Löchern 














4, 
: . r i I/ 
x,y = rechtwinklige Koordinaten | 
J ei 60660000000 
r,p —= Polarkoordinaten  " 
z=r+i Y [ERSRERASRAHRRRRDDBALZDLLLERBALLTE. 
e—Rr ntialf ti Zc 
PD —= Potentialfunktion ee — 
2c= Normalspannung im Unendlichen 
'= Spannungsfunktion | 
Ga, — maximale Spannung am Lochrande ER. 
ar Hr u 
54 — durehsehnittliche Spannung im Steg. Bild 1 bis 3. 





I. Die gezogene Blechscheibe mit einer Lochreihe. 
1. Der Lochdurchmesser ist klein im Verhältnis zum Lochabstand. 


a) Die Spannungsfunktion. Für das ungelochte Blech lautet die Spannungs- 


funktion: 
ee (0 +) bei allseitigem Zuge . . . ». .». 2:2...) 
F=c.? bei Zug in der y-Richtung . . . . 2... 
Zu diesen Lösungen nehmen wir periodische Glieder hinzu, die in den Punkten z=0, +, 
+27... bestimmte Singularitäten aufweisen; 21 ist hier gleich z gesetzt. Beim Bleche mit 
einem Loch wird bei allseitigem Zuge das Glied 2ca?-AInz hinzugefügt, bei Zug in der 


. i : Eu 1 1 i 
y-Richtung die Glieder ca’ -Rlnz+ zZ ca -RG—5:ca?:-y-R—. Ganz entsprechend 


2 
- 


verwenden wir hier periodische Glieder ®, und y®,. Wir machen die Ansätze 


F=c@+y)+P,+yP, bei allseitigem Zuge . . . 2..2...ß8), 
F=c’+®P,+yb, bei Zug in der y-Riehtung . . . . . (4). 
Hierbei nehmen wir: 
: 2 0° : 
P,=a,Rlnsinz+a,R N „insinz+a,R dp ES ee ° 
n | Ö o° | 
yP,=y|b,R Insinz+5,R, ;Insinz-+ ---) (6) 
yY#, 919, oy Hd; oy FT Er I )). 
Die noch wählbaren Beiwerte a,, @,,..., b,,b,... sind für die Lösungen nach (3) und (4) ver- 


schieden. Die Funktionsglieder sind periodisch und symmetrisch zur x- und y- Achse. 


XUM 








XUM 
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’ 


Am Lochrande muß für r=a sowohl F=const als auch —— gleich Null werden. Da 


£ Y 


a a E ; Be oF. 
] 2 als klein angenommen ist, so entwickeln wir F und —— in Reihe 
. ( 


2a 


das Verhältnis 
nach Potenzen von 


Die Reihenentwicklung von a,: Rn sin z lautet: 


/ 1 1 
-RAinsinz=a,NR|lnz -. 2... 
a,-Alnsin oR| 6 180 
Kr (T). 
a (in r r? -Ccos 2« r’-cos4« ae ) 
”r a den 
Für die übrigen Glieder erhalten wir die Reihen durch entsprechende Differentiation: 
e,. ur 1 ie an f 
u«,AsZ;insinz= a, 7 4, |; C0S&c = 7” C082« ar c0Ss%«( 2 7 
od :""sin?z ij 75 777189 + 
a,R 0" Insina= — 2a R| L - 2 | 
a - 1 |sin’ze sinie] 
(9), 
8 | 
=—_D9 | 08 — nnd rin 
=— 2a cos4@ — 5; r’cos2« -r'cos4 . 
a Me E Ts 4 | | 
4 MR Ei EEE ER Burg | 
b,R j,Insinz=b, Rietgz=b, [sin? q +yzr’sin2p+ jr" sing sin39-+:-::| . . . .. (10), 
0° & , cos2z 
b,R;— Insins=26, Ri -—; 
70 sin’z 
(11). 
| | 
-91 > si sin 8 Re "sin? Agi gi By Er. 
2b, 5 ing sin3p+7z r"sin’p+ zgg r'sinpsindp+ | 


IF 
R 0 2 n c . de 
Diesen Ausdruck setzen wir in den Ansätzen von F und __ ein und verlangen, daß möglichst 
f - 
) 


hohe Potenzen von verschwinden. Damit ergeben sich die Konstanten für allseitigen Zug: 





a=—2ca? year; 8 
’ 6.= 3 ca* & ca°: 
d, 3 ca +7zea®; * 
N BD: Ze. 
am 2 ie 


Mit diesen lautet die Spannungsfunktion: 


B j BAD 
F=cr—c ((e +7 a°| R In sin z +| ya+7 a‘) N 














sin?z 
. 2 ) (12). 
’ & € ß ı fe) 4 cos2 
A ne 3 En 2 E a) Minto Be er 7 2\ 
gg R sin?z Fr] +| 3° r 9 yRicige 5"9 Ri sin?z] 
In gleicher Weise erhalten wir für Zug in der y-Richtung: 
9} 5} 
ze Sn = 
= —R8 yeca r 9.ca 9 59 
b.=9 u gt 
f f B; ca zeatea; 
— a 4 a‘ 6 z 8. 
d,= > ca’ + 6 ca 3964 . f 
B= ce. 
FEN 10, 
wegen"; 
’ np? € I(o 2 4 4 4 IR 31 ( 4 13 6 11 815 1 
F=cx— 3 (2« + za ga Rlnsinz+ a’ + za jz« Renz 
1 > 3 4 104 un. 
2 € 0 4 cos2] 
ART 109 . Ba 4a? u % 6 icte2 — — a yRi— . 
30° sin’z sin’z + er 5 2° \yRicte 5" u. sin?z/ 
21° 
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Durch Reihenentwicklung kann man sich leicht davon überzeugen, daß F bis zu den Gliedern 


>F 
vierten Grades einschließlich und _ _ bis zu den Gliedern dritten Grades die Randbedingungen 
on 
erfüllen. 
b) Die Spannungserhöhung. Aus den oben ermittelten Spannungsfunktionen er- 
halten wir die maximale Spannung am Lochrande mit r=a und g=%: 


’ Ss @ 16 
für allseitigen Zug: 1 FE 7 e| + a'+:-:) (14), 
EEE 0. £ 
für Zug in der y-Richtung: „= 0r= el6+ za+ “A (15). 


Ist die Lochentfernung nicht gleich x sondern 21, so ist jeweils für a der Ausdruck 
na 
37 zu setzen. 
Die durchsehnittliche Spannung im Steg beträgt: 


1 




































a=2e (16). 
1 er 
I 
Dies liefert die Spannungserhöhung 
f 6) 2, 2 4 4 ] 
2 a\| 2 /n a ı /n a I ze sr h u 
an 2 [1 — -\1+53(5) | -| >) . +:.-| für allseitigen Zug . . (17), 
a I)» s\ı83/\3 52a | - 
f \ ” 4 4 1 
r, a 2 ai’. I A : . 
max — 3-[1 ) l+- | = (—) +...) für Zug in der y-Richtung . . . (18) 
Od l ' a I | B 
2 u © ” | ÜTLLILITITSLILTEITIELTERTELETESTT 
78 3l 00000000000 |5 z6\ se | 
3 E | a 
3 h = | | 
x 6 at T IPTERTTPPRTTTITITTERTRTTTETERRTTN) 
S IS] S_ I 
sr Se | 
PA I nn 2 
| 
Zu a 
I 
| 
0 0 00 08 04 09 an —%0 TE; “0 0 4 4 
WW 08 0m 08 075 7 00 05 02 YYl— 0008 } ”% 7 r r 
Bild 4. Bild 5. 


Diese Ausdrücke geben für a/l>0 den Wert der Spannungserhöhung und den 


7: 
weiteren Verlauf bis zu den Gliedern fünften Grades einschließlich. In Bild 4 und 5 sind 
diese Schaulinien dargestellt. 


2. Die Stegbreite ist klein im Verhältnis zum Lochabstand. 


s R a ’ Ro 
Ist das Verhältnis von Lochradius zur Stegbreite groß, also ] <1, so kann man die 


Spannungserhöhung auf eine andere Weise bestimmen. H. Neuber') hat mittels angepaßten 
Koordinaten für ein Blech mit Hyperbelkerben nach Bild 6 eine Lösung gefunden. 


Ist der Steg im Verhältnis zum Krümmungsradius klein, kann man dieses Ergebnis ohne 
weiteres auf den vorliegenden Fall des Bleches mit periodischer Bohrung anwenden. Dabei 
spielt es keine Rolle, ob das Blech allseitig oder nur in der y-Richtung gezogen ist. Der 
von Neuber gefundene Wert für die Spannungserhöhung ist: 


MN l | b 
max __ ji I b l b 


f Od I a b 7 bh 
1 „arctg | re Ib 


(19). 


1) H. Neuber: „Kerbspannungslehre‘“. Berlin 1937. 
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Dieser Ausdruck läßt sich leicht umformen und in Reihe nach entwickeln. Man erhält: 













































Om _4,2,5,8 (br 80 (by 
- — oe © = ER Br ”) 
: Fa =*| | +. An Dre a + ME 
04 3 ı "urn 
We) 2 
[ LELTTITTITETTETETTV TOTEN 
N - > 
a > E CD) (A E 
nach Poschl = 58) z 
16 01.23 = _ 2l - = 
SS ei = 
a; “ 7% “inmmammn ÄMUEABABALARAEBLDBLGABABDRLARDRD BEE] 
3 : 
E nn EB Bild 8. 
7 E OO 3 d 
|p f = = 
' > It EITIPITTeTeTTT £ 
Bild 6. -00 97 02 03 04 5 06 97 08 ) 10 Bild 7 (links). 
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Das Ergebnis dieser Lösung ist in Bild 7 aufgezeichnet. Für große Lochdurchmesser 
kann man auch eine von Th. Pöschl?) angegebene Lösung für das unendliche Blech mit 
zwei kreisförmigen Löchern nach Bild S benützen. Dann lautet die Spannungsfunktion: 


u . : ah R 
e+y7+—:-y[lIn<a?+(y+e?>—In<2?-+(ı e)?’ > \ a 
| pr Yy | Yy N 


Hierbei ist ®=(2a+b)b. 

Der Verschiebungszustand dieser Lösung ist mehrdeutig. Er entsteht dadurch, daß für 
positive & ein schmaler Streifen aus dem Steg herausgeschnitten und die Scheibe gewaltsam 
wieder verbunden wird. Es handelt sich also um eine Vorspannung der Scheibe. 

Für die Punkte der y-Achse erhalten wir: 


BER; 
| Re 2 = Iy+2 , In(y-+e) - In(y 9) EEE 
und 
N ta? e? ” 
0x = ee | Kae A ar AD). 
X,‘ 0 Pi (e? y°)) (29) 
Für y=+b ergibt sich als maximale Spannung 
? 53 
Omas = (x, y=d =& 7 De : ae ee re A 
Durch Integration von o, im Bereich y= —b bis y=+b erhalten wir die gesamte im 
Steg übertragene Kraft und daraus die durchschnittliche Spannung 
I 2 Ja?h] 
p| ‚0 sad8| en 
um |b- IIn(e+b) — In(e— b)}-+ — ’ 25 
Od 7 1? rZi ( ) ( )y Pr p2) ir : £ 6 


l 


ir £ a n Mn i i ’ . : 
Den Ausdruck für die Spannungserhöhung -""* entwickeln wir wieder in Reihe nach 
b/l und erhalten 
Om a 5b u /85%. 99 b\’ 
ae Br a a ER _ ' 
Od 3 l 15 ® 3780 l 


Die Näherungen, die aus den Lösungen von Neuber und Pöschl gewonnen werden, 
stimmen also einschließlich der linearen Glieder überein. Sie ergeben somit genaue Werte 
für die Spannungserhöhung im Punkte b/1>0 und für die zugehörige Tangente im Schau- 
bilde. Die Krümmung dagegen ist unsicher. Bild 7 ergibt einen Vergleich der beiden Lösungen 
von Neuber und Pöschl. 

3. Ergebnis. 

Aus der Reihenentwicklung für kleine Werte a/l und aus Endwert und Tangenten- 

richtung für kleine Werte b/I=1--a/l kann man mit ziemlicher Sicherheit den Verlauf der 


I max 


Schaulinie für die Spannungserhöhung angeben. 


d 
Bild 9 zeigt das Ergebnis für allseitigen Zug, Bild 10 für Zug in der y-Richtung. 


2) Pöschl: Z. angew. Math. Mech. Bd. 1 (1921), S. 174. 
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Bild 9. Bild 1. 


II. Der Streifen mit einem einzelnen Loch. 


1. Der Lochdurchmesser ist klein im Verhältnis zur Streifenbreite. 

Wir gehen von einem Ansatze nach Gl. (4) wie bei der in y-Richtung gezogenen Ebene 
aus. Weiter bilden wir uns die Glieder der Funktionen ®, und yP,, die im Koordinaten- 
anfangspunkte wieder verschiedene Singularitäten aufweisen, am Rande des Streifens aber 
die Spannung Null geben. Wir nehmen hierzu die entsprechenden Glieder des in y- Richtung 
gezogenen Bleches und schneiden einen Streifen heraus (Bild 3). Der Rand ist nur durch 
Normalspannungen belastet, da die Schubspannungen wegen der Symmetrieeigenschaften des 
unendlichen Bleches verschwinden. Die Normalspannungen beseitigen wir dadurch, daß wir 
dem ungelochten Streifen einen geeigneten Spannungszustand überlagern. Zu diesem Zwecke 
untersuchen wir zunächst den ungelochten unendlich langen Streifen. 

Die Glieder ®, und yP, des unendlichen Bleches liefern eine Airysche Spannungs- 
funktion, die durch folgende Glieder und Spannungen gekennzeichnet ist: 








a } , 1 
F*—#RInsinz, = 
sin?z 
u | .- [| 4 6 
Fra -—. "_R|— et, 
sin?’z sin’z  sin?z| 
TEE .@: *E%* IR >: cos? 
Fr_ yRictige, er = + Ri; 
‘ sin’z sin’z 
it d g rate N l 
Am Streifenrande erhalten wir mit e=+l!l=+ — die Normalspannungen p} = Col’y’ 
E yrgy #*E% Fü liese Randbelas N oll ir le Sn: os 
pP : Sof? y u: Cory‘ ür diese Randbelastungen wollen wir nun den Spannungs- 


zustand eines ungelochten Streifens bestimmen und benutzen die Methode der Darstellung der 
Randbelastung durch Fourier- Integrale’). 
Für die Spannungsfunktion machen wir den Ansatz: 
= 


Ft [A Cojix + Bir Sind x] cosiyd/ 
) 


7). 


Die Konstanten A und B bestimmen wir so, daß auf dem Rande die Normalspannung 
den vorgeschriebenen Wert annimmt und die Schubspannung Null wird. Es ergeben sich 


2(&inAltAlCofil)2 \ ’ 
A= » (&insii teil 7 \p@eosätat, 
0 
2 Sinil 2 2 h 
B= 5 &imaaı rain 7 \pMecositdt. 


1 . 


0 
p (t) bedeutet hierbei jeweils die vorgeschriebene Randbelastung 


| reg ie PO 


| 3) Vgl. Seewaldt: Die Spannungen und Formänderungen von Balken mit rechteekigem Querschnitt. Berlin 19237 
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Damit lautet die Spannungsfunktion 


v. 


IV 


uch nuhei ab ER AN 2 I. 
Sinil ride [nenärarlas 


m | 2 \ Gojircosiy 


I 


= 2Sinil I\ »(neosätdtlai 
(Sin2/il+2ID2# \)*‘ 54 di 


0 v0 


r \ ix ZSiniwcosiy 





% 


Der Wert der Integrale \p(f)-cos4t-dt läßt sich für die gegebenen Randbelastungen 


0 
nicht in geschlossener Form angeben. Die Auswertung erfolgt deshalb mit Hilfe der Simp- 
sonschen Regel. Dies ist ohne weiteres möglich, da der Wert der Integranden für große t 
sehr rasch gegen O geht. 
Um den Rand des als klein angenommenen Loches zu untersuchen, wird F in Reihe 
nach r und y entwickelt. Wir wählen folgende Form, in der sich die Konstanten am ein- 
fachsten bestimmen lassen. 


F=a,+a,:r?-cos2p+a,:r’-cos4p+::- . 

ER EA: 
+r.cospl[a,- r-cosp+ta,-r?-cosd3p ++: : | 
Mit den Reihenentwicklungen 
Gol/ixzcos/iy=1+ 5 r’cos2yg * 'costp+ root .... (30), 
- i- 

ee ! er we ! m Er i 

2 Smn/ixcosiy=4#’r’cos2p + g 7" c0S p cos 397 + 120’ cospcosdöp+t:-: . (31 


lauten die Konstanten «a,, a, usw. wie folgt: 


2 7 92 (Einil+AlCojil { 6% Do. 
a, = 7 (Sinailr2iD® \p(l eositdtidA, 
2° 2(Sin äl+AlCofAl) 2 | en 
a..=-— = \ (Sindil+2iDE 3 \p{Weositdt; di, 
u 0 
= f o : 
2 7 12 (SinZl+A1lCojiD 4 
.. 7 \ \ (Sin2i1 +24) 3 \r@Weosärar) ”. 
0 1) 
ce os 
2 7 g2(Zinil+AlCojil) 4 3 Ka 
( ., \| Ein2il+2iN2 720 \pWeosätdt; nis 
0 m 
 f >) 31 : 
2(f 2 Sinıl \ 
,=— \\emsirz2in:} \v@ )cositdtidä, 


20 2Sinil »z Zug pr 
4=z | \(Cmarn Fame 6 \p@eosätdiidä, 


0 u 


9} P 9. 2 26 r 
2[(j 2Sinil 4 s er 
— en - —— f) cos bs 
-\ \(Sin241+ 2307: 120 \r )cosAtdtr di 
0 0 
Für p= "halt ‚ir durch Auswer : init Hilfe der Si onsel Regel 
ür p’= Coi: t erhalten wır durch Auswertung ınıt Hilfe der Sımpsonschen hege 
foleende Werte: 
a; = — 0,3779, a=002753, © - 0,00257 , 
a= 031136, a =0,03377, «= 0,00254. 


Auch hier kann man sich bei der Berechnung auf das Intervall von /=0 bis »/=5 
beschränken, da der Wert des Integranden für größere A sehr rasch gegen Null geht. 
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Für die Randbelastung können wir die Konstanten auf eine einfachere Art durch 


1 
Co’ t 
METER IN .: 
Differenzieren der Reihe für Gojz; gewinnen. 
1 2 1 3 


Mit Soft 3 Coft " 6 dr Koft 
l 


folgen die Konstanten für die Randbelastung 


p**(h) ori leicht aus denen für p*(#). 





-_2 „,21 . »_# .,68 5 ._2 »,65 . 
a, 3 a, + 6 :@,, Ad, De da, 6 üs a, ii a-+ e Os; 
„_2 -,32 ARTOHEN.. DEF 2... 
a 3.4 a htm 
a; "= — 0,30672, "= —0,03119, as* 0.003883 . 
a = 02434, a"= 0,0808, a}"= 0,00887. 
Die Konstanten für die Randbelastung p %}* a werden wieder durch Auswertung 
mit der Simpsonschen Regel erhalten. 
a, ** = —0,06414, 0,0062, a; "= VOOLTH, 
a**— 00895, a =0,M110, a5"* = — 0,00084. 
Die Bestimmung der Konstanten a}, a;*, a,” erübrigt sich, da diese bei der Spannungs- 


berechnung ohnehin herausfallen. 
Damit lauten die Reihen für die Spannungsfunktion, wenn wir die konstanten Glieder 
gleich weglassen: 


1 
1. Randbelastung ;;;: 
EN GCoj?t 
F*- 0,37749 - r?- c0s2gq 0,02753. r’- cos4tg 0,00257 - 7°. CoS6y ... | 
(32) 
+r:cosgp [031136 - r- cosp + 0,03377 - 7? cos3p +0, W254 r?-cosdp +] 
1 
5) 5 eo em R 
2. Randbelastung Soft’ 
F** = — 0,30672 : r?- cos2g 0,03119 : r!- cos4q VOOBSS- Fr" -cCoSs6p —: - - 
(33) 
+r:-cosp[0,24134- r- cos + 0,03098 - r?- ecos3 7 + 0,00387 - r?- cos5p +: -] 
tIgt 
3. Randbelastung <-;.;: 
Goj?t 
F*** — — 0,06414 - r? - cos 29 + 0,00692 - r! - cos to +0, 00174: r°-cosbp +: - 
(34). 


+ r- cos [0,08935 - r - cosp + 0,0110 - 7? - cos3g 0,00034- r°-cos5gp-+- Jj 


Für die Erfüllung der Randbedingungen am Lochrande ist es zweckmäßig, die Reihen 
in folgende Form zu bringen: 


F=a,+ (a,+4,): r?60s2p + (a + a) r!- 60s4p+(a, + a,)- 7%: 6086P+---) 
: , . 2er (35). 
+r-sınpla, r-sing+a,:r*-sin3p-+ta,:r’-sindp-+--- | 
N ? 3 / 5 
1 
. Ri lastung 5; : 
1. Randbelastung Gojet 
F* = —0,06613- r?- cos2« 0,00624 : r?- cos 4q 0,0003: r® -cos6p — 
ü i i (36) 
DD), 
+ rsing [0,31136 - r - sino + 0,03377 - r? - sin 30 + 0,00254 : r’- sind« SR 
/ / 7 / 
1 
2. Randbelastung — ——;: 
> Goj*t 
F**—= — 0,06538 - r? - cos 29 — 000021: r" cos4gq 0,00001 - 7° - cos 6g ar | 
i s : FUN (37), 
+r- sing [0,24134 - r- sing + 0,03098 - r - sin3g@ + 0,0387: r°- sindp +» J 
tIgt 
3. Randbelastung <-;,;: 
dbelastung Sofrt 
F*** — 0,02521 - r? - cos 29 + 0,00802 » 7? - cost +0,00140 - r" -cos6p +: -- 
A i s a (38). 
+r- sing [0,08935 - r - sing + 0,0110 - 7’ - sin3@ — 0,0034 - r?- sindg + -- -] | 


Damit ist der Spannungszustand des Streifens für die gegebenen Randbelastungen gefunden. 
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Die Spannungsfunktion für den gezogenen Streifen mit einem Loch 
Wählen wir für die Spannungsfunktion für die Scheibe mit periodisch auf der -Achse an- 
so erhalten wir einen Spannungszustand, bei dem dis 
Scheibe in der y-Richtung gleichmäßig gezogen ist. Zur Erfüllung der Randbedingungen am 
Loch nehmen wir wieder periodische Funktionen mit Singularitäten im Koordinatenanfangs- 
punkte hinzu. Wir erhalten geeignete Funktionen, wenn wir zu den bei der unendlichen 
Scheibe verwendeten Gliedern ®, und y®P, die entsprechenden Spannungsfunktionen nach 
Gl. (36) bis (35) addieren, wobei wir diese mit Konstanten multiplizieren müssen, damit der 
Streifenrand spannungsfrei wird. Die auf diese Weise entstehenden Funktionen haben außer- 
dem im Koordinatenanfangspunkte die gewünschten Singularitäten. 


F=N#Insinz gibt im Selmitt «= +1 die Randbelastung p* 5... Wir überlagern 
bol’y 
also in diesem Falle den Spannungszustand F* nach Gl. (36) und erhalten damit die Spannungs 
funktion F,, wenn wir für R In sin z die Reihenentwicklung nach Gl. (7) setzen 
F,=1Inr 0.1054: r?- cos?2aq VOLISO- 17. costg VMOB2:- Tr -Ccosbg 
. om R - (39) 
r- sing [031136 - 7 - sing + 0,03577 - 1° - sindg 000254 Fr’ sın dr 
F=-#R -—— und F yNictez liefern die Randbelastungen 
in? 2 ’ 
| 6 . 2 2ylIgy - 
B°° ETES en bzw » Ben — ._ - 
kol’y bol'y Goj’y 60)? 


Wir überlagern die Spannungszustände nach Gl. (36) bis (35), so daß der Streifenrand 


spannungsfrei wird, entwickeln F in Reihe nach Gl. (8) und (10) und erhalten damit folgende 
Funktionen: 
| | 
. | 


Fr, -Cos2g 019443 - r? cos?ı« V.O36S0 - r?- coste 0.000142» cos6en.. 
3 r? / / f j j (40) 
FHr-sing [0,20260 - » - sing + 0,05080 » 7° - sindg + 0,01306 - 7° - sindg He] | 
F', 018268 - r? - cos24 - I.IO356 - r? - cos 27 -OMW2S6 : r® - cosbg .. | 
1 1 « 4 . 919.3... op  - > . , (41) 
rr-sıng -sıny 011069: r-siny 004312 r°sın 39 0,0056 SIndg . | 


Des weiteren nehmen wir noch die Glieder F, und F, hinzu: 


.) *) .) 


ö 5) > 
F i R —— Sue „ cos hu . , R (+2), 
sın" 2 sın 2 r 
cos2 4. wer 
;snysm5p—+ (43). 


F yNRi 
k J sın"z 


Da diese Glieder in der endgültigen Spannungsfunktion mit Faktoren auftreten, die sehr 
hohe Potenzen von a/! enthalten, können wir hier auf die Überlagerung der Spannungs- 
zustände nach Gl. (36) bis (38) verzichten. 

Mit den auf diese Weise gewonnenen Funktionen wollen wir die Randbedineunzeen am 
Lochrande erfüllen. Damit dieser spannungsfrei bleibt, muß für r=a wieder F=const und 


oF i b » ww 
= gleich Null sein. Wir kombinieren derart, daß die Randbedineungen am Lochrande 
( 
möglichst erfüllt werden. 
Wir machen hierzu den Ansatz: 
F=cx’+a,F,+a,F,+a,F,+a,F,+a,F, . Ehe rar A 


4 


Die Koeffizienten bestimmen wir aus den Randbedingungen für r= «a, also F= const 


oF i , ANTER ' as 
und y- —(), und zwar verlangen wir wieder, daß möglichst hohe Potenzen von (a/T) fortfallen 
07 


Wir erhalten 


a ca? 008996 - ra? + 0.2434 1: ca! 
a,—=1/2- ca +0,42091 - ca® + 0,41555 - ca®, 
a,=2ca? + 1,6364 - ca’ + 1,63572 - ca®, 

a, VO4516 ca’, 


a, 0.36126 - ec a®. 
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Die Spannungsfunktion lautet somit: 
F= ex? (ce a?+0,08996 ce a* — 0,24341 c a®) (Inr — 0,10054 7? cos2 9 — V,011S0 r! cos 4q 


0,31136 r? sin29 — 0,03377 r’sing sind...) 


1 Bis: 
is ca* + 0,42091 ca® + 0,41555 c a® 
„cos2p - 0,1943 r? cos2p — -- +0,20260 vr? sin’ ++) 
y2 
+ (2 ca? + 1,68364 ce a + 1,63572 ce a®) (0,18268 r? cos2 9 + 0,00356 r? cos4p + -- 
- sin’q 0,11069 »* sin? g 0.04312 r? sing sin3® —:::) 
Dun 3 Bro: 
-004516 ca | ‚ecostg a cos24q FE 
F ” 
sn RS RE ER 
0,36126 e a*| „a Sin p sın 37 + gsinipr-- | 


Dieser Ansatz erfüllt die Randbedingungen am Lochrande bis zu den mit «a° behafteten 
Gliedern. Die maximale Spannung erhalten wir daraus wie folgt: 
oF 


Ina = 9,9=0 =, (9 RER ee lee AN: 


Die durchschnittliche Spannung beträgt wie früher: 
On=2ec- BEER, SEE NEE Ba Ze a 9 m 5 2 


Aus diesen beiden Gleichungen ergibt sich die Spannungserhöhung 


9, PEN ri ) .f 
ax — (3 + 1,77361 a? + 2,8066 - a +: (l-afl)). 
Od 


ö a “ ee de e Bo. ra 
Nun setzen wir für die Streifenbreite 27, an Stelle von a den Ausdruck >] |, und 


erhalten 


( 
max _ 3 _30/1-+4,37622 (a|D)? — 4,3762 (a1 | | 
Od y F ; e & ® £ : (48). 


+ 17,10552 (a/D)* — 17,10552 (a/[D’ + - - 


Damit ist die Spannungserhöhung beim Zugstab bekannt. Die Ergebnisse sind in Bild 15 
dargestellt. Die Entwicklung hat jedoch nur Gültigkeit, falls a/l<1 ist. Ist dagegen b/I<1, 
so muß ein anderer Lösungsweg beschritten werden. 


2. Die Stegbreite ist klein im Verhältnis zur Streifenbreite. 


Ist der Durchmesser des kreisrunden Loches im gezogenen Streifen fast so groß wie 
die Streifenbreite selbst, bleiben seitlich nur sehr schmale Stege übrig. Bezeichnet man die 
Streifenbreite mit 21 und die Zugbelastung mit e, so wird in den beiden Stegen eine Längs- 
kraft P=c-1 übertragen. Diese wirkt in jedem Stege etwas exzentrisch. Die Stege bilden 
den elastischen Teil, während die Streifenteile in einer gewissen Entfernung von dem Stege 
nur geringe Spannungen und Formänderungen haben. Betrachtet man den starren Teil mit 
dem daran anschließenden Stege bis zum Querschnitt an der schmalsten Stelle desselben, so 
muß die Längskraft P so wirken, daß dieser Querschnitt eben und parallel zur Anfangslage 
bleibt. 

Wir geben nunmehr ein Spannungsproblem an, das dieselben örtlichen Verhältnisse auf- 
weist und das theoretisch behandelt werden kann: 

Eine Halbebene wird in der y-Richtung gezogen. Nahe am Rande befindet sich ein 
kreisrundes Loch. Aus dem so entstandenen Stege wird ein schmales Stück durch Schnitte 
parallel zur - Achse herausgetrennt und der Steg danach gewaltsam wieder verbunden. Durch 


k den Zug der Halbebene und das Zusammenfügen des Steges entsteht in letzterem mit großer 
H Annäherung der gleiche Spannungszustand wie in den Stegen des gezogenen Streifens, falls 
{ diese schmal genug sind (Bild 11). 

| Zur Ermittlung der Spannungsfunktion wird zweckmäßig der Weg der Inversion oder 
i Umkehrung beschritten ’). 


j #4, Z2. B. C. Weber: Z. angew. Math. Mech. Bd. 2 (1922), S. 185 bis 187. 
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zu j : n : e F(r,9) Br 
Setzt man in der Airyschen Spannungsfunktion = F 
a’ y f x 1 
C=-73 35 Y= ; 5: )ZW. Y=@_; =, 
2? +y® „ x’? + y° / / r 


so erhält man eine neue Airysche Spannungsfunktion F’. Jedem Punkte der Funktion F 
entspricht ein Punkt der neuen Funktion F’, der auf einem gleichgerichteten Strahle liegt, 
dessen Entfernung vom Umkehrpole jedoch den umgekehrten Wert hat. Beliebige Kreise 
gehen wieder in Kreise über. Durch nochmalige Umkehrung erhält man wieder die ursprüng- 
liche Funktion F'. 

Bezeichnen o,, o,; und r die Spannungen des durch die Funktion F gegebenen Spannungs- 
zustandes mit dem Umkehrpol als Koordinatenanfangspunkt, ebenso o}/, 0; und r’ die durch 
F’ gegebenen Spannungen, so besteht für entsprechende Punkte beider Funktionen folgende 
Beziehung: 


er m“ 
o,=r"?o,+2|F’— r’ Ir ‚ g=r"oa+2|F’—r' 7 z - 
Hieraus folgt: Die Hauptspannungen gehen wieder 
in Hauptspannungen über. Den Hauptspannungs- 
linien (Spannungstrajektorien) der Funktion F ent- 
sprechen ebenfalls Hauptspannungslinien der um- 
gekehrten Funktion F’. Ein gleichmäßig belasteter 
bzw. unbelasteter Rand gibt bei der Umkehrung 
wieder einen gleichmäßig belasteten Rand. 

Bei passender Wahl des Umkehrpoles € gibt 
die Umkehrung der Halbscheibe mit einem Loch 
am Rande eine durch zwei konzentrische Kreise 
gebildete Ringfläche. Bild 11 zeigt die Halbebene, 
Bild 12 die Ringfläche. 

Der Punkt B habe in der x-y- Ebene vom Rande der Halbscheibe den Abstand 1/b, 
wobei b die Stegbreite AM bedeutet. Wählt man jetzt den Umkehrpol € in der Entfernung I 
vom Rande, so erhält man folgende Beziehungen : 








Bild 11. Bild 12. 


1 E; Fr XB l b ; 

2. BTL u 1b’ 
1 

1: 2ı=1-b; FE 75 

H LO; 17 0 

M: Xu > 1 4 ey - 1 R 


Der Rand der Halbebene ergibt also den inneren Kreis, während der Lochrand in den 
äußeren Kreis übergeht. 

Durch eine Verschiebung des Koordinatensystems in der «’-Richtung um u’ = 1/2 läßt 
sich weiter erreichen, daß der Mittelpunkt der Kreise mit dem Koordinatenanfangspunkt zu- 
sammenfällt. Die Transformationsformeln sind: 


Der Spannungszustand der Scheibe setzt sich bei einer Zugbelastung im Stege b aus 
zwei Anteilen zusammen: 
1. Anteil, Fall a: im Stegquerschnitt greift nur ein reines Moment an. 
2. Anteil, Fall b: im Stegquerschnitt greift nur eine Längskraft mit dem Angriffs- 
punkte in M an. 


Für die Spannungsfunktion F und deren Ableitung nach der Normalen ergeben sich 
daraus folgende Randbedingungen, wenn man beachtet, daß die beiden partiellen Ableitungen 
der Spannungsfunktion nach x und y die Summe der y- und «-Komponenten aller Kräfte 
zwischen einem beliebigen Anfangspunkte A und einem Bezugspunkte S des Randes bedeuten, 
während die Spannungsfunktion selbst ihr Moment darstellt. 


m 


. i ( =. 
Im Falle a und b ist auf dem Lochrand F=0 und es gesetzt, so daß dieser 
on 


- a u O4 te 
spannungsfrei ist. Auf dem Rande #=1 müssen im Falle a F=(C, und >. —0, im Falle b 
5 n 
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dagegen F=0 und sr Ü, sein. (', ist hierbei die im Steg übertragene Längskraft. 
- 
Aus dieser ergibt sich die durchsehnittliche Spannung im Steg 0, =— : 
j ) 


Diese Randbedingungen ändern sich für die durch die Inversion erhaltene Spannungs- 
funktion F’(#’, y’) bzw. F', (#,,9,). Der Rand »—=1 wird zum inneren Kreis, r=r;, während 
der Lochrand in den äußeren Kreis, r=r,„, übergeht. Für F, ergeben sich damit folgende 
Bedingungen: 


i ER OR, 
Fall a: 1.r=r;: F,=—0,(csp, +1), ey C,(cospy, +). 
1 
YF 
( 
> mw „ 1 
2. r=r.:F,=0 Ir, 0 
\F 
y , ( j7 
Fallb:1i.r=r:F,=0, 7—: C, 
or, . 
3F 
( 
BI u 1 
2: — 08, u 0. 


1 


Von der Richtigkeit dieser Bedingungsgleichungen kann man sich dadurch überzeugen, 
daß man Spannungsfunktionen wählt, die den oben für das & — y System angeschriebenen 
Randbedingungen genügen, und bei diesen die Umkehrung vornimmt und zwar: 

Für Fall a: ie ER 

Diese Funktion genügt den gestellten Bedingungen, da ja für z-=1 F 
oF oF a = en . 

3,0 ist. Nimmt man die Umkehrung vor, so erhält man: 
( 


—=0, und 


or; 
TERROR BEN i Ei. 
F 0, „=Uur,; und F=\n+2,+7)6, 
OF, ) ( 
= (Zr, (COS 9o,) 5 
or, EN 


Der Rand #=1 geht in den Kreis r; =>; über. Damit lauten die Gleichungen 


of ' 
r, ( tie, "o9- und h) 7 C, (cos Y%7r | 
| 
f 1 5) (cosYp,+ 1) 
Für Fall b: F=(,(—]), 
. . Ö o . . 
ür 2=1 ist also F=0 und \ ‚. Die Inversion ergibt 
ol) & 
F'F=F:.r"=(,-(x —r’”) 
daraus 
B 1 1 
F=C(2+35 N; I, f 
E 1 
F,= c,(r rc)» 
OF, 
0,233 
O r, 2 1 
i ; a le ; 1 
Auf dem inneren Kreis ist also mit nr, =5 
F, 2 VÖ, 
OF, ’ 
= E 
or, . 


Aus den Randbedingungen lassen sich folgende Schlüsse für den Aufbau der Spannungs- 
funktion ziehen: Im Falle a besteht F, aus rotationssymmetrischen Gliedern und solchen, die 
cos, enthalten. Im Falle b dagegen setzt sich F, nur aus rotationssymmetrischen Gliedern 


zusammen. 
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Gewählt werden daher folgende vollständige Ansätze: 


a FF =a+a,r+ta,riilnr;+a,inr +a,r, cosyg, +a,ricosg, 
Li 
cosq, d (49). 
+ a- +a,r,cosy, nr). 
7 


1 
b)F=b, +b,r.+b, br +b,lhr : . -. .. 2220. (50) 


Die einzelnen Anteile sind Potentialfunktionen oder mit r* multiplizierte Potential- 
funktionen. Die Differentialgleichung JA IF, —=0 wird also erfüllt 


Um den endgültigen Spannungszustand zu erhalten, muß man Fall a) und Fall b) 
überlagern. Zur Bestimmung der Konstanten erhält man also 12 Gleichungen aus den Rand- 
bedingungen. Eine weitere ergibt sich aus der Forderung, daß dem Verschiebungszustand 
nach der Überlagerung ein Spannungszustand entsprechen soll, der mit dem der gezogenen 
Scheibe übereinstimmt. Die mehrdeutigen Anteile der Verschiebungen in der y- Richtung 


1 | n 
I=7 \(es+ = o.)dy haben folgende Form: 


Y Y : 
V®, arcig —; v,=xc.-arcig = " x — 2) -arctg 
u - E X x 


i A r 4 Be 
Diesen entsprechen folgende Verschiebungszustände: Arc tg ist der Winkel zwischen 
.Z 


der &-Achse und der Verbindungsgeraden des beliebigen Punktes (x, y) mit dem Koordinaten- 
anfangspunkte. v, ist also für die Punkte der »-Achse Null. Geht man jedoch im Sinne wachsender 
Winkel 9 um das Loch herum, d. h. vergrößert man @ von O bis 27, so kommt man wieder 
auf die Punkte der positiven z-Achse. ®, hat nun jedoch den Wert 27. Die Scheibe schiebt 
sich also um den Betrag 27 übereinander. Dieser Spannungszustand kann so entstehen, daß 
ein schmaler Streifen von der Breite 27 aus dem Steg herausgeschnitten und die Scheibe 
gewaltsam wieder verbunden ist. Es handelt sich also um eine Vorspannung der unbelasteten 
Scheibe. Die Spannungsverteilung im Steg stimmt mit 


der der gezogenen Scheibe überein (Bild 13). hy, 
v, dagegen liefert mehrdeutige Verschiebungen, die 
proportional der Abszisse x sind. Es entsteht dadurch 
eine keilförmige Aufspaltung vom Koordinatenanfangs- eh n 





punkt aus. Diesem Verschiebungszustand entspricht ein 
Spannungszustand, der nicht mit dem des vorliegenden Bild 13. 
Falles übereinstimmt (Bild 14). 

v, liefert mehrdeutige Anteile, die proportional 
(e 2) sind. Die Aufspaltung berührt also das Gebiet der Halbscheibe nicht. 








Bild 14 (rechts). 


Hieraus ergibt sich eine weitere Bedingungsgleichung für die Konstanien des Ansatzes, 
da die Summe aller ®,- Anteile Null sein muß. 


Es ergibt: 


ö F Y 
F,=Inri, „= —8Sx.arctg 7, 
A . \ Yy 
Pr =ernlbor, p,-: Sr. ac, 

2 a 
; R Y 
Fr, 89, lar;, 0,42. anig”. 
Fon 


Damit keine keilförmige Aufspaltung auftritt, müssen also die Konstanten folgender Be- 
dingung genügen: 


er rrtet+ea) - ee .. ... rn es WB. 


Die Ermittlung der Spannungserhöhung wird im folgenden durchgeführt. Dabei ist es 
nicht nötig, die Konstanten der Spannungsfunktion einzeln zu bestimmen. Das Verhältnis 
I max BIER ° . . . - . . ne. . 

u läßt sich vielmehr aus den Bedingungsgleichungen für die Koeffizienten unmittelbar 
Id 
gewinnen. 
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Die Bedingungsgleichungen lauten: 


a) aus den Randbedingungen: 


l 
u. +, N rs, ar +rahbri=>C .....:. 0.2.2.0. 188), 
e., ra. ts hr, tar, =®.’. ... re ° 
) 
2r;a,+ta,(2r;inr; +2r,)+ = AT rer 
rg 3 ı i ı Y; I) 1 
q 
») 
2. +a,2r.iar, +2r)+ 7 SEITE 
a 
i cos Yo, . B -6 
d,1;c0oSp, t+a,r} cosy, + a, z +, so hri=—C,. . . . . . (CM, 
i 2 
i cos y, £ nei 
d,F4c0oSsp, ta, r, cosp, +4; ie + a, r.Cosp, In e a a A 
a 
2. cos Yp, 9 ) 1 -o 
a,cosp, ta,3r}cosp, —d,- „ t+ascosg, Wi FIDEL, . .., ., 68, 
q 
PR‘ cos Yı 2 »\ ri 
a,cosp, ta, 3rc0osyp, —A, a tWcosp, ma +Dd=0 .. ..: 
a 
rt anne ..:. . re 
EIER EIEBESE. : .:, 2 2 er ec 
) 
25,1; +b, Zr; Inr}?+2r)+ 2 nn ES UER ee Pe ee 4° 
ı 
> 
2b, ru +b, 2 r,„Inr, +2r.)+ E ee N ee Tr 
a 
b) Aus den Verschiebungen nach Gl. (51): 
+) ri, Hi) -2u,=0 .... 2002. (A) 
r N eo 0, P PT, . . . . 
Zur Berechnung der Spannungserhöhung -""* benötigen wir die Spannung o, im 
e Id ö 


„F 
Ei ’ ’ i Ri ' 
Punkte B(«=r„,9=0), im folgenden mit o, bezeichnet. Da im Punkte Bo,, = dt gleich 


Ji 
Null ist, wird 
2, oe 
— 


0, da? Toy} =AF, (65). 
Die Potentialfunktionen liefern demnach keinen Anteil zu o,. Damit wird 
o,=4Ala,#x, ri +a,&, mr) + (a,+b,) r!+(a,+b,)-riInr; 
(66). 


t r 
0,=8r.:a,+ or +4(a, + b,)+(a,+b,) (4Inr, +8) 
a 
Wie man aus o, die Spannung am Lochrande im #-y-System berechnet, wird später 
gezeigt. 
Aus den Gleichungen werden nunmehr sämtliche Unbekannte außer o, und C, beseitigt. 
Durch schrittweises Eliminieren erhalten wir schließlich: 
0, 4 3 1 2 
= i — ha a a ale 2 er 
2 l „? ] Se . Ya Ta 
nn nz |ra 


Mit €, ist die durchschnittliche Spannung im Steg gegeben: 


a RE u Er 


Die maximale Spannung am Lochrande wird nun wie folgt aus o, berechnet: 


3? F 
On“ a ea ar Sr 
y—ı 


Ö x” 
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or Be a5 Set; 
ferner ist F„_u—r? F’| ,‚ae=->— : und mit F,-o=r? F’(—-) ist 
x r x r r 
0° 47 0 | e P* l } 
j ar 
“max 99 or? | r/} 
(70). 
WESER 6... Be 
"max N r or r 01 
N 


Da nun auf dem äußeren Kreisring sowohl F’ als auch \y' gleich Null sind, wird 
( 


Ir 1 " 
oO jr: »- au EEE } i 6 A 
Am Lochrande ist r=1—b 
| 
Ö _ Ö ») 
i (1 b) ; 1,27. 
Mit Gl. (67), (68) und (72) ergibt sich die Spannungserhöhung 
I max 6, r b (73) 
o Ad-» TC, ‘3 
und mit Gl. (67) 
Een 2 l [2 T 
m ar = “ - (14). 
Od Ya In2r,„ (1—b) 


Diese Gleichung ist noch nicht dimensionsrichtig, da im Bild 11 die Länge I eingefülhrt 
ist. Zur Richtigstellung der Dimensionen wird folgende Beziehung benutzt: 


ET DZ . (5), 
ferner 
| b I 1+b hr 
Ya 5) we h 2 2 = 1 u . . . . - r . ö r ’ . (6). 
Damit geht Gl. (74) über in 
N. FR tb 1 BER 
= 2° (“i). 
Gl (1 b)° 1+b 
In me 
Nach Gl. (75) ist (1—b)”’=2ab, also 
RR 2 1 2y(2a+b)b 1 ” 
1 ni (19). 
Od a I I+b a ' l+ b: Z2a+b)b 
ib "| b:y(2a+b)b 
ei Zu ER : ; SE _ y “ 2 
Gl. (78) läßt sich in Reihe entwickeln, da ja z l ist. Wir erhalten mit a=1-—b 
On _a 2 b ,3/b\ , 536 b \ 70 
stehe... 


Mit Gl. (79) ist der Wert der Spannungserhöhung für große Lochdurchmesser bekannt. 


3. Ergebnis. 


Zusammen mit dem Ergebnis für kleine Lochdurchmesser können wir nun den Ver- 
lauf der Spannungserhöhung . 
I max 
Od 
in Abhängigkeit von a/l bzw. b/l angeben. 
In Bild 15 ist das Ergebnis zusammengestellt. 
Für kleine Werte a/l ist die nach Gl. (48) 
gefundene Spannungserhöhung eingezeich- 
net (gestrichelt). Für kleine Werte b/l ist 
der Endwert und wahrscheinlich auch die 
zugehörige Tangente richtig (gestrichelt). 2/7-2% 
Dazwischen läßt sich der vermutliche Ver- 
lauf einzeichnen. 413 Bild 15. 
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Über die Anwendungen der allgemeinen Bewegungs- 
gleichungen starrer Körper in bewegten Bezugssystemen”. 
Von K. Magnus in Darmstadt. 


Bei der Berechnung von Kreisel- und Pendelygeräten muß man häufig die 
allgemeinen Bewegungsgleichungen starrer Körper verwenden, bei denen 
die Bewegungen der zu betrachtenden Körper ron einem selbst irgendwie 
bewegten Bezuyssystem aus beurteilt werden. Unter Anwendung des der 
lheorie starrer Körper besonders angepaßten Hilfsmittels der Motorrechnung 
werden die verschiedenen Formen dieser Gleichungen und die Zusammen 
hänge mit ihren entsprechenden Vektor- und Skalarfassungen behandelt 
Nach einer Erweiterung der Gleichungen auf Systeme starrer Körper werden 
die Ergebnisse an den Beispielen eines physischen Pendels und eines kräfte- 
freien Kreisels mit massebehafteter kardanischer Aufhängung veranschau 
licht. Da jedoch die strengen Gleichungen nur in ganz wenigen. besonders 
einfachen Sonderfällen zu lösen sind, ist es wichtig. daneben Näherungs- 
methoden zur Verfügung zu haben, mit deren Hilfe auch in schwierigeren 
Fällen noch Aussagen getroffen werden können. Die Grundlagen dieser 
Näherungen sowie ihre Ableitung aus den strengen Gleichungen werden 
gezeigt. 


1. Einführung. 


Bei der Berechnung von Pendel- und Kreiselgeräten ist man häufig vor die Notwendigkeit 
gestellt, die allgemeinsten Bewegungsgleichungen für den starren Körper relativ zu einem ge- 
führten Bezugssystem anzuwenden, um das Verhalten der Geräte auf bewegtem Fahr- oder 
Flugzeug zu bestimmen. Die allgemeine analytische Form dieser Gleichungen ist seit langem 
bekannt und wird z. B. von Routh [1]') sowie Herglotz [7] in impliziter Form gegeben, bei 
der die Impulse sowie ihre zeitlichen Ableitungen im Führungssystem auftreten. Sämtliche 
Koordinaten der vorkommenden Vektorsysteme (der Kraft, der absoluten und der relativen 
Geschwindigkeit) werden dabei im Führungssystem angeschrieben. Man kommt bei der Ab- 
leitung dieser Gleichungen mit zwei Systemen, dem Inertialsystem und einem Führungssystem 
aus, das sich gegenüber dem Inertialsystem in bekannter Weise bewegt. Die noch einfachen 
dynamischen Gleichungen verlieren jedoch an Übersichtlichkeit, wenn man die Ausdrücke für 
die Impulskomponenten als Funktionen der Geschwindigkeiten des Körpers einsetzt und so 
die expliziten Differentialgleichungen für die Komponenten der Geschwindigkeit ausrechnet. 
Diese Arbeit ist von Heun [2] geleistet worden, der die Vektorform der Gleichungen für den 
Fall angibt, daß die zeitlichen Änderungen der Körpergeschwindigkeit ebenfalls vom Führungs- 
system aus beurteilt werden. Der in den Gleichungen auftretende Trägheitstensor besitzt 
daher zeitlich veränderliche Glieder, so daß die dynamischen Gleichungen durch eine Diffe- 
rentialgleichung für den Trägheitstensor vervollständigt werden müssen. 

Eine sehr durchsichtige und begrifflich klare Darstellung der allgemeinen Bewegungs- 
gleichung hat R. v. Mises [3] gegeben, indem er das für die Theorie der starren Körper be- 
sonders geeignete Hilfsmittel der Motorrechnung anwansdlte. Die v. Misessche Darstellung 
wurde von Winkelmann und Grammel [4] übernommen und in bezug auf die explizite 
Ausrechnung weitergeführt. Die Beurteilung der Impulsänderungen wird dabei von einem 
körperfesten System vorgenommen, in dem der Trägheitsmotortensor zeitlich konstant ist. 
Die in Motorform hingeschriebenen Bewegungsgleichungen sind zwar unabhängig von der 
Wahl des Koordinatensystems, sie erweisen aber ihren vollen Wert nur, wenn die Kompo- 
nenten bei der Zerlegung sämtlich im körperfesten System genommen werden, da nur dann 
der Trägheitstensor zeitlich konstant ist. Es müssen also auch die Komponenten der Führungs- 
geschwindigkeit im körperfesten System genommen werden, was von vornherein nicht möglich 
ist, da die Bewegungen des Körpers ja erst bestimmt werden sollen. Man hat also zu den 
dynamischen Gleichungen noch die Transformationsgleichungen für die Komponenten der 
Führungsgeschwindigkeit hinzuzufügen. 

Diese Schwierigkeiten umgeht Hölder [6] dadurch, daß er die dynamischen Gleichungen 
wieder im Führungssystem anschreibt. Er erzielt durch Verwendung Plückerscher Koordi- 
naten starke formale Vereinfachungen, jedoch sind seine Gleichungen für die Anwendung nicht 
so durchsichtig wie die Motorschreibweise. 

*, Diese Arbeit wurde von der mathematisch-naturwissenschaftlichen Fakultät der Universitat Göttingen als 
Habilitationsscehrift für das Fach der angewandten Mechanik angenommen. 


1) Die in eckigen Klammern stehenden Zahlen verweisen auf das Sehrifttumsverzeiehnis am Schluß der Arbeit, 
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Im Anschluß an die von Winkelmann und Grammel gegebene Darstellung [4] sollen 
im folgenden insbesondere die Möglichkeiten der Zerlegung der allgemeinen Motorgleichung 
in Vektorkomponenten dargelegt werden, wobei sich Gelegenheit bieten wird, auf die ver: 
schiedenen Formen der allgemeinen Motorgleichung selbst hinzuweisen, die es ermöglichen, 
die Ausgangsgleichung dem jeweils vorliegenden Problem anzupassen. 

Wenn auch v. Mises die Zerlegung der allgemeinen Motorgleichung in ihre Vektor- 
bzw. Skalarkomponenten als eine rein schematische Arbeit bezeichnet, so haben doch bereits 
Winkelmann und Grammel sowie auch Hölder darauf hingewiesen, daß die Zerlegung 
neben schematischer Rechenarbeit besondere Vorsicht und Gedankenarbeit erfordert. Durch 
eine Betrachtung der Zusammenhänge der in den verschiedenen Systemen genommenen zeit- 
lichen Ableitungen der Geschwindigkeitsmotoren sowie ihrer Vektorkomponenten und Dar- 
stellung in Tabellenform soll die Zerlegung der Motorgleichungen nunmehr soweit vorbereitet 
werden, daß die Zerlegung tatsächlich nur noch schematische Rechenarbeit erfordert. Für 
eine von der Darstellung bei v. Mises und Winkelmann-Grammel nur wenig abweichende 
Form der Motorgleichung soll die Zerlegung in Vektorkomponenten ganz allgemein durch- 
geführt werden. In einigen praktischen Beispielen wird die Ausrechnung bis zur Angabe der 
Bewegungsgleichungen in Skalarkomponenten erfolgen, um so die Anwendbarkeit der all- 
gemeinen Gleichungen in Spezialfällen zu zeigen. 

In Erweiterung der bisher vorliegenden Ergebnisse soll die allgemeine Bewegungsgleichung 
auch bei Beurteilung der Impulsänderungen vom Führungssystem aus als Motorgleichung an- 
geschrieben werden. Dann fällt zwar der Vorteil eines zeitlich konstanten Trägheitstensors 
fort, jedoch gewinnt man auf diese Weise eine Form, die als Ausgangsgleichung für die 
technisch wichtigen Näherungsrechnungen von Bedeutung ist. Die Notwendigkeit derartiger 
Näherungen, die auf eine Beschränkung der Betrachtungen auf „kleine Schwingungen“ hinaus- 
laufen, ergibt sich aus der Forderung, auch bei komplizierteren Geräten noch Aussagen zu 
treffen, denn die strengen Gleichungen bieten im allgemeinen unüberwindliche Integrations- 
schwierigkeiten. Desgleichen erfordern die Belange des Praktikers eine Erweiterung der all- 
gemeinen Gleichungen auf Systeme von starren Körpern, wie sie bei kardanisch gehängten 
Pendeln und Kreiseln sowie bei Mehrkreiselgeräten häufig vorliegen. 


2. Die verschiedenen Formen der Bewegungsgleichungen in Motorschreibweise. 


Nach dem Vorbilde von Winkelmann und Grammel bezeichnen wir im folgenden 
zeitliche Ableitungen im Inertialsystem durch einen darübergesetzten Punkt, im Führungs- 
system durch einen Ring und im körperfesten System durch einen Stern. Die bei den Betrach- 


tungen benötigten Systeme bezeichnen wir entsprechend mit 2” (Inertialsystem), I (Führungs- 
system), 2'* (körperfestes System). Den Motor der Relativgeschwindigkeit (I* gegen 2) be- 


zeichnen wir mit ®, den Motor der Führungsgeschwindigkeit (2° gegen 2°) mit $. Die 
Absolutgeschwindigkeit des Körpers (2* gegen I) sei B=-W’+%. Das System der auf den 
Körper wirkenden äußeren Kräfte sei durch den Motor ® repräsentiert. Mit dem Trägheits- 
motortensor T kann man nun das allgemeine mechanische Grundgesetz schreiben: 

a ae ae 
Um die Tatsache auszunutzen, daß T in 2* konstant ist, wird (1) vermöge der Transformations- 
formel ?) 

A-A+[PA 
in 
TO+HBTB)—=-P 


umgeformt. Daraus folgt wegen T=0: 


TÖ+H[BTBD]I-P. 


Mit B=9%+ 7% bekommt man damit die Bewegungsgleichung in der Form: 


TO +H[B(TBO)+HTZ L FTHHPTH + TFTP)] = B ir Fr 


Dabei bilden die ersten beiden Glieder den Anteil der Relativbewegung, die mittleren beiden 
Glieder den Anteil der- Führungsbewegung, die letzten beiden Glieder der linken Seite geben 
den Coriolisanteıl. Die Grundgleichung (2) läßt sich noch in andere Formen bringen, wenn 
man die zwischen den verschiedenen zeitlichen Ableitungen bestehenden Beziehungen berück- 
sichtigt. Man bekommt für die Ableitungen von ® unter Berücksichtigung von [VW]=0: 


2) Siehe [3], S. 195. 
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* Oo, 
NL RE ee er ee 
Er, 
°© . 
ee . (3e). 

Entsprechend für die Ableitungen von unter Berücksichtigung von [5 5] =0: 

\ 5 r - 353 

. ° 
a EN a en a ee 2a 2 


Ban aa 
Be ie 


Je nachdem, in welchem System die zeitlichen Änderungen der Relativ- und Führungs- 
geschwindigkeit beurteilt werden sollen (bzw. gemessen oder gerechnet sind), kann man die 
Beziehungen (3) bzw. (4) in (2) einsetzen und kann so zu insgesamt neun verschiedenen 
Formen der Grundgleichnng kommen, von denen hier nur noch die beiden wichtigsten an- 
gegeben werden sollen. Will man die Änderung der Führungsgeschwindigkeit von 2° aus 
beurteilen, so kommt man wegen (4c) zu der von v. Mises und Winkelmann-Grammel 
angegebenen Form: 


TEHB(TO)HTFHETDIH FTD HFCTCSHTMETD-P 0, 


bei der gegenüber (2) noch ein Glied im Coriolisanteil hinzutritt. Will man schließlich die 
\ 


Geschwindigkeitsänderungen in 2” haben, so kommt man wegen (3b) und (4a) zu: 


oO 


TE+HBS(TBHTTEHSETCDIHFTDHETÖI-P. 0.0.00... 06). 


Diese Gleichung stimmt im Aufbau mit (2) überein, nur daß an Stelle der Sterne jetzt 
Punkte stehen. 

Die ursprüngliche Transformation der linken Seite von (1) auf das System 2°* brachte 
den Fortfall der Glieder mit den zeitlichen Ableitungen des Trägheitsmotortensors mit sich. Da- 
durch ließen sich die weiteren Formeln in besonders einfacher Weise ableiten. Es ist jedoch 
für manche später zu besprechenden Fälle zweckmäßig, auch eine in 2° angeschriebene 
Motorgleichung zu haben, weil die bei der Zerlegung der Gln. (2), (5) und (6) notwendige 
Transformation der Komponenten von & nach 2* umständlich ist und die Anwendung der 
Gleichungen erschwert. Andererseits würde man durch das Anschreiben der Komponenten- 


gleichungen in 2° oder 2” den Vorteil eines konstanten Trägheitstensors verlieren. Wir ge- 


winnen die gewünschte Motorgleichung durch Transformation von (1) nach 2° vermöge der 
Transformationsgleichung: 

® ° 

A=-A+[FA 
zu 


_o 
TEHFTBI-P 
oder: 
u 
TU’ +HTE+HSITPHFTDI=-P : : : : 22222 M. 
Es sei beiläufig bemerkt, daß man aus dieser Form sofort wieder auf die allgemeine Motor- 
gleichung (5) kommen kann, wenn man die Ableitungen der Impulsanteile wirklich ausführt: 
© o o 
Td=-TPH+T®P, 
T5=-T% +T3: 


Unter Berücksichtigung der für die Ableitungen von Motortensoren geltenden Rechenregel: 





T=[®T] - [TB] 9) 


sowie unter Berücksichtigung der Regeln für dyadische Produkte zwischen Motor und Motor- 
tensor ?): 


TO -(OT|O)- (TB) -|P(TPO))], 
TF- (BOTH) -(TB]H =IS(TH]+ TIER) 


3) Siehe [3], S. 175 ff. 


4, Siehe z. B. [3], S. 196. 
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bekommt man aus (7) die Motorgleichung: 
o je} 
TE’ +[P(TP) TEST HP TAI HF HTIEP)= BP .. 08, 
die wegen (3a) mit (5) identisch ist. 

Wir wollen Gl. (7) mit Rücksicht auf spätere Anwendungen noch durch Einführung 
der Impulsanteile der Relativbewegung $=T® und der Führungsbewegung Y—=T% auf die 
einfache Form bringen: 

je} je} 

BI. . . re ER 
Die Komponentenzerlegung dieser Gleichung führt wieder auf die anfangs bereits erwähnten 
impliziten Formen der allgemeinen Bewegungsgleichungen, wie sie in einigen Lehrbüchern zu 
finden sind. Nur ist dort meist an Stelle der Impulsanteile % und %, mit dem Absolut- 
impuls $+%, und dem Relativimpuls % gerechnet. 

Damit ist das im folgenden gebrauchte Formelmaterial abgeleitet, so daß wir diese Über- 
legungen abschließen können, um uns weiterhin der für die Anwendungen wichtigen Zerlegung 
der Motorgleichungen in Komponenten zuzuwenden. 


3. Die Zerlegung in Vektorkomponenten. 


Der Motor ® der Relativbewegung habe die Resul- sa , /B 
tantkomponente u’ (Drehung von 2* gegen 2°) und bezüg- I \ Pa 
lich des Anfangspunktes 0° von 2* die Momentenkompo- \ is / 
nente v’-v’ ist also die Relativgeschwindigkeit von 0’ in z “ uch 
2°, Der Motor der Führungsgeschwindigkeit % habe die \ \ 
Resultantkomponente u (Drehung von I°° gegen 2”) so- \ \, 
wie die Momentkomponente d ebenfalls bezüglich 0’. v BEE Fr 
ist demnach die Absolutgeschwindigkeit desjenigen in | I u 
2° festen Punktes, der augenblicklich mit 0° zusammen- a er 3 
fällt. Des weiteren benötigen wir für die Zerlegung die er 
Matrizendarstellung des Trägheitsmotortensors. Nach Y /z ws: 
v. Mises°) hat man: Bild 1. 
m 0 0 0 mc mb 
0 m 0 ME 0 ma 
0 0 m mb ma 0 
T= ee 
0 mc mb I; Dr PD, 
me 0 ma PD, I, PD, 
— mb ma 0 Dy D, 7 








Darin ist m die Gesamtmasse des Körpers, a, b, e sind die Komponenten des Vektors tr, 
vom Ursprung 0’ zum Körperschwerpunkt, © sind die Trägheitsmomente, ® die Deviations- 
momente des Körpers bezüglich der Achsen A, B, C von 2*. Das rechte untere Viertel der 
Matrix (10) ist das Schema des Vektorträgheitstensors 


| 9, D- Dp 
T=! — Dr Op D, 
| -D) -P, a7 


den wir im folgenden häufig brauchen werden. 

Der eigentlichen Zerlegung soll nun eine allgemeine Betrachtung vorausgeschickt werden, 
die die Zusammenhänge der Vektorkomponenten der verschiedenen zeitlichen Ableitungen 
betrifft. Aus jeder Motorgleichung lassen sich zwei entsprechende Vektorgleichungen folgern. 
Aus der Gleichheit zweier Motoren folgt die Gleichheit ihrer Resultant- sowie ihrer Moment- 
komponenten, wobei sich die Momentkomponente d, bezüglich eines Punktes 0 aus den „physi- 
kalischen Bausteinen“ des Motors, seiner „Öffnung“ u sowie seiner „Länge“ dv, nach 

a Be a 
zusammensetzt. - Dabei ist r, ein von O nach einem beliebigen Punkte der Motorachse ge- 
zogener Vektor. Bei der Betrachtung der Vektorkomponenten von in verschiedenen Systemen 
abgeleiteten Motoren ist jedoch besondere Vorsicht geboten, da Motoren bei zeitlicher Ver- 
änderlichkeit im allgemeinen nicht in ihren Momenten übereinstimmen. Z. B. folgt aus der 
Gleichung 3, (siehe (4a)) zunächst nur die Gleichheit der entsprechenden Ableitungen für 
Länge und Öffnung: 


5) Siehe [3], S. 170. 





i 
H 
j 
i 
i 
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AU, 
VD, zu D,, 


nicht aber die Gleichheit der Momente. Vielmehr stellt man durch Differenzieren und Sub- 
trahieren von (11) leicht fest, daß 
ee A) Pr 


im allgemeinen von Null verschieden ist. Wir wollen in (12) noch x 
achten zu diesem -Zweck die Transformationsgleichung: 


„ eliminieren und be- 


Wet +[ur], 
die wegen (11) auch 
u d=/[ur]=p, od 


0 


geschrieben werden kann. Setzt man das in (12) ein und berücksichtigt, daß die Vektoren v, 
und u stets gleichgerichtet sind, also [vb, u]=0 ist, so hat man: 


d,=Dd,+[ud,]. 


. OS . r . 
Aus der Motorgleichung $=% folgen also «die beiden Vektorgleichungen : 


ng a ee NA 
ud, +[uv,] J 


In gleicher Weise wollen wir die den Motorgleichungen (3a) bis (4e) entsprechenden Vektor- 
gleichungen ermitteln. Die Resultantkomponenten folgen dabei einfach durch schematisches 
Anwenden der Rechenregeln für Motorprodukte®). Die Momentkomponenten dagegen erfordern 
Aufmerksamkeit. 

Aus (4b) folgt zunächst: 


du =d,+[Wv,]+liul . .» . . . Der yan Ar, Fi 
denn die schematische Zerlegung darf hier nur auf die „Länge“ angewendet werden. Anderer- 
seits kommt durch Differentiation von (11) in 2’ und 2°: 


d—=d,+li,ul+ [vu]. 
++ Ti]. 
Durch Subtraktion dieser beiden Gleichungen folgt unter Berücksichtigung von (14): 


dv, —-d,—= [u’d,)+ [vu] + [i, - t)u]) + [r, ltü -— W]. 
Darin ist nun: M 
u u=[u’ u] (Gleichung der Resultantkomponenten von (4b)), 
ti, - L—=[(u + w)r,) (Beziehung zwischen den zeitlichen Ableitungen eines 
Vektors in Systemen, die sich mit u + u’ gegen- 
einander drehen), 
d„,=d,+[ur)] (folgt aus (11)). 
Durch Einsetzen dieser Beziehung folgt schließlich 
. * 
d,=Dd,+[wd,]+[v} ul+[uD,]. 
wobei von der Beziehung 
Wirte . ....... 0. MB, 
Gebrauch gemacht wurde und wie früher [u r,]=v,— dv, gesetzt wurde. Wegen v,„|u fällt 
dann das Glied [v, u] heraus. 
In ganz entsprechender Weise folgt aus (4e) zunächst: 


d,—=Dd,+[u’v,]+[v/ u]. 


Dies hat man wieder in eine Gleichung einzusetzen, die sich durch Differentiation von (11) 
in 2° und 2* und nachfolgende Subtraktion ergibt. Unter Berücksichtigung der Gleichung 
der Resultantkomponenten von (4c), der allgemeinen Beziehung der zeitlichen Ableitungen 
des Vektors tr, in 2° und 2*, sowie der bereits verwendeten vektoralgebraischen Identität (15) 
bekommt man schließlich: 

dv,=d,+[wd,]+fp; u]. 


6) Siehe [3], S. 163. 
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In ganz entsprechender Weise bekommt man die Momentengleichungen von (3a) bis (3e). 

Die zunächst nur auf die „Länge“ angewendeten Motorgleichungen hat man in die aus 
v=pd,+ [ru] 

durch Differentiation in verschiedenen Systemen folgenden Beziehungen einzusetzen und unter 
Berücksichtigung der Gleichungen zwischen den zeitlichen Ableitungen von rt’ sowie der 
Gleichung der Resultantkomponenten weiter auszurechnen. Wegen der Beziehung (15) und 
w’||v,, vereinfachen sich die Endformeln wieder erheblich. Wir stellen die Ergebnisse zu- 
sammen: 

Aus den Motorgleichungen (3a) bis (4e) folgen die Vektorgleichungen: 

l. Die Gleichungen der Resultantkomponenten: 


= . A u ner a ee Se 
N See "7, 2. .(16b), 
“=W—-[un]) . :.... Fe 
ee er ei + 
a=8 —[uH]) . ...: 22 2 20. ... (16e), 
eh un] :. 2.2.22 2 2 220 » (i6h 


2. Die Gleichungen der Momentkomponenten: 


% ‚ y 

s=d—-iu®]) ....... |; ° 
*% .r DR ‚ ‚ 

’=v-[Wv]+[Wo)+[Vul. . . een Ü 
v’=v’+[lWvl +{Vul . 2. 2 2 2 2 nenn... (17e), 
0) Pe (17d), 
bed + +WwWoul+ul. . 2 22202020202. (re), 
bed +wWol+lul. . » 22 2 2 222... (M7N. 


Es gilt nun zu entscheiden, welche Momentkomponenten bei der Zerlegung der Gleichungen 
von Abschnitt 2 für die verschiedenen zeitlichen Ableitungen der Motoren ® und & ein- 
zusetzen sind. Dazu beachten wir, daß die Motorgleichungen selbst in allen drei betrachteten 
Systemen 2°, 2° und 2* gelten, da sie ja eine von jeder speziellen Wahl des Koordinaten- 
systems unabhängige Formulierung darstellen. Auch die aus ihnen abzuleitenden Vektor- 
gleichungen sind noch unabhängig von der Walıl der Achsenrichtungen, nicht aber von der 
Wahl des Koordinatenanfangspunktes, da die Momentkomponenten ortsabhängig sind. Auch 
die Vektorgleichungen gelten nach Verfügung über den Koordinatenanfangspunkt — noch 
in jedem beliebigen System. Jedoch bedingt die Verschiedenheit der zeitlichen Ableitungen 
des Vektors r, (bzw. tr’) bei der Zerlegung der Motorgleichungen verschiedene Momentkompo- 
nenten desselben Motors, je nachdem, in welchem System die Zerlegung vorgenommen wird. 
je] 

Zum Beispiel hat 5 in I°* die Momentkomponente d, entsprechend hat $ in 2° die Moment- 
komponente d, und % in 2° die Momentkomponente d. Wird dagegen die Zerlegung eines 
abgeleiteten Motors in einem anderen System vorgenommen als die Ableitung, so folgen die 
Momentkomponenten zwangsläufig aus den zuvor abgeleiteten Beziehungen (17a) bis (17f). 
In Tafel 1 sind diese Verhältnisse für die drei möglichen Zerlegungen nach 2*, 2° und 2° 
dargestellt. Neben den Momentkomponenten selbst sind noch die aus (17a bis f) folgenden 
Zusammenhänge mit den im Zerlegungssystem genommenen Ableitungen desselben Momenten- 
vektors dargestellt. Mit Hilfe der Tafel 1 kann nun die Zerlegung einer beliebigen Motor- 
gleichung — nachdem man sich für ein Koordinatensystem entschieden hat — auf tatsächlich 
nur noch rein schematischem Wege durch Anwendung der Rechenregeln der Motorrechnung 
erfolgen. Gleichzeitig dienen die in Tafel 1 mit angegebenen Zusammenhänge der ver- 
schiedenen zeitlichen Ableitungen dazu, den Wechsel in der Beurteilung der Veränderlichkeit 
der Vektoren noch nachträglich in der Vektorgleichung vorzunehmen. 

Wir wollen nun dazu übergehen, die zuerst abgeleitete Motorgleichung (2) in ihre Vektor- 
komponenten zu zerlegen. Um mit konstantem Trägheitsmotortensor (und damit auch mit 
konstantem Trägheitsvektortensor) rechnen zu können, entscheiden wir uns für eine Zerlegung 
in &* und haben folglich die Vektorkomponenten aus den ersten beiden Spaiten von Tafel 1 
zu verwenden. 

Durch Anwendung der Rechenregeln für das Produkt eines Motortensors mit einem 
Motor’) findet man, daß der Motor des Relativimpulses TW die beiden Vektorkomponenten 


7) Siehe [3], S. 171. 
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T®’ = im (v’ +[u’r,)); m[t,v’]+ Tu’) 
hat. Entsprechend hat man: 
T5 = im(+[ur,])); m[r,v]+ Tu}, 
Td (m +): md] + Tin, 
T5=m6+l[in); mfe5]+ Thy. 
Für die motorischen Produkte bekommt man unter Berücksichtigung der Beziehung: 


[a [6 01] + Is WW) + [0 ain]]=0 
die Vektorkomponenten: 


[V(TB I] = Im d[u’ dv] + [ui [u r,] 
[FITH] = Imifuo] + [ufur,]); mfr,[uo]] + [u(Tu)]}, 
[FTP] = mia] + uf, ] PD: mefufr,v’]] + pv/] + [ofu’r,]) + [ul Tu’)]}, 
[TH] = tm (uw v] + [uf]: mifa’frs0]] + [Po] + [v’ [ur] ) + [u (Tu]%. 
Die Vektorkomponenten des Kraftmotors seien 
BP=t}:M. 


Nun bekommt man durch Zusammenfassung der ersten und der zweiten Vektorkomponenten 





); m[r,[w’v’]] + [w’(Tu’)]%, 








jedes der in (2) vorkommenden Glieder die allgemeinen Vektorgleichungen: 


mid + vr) + [ww] + [fr] Hm + fr] + fuo] + fufur,]]} ) 


‚ ’ 1 ’ ’ N . (18a), 
mtv] + [ur] + [wWfur,)] + [ufu’r]]} = PB | 
im (fr) + [sw v’]) + Ti’ + (Tuw)]) | 
+ {mio + [tsfuo])) + Tü+[u(Tuw)]) . (18b). 


+ tmi[r,[wW0]] + [rs fur] ) + u (Tw] + u (T7W) =-M 


Durch die Zusammenfassung in geschweiften Klammern sind die Anteile der Relativ- und 
der Führungsbewegung sowie der Coriolisanteil hervorgehoben. Ein Vergleich der beiden 
Vektorgleichungen (18) mit der ursprünglichen Motorgleichung (2) läßt die großen Vorteile 
der Motorschreibweise besonders deutlich werden, denn die Verwendung der Motoren erspart 
nicht nur viel Schreibarbeit, sondern gestattet vor allem die allgemeinen Gleichungen in einer 
kaum noch zu übertreffenden Kürze abzuleiten. Man vergleiche damit die direkte Ableitung 
der Vektorgleichungen z. B. bei Heun [2]. 

Wir ersparen uns die Zerlegung auch der anderen Motorgleichungen des Abschnittes 2, 
zumal diese mit Hilfe der Tabelle und nach den Regeln der Motorrechnung leicht durch- 
geführt werden kann. Jedoch wollen wir an Hand des Beispieles von Gl. (185) eine Um- 
formung vornehmen, indem wir die Änderung der Führungsgeschwindigkeit nunmehr von X'° 
aus beurteilen wollen. Dann folgt aus der Tabelle: 


= u-+ [uw], 
v—=d+[up] + [ou]. 
Geht man damit in (18) ein, so fallen verschiedene Glieder heraus, ferner läßt sich 
[u fur,]] +fu[u’r,]] + [[uu’]t;] = 2 [u fr; u’]] 
zusammenfassen, so daß die Vektorgleichungen entstehen: 


md’ + ir) + [Wv]+ Wr] + mio + [ur] + [up] + [ufur,]]} | 


’ ‚ N . . (19a), 
+2mt[fuv’] + [ufu’r,] = PB | 
sm[rsd] + [rsfuv’/]) + Ti’ + [u Tu’)]) 
imifrsv] + [sup] ) + Tu+ [u(Tw)], .(19b). 


omfe, [uw] + u FW] + +) 

Diese Gleichungen sind mit den von Winkelmann und Grammel*) angegebenen voll- 
kommen identisch, wenn man die dort vorgenommene Umformufg des Trägheitstensors im 
Coriolisglied beachtet. 


+ 


3) Siehe [4], S. 448 
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Selbstverständlich kann man (18) und (19) noch weiterhin umformen, wenn man die 
Beurteilung der Geschwindigkeitsänderungen von anderen Systemen aus vornehmen will, wir 
wollen jedoch darauf hier nicht weiter eingehen. Ebenso wollen wir uns nicht mit der all- 
gemeinen Zerlegung der Vektorgleichungen in ihre entsprechenden Skalargleichungen auf- 
halten, da diese Umrechnung in einfacher Weise nach den Regeln der Vektorreehnung er- 
folgt. Wir werden ohnehin bei den später zu besprechenden Beispielen derartige Zerlegungen 
in Spezialfällen durchzuführen haben. 

Im Hinblick auf die vielseitigen Anwendungen soll jedoch an dieser Stelle noch die im- 
plizite Form der Motorgleichung (9) in ihre Vektorkomponenten zerlegt werden. Die beiden 
Vektorkomponenten des Impulsmotors \$ seien ® und \}. Der Resultantvektor © ist der 
Impuls der Fortschreitgeschwindigkeit, der Momentvektor \} ist der Drehimpuls. In den 
vorher benutzten Größen ausgedrückt hat man: 

$=m(v’+[u’r,)). 
y=nml[ı,v]+ Tu. 
Ganz entsprechend seien die Vektorkomponenten des Führungsimpulsmotors $;: ©; und 97. 
| \ 
Für sie gilt: 
$;=m(v+[ur,)), 
Yy=mflıv]+Tu. 
Die Zerlegung von (9) ergibt nun die beiden Vektorgleichungen: 
su +8 +uBl=P. . . . .. + 2 22: (MDR), 


sts) +PpSI +Yy their) t+PSH)I=SM. . . . . . . ..(20b). 


4. Grenzfälle. 

Um uns mit den bisher abgeleiteten Formeln vertrauter zu machen, wollen wir im folgen- 

den einige bekannte Fälle herausschälen, die sich aus den allgemeinen Gleichungen ergeben. 
Zunächst sei der Grenzübergang zum Fall des Massenpunktes betrachtet. Wir wählen den 
Massenpunkt zum Ursprung des 2*-Systems und können folg- 
lich r,=0 setzen. Die Achsenrichtungen von 2°* können denen 
von 2° parallel gelegt werden, dann ist uU’=0. Zur Berech- 
nung genügt in diesem Falle die Gleichung der Resultantkompo- 
nenten z. B. in der Form (19a). Es bleibt wegen 1, = u’ —0: 





N 


mid +0 +[up] +2fuv]}=P. 





Dabei ist d die Geschwindigkeit desjenigen in I° festen 

Punktes, der augenblicklich mit dem Massenpunkt zusammen- 

2 2 ° a Fr 

fällt. Will man stattdessen lieber mit der Geschwindigkeit 5 

b, des Ursprunges 0 von 2° rechnen, so hat man zu zetzen: 3ild 2. 


vb=d,-+[ur]. 
Darin ist r der Vektor von O nach m (Bild 2). Durch Differentiation in 2° folgt daraus: 
db=bd,-+[ur], 


denn t=0, da r in 2° fest ist. Durch Einsetzen folgt nun: 


md +v+[un)]+für+fufur]] +2 =P. . ..... 021). 
Darin kann auch an Stelle von d’ nun vd’ geschrieben werden, da wegen W=0 d’ —Y’ ist. 
Dann werden also sämtliche Ableitungen in 2° genommen, so daß man diese Vektorgleichung 


bei Zerlegung in Skalarkomponenten am vorteilhaftesten nach &° zerlegt. (21) ist die be- 
kannte allgemeinste Form der Bewegungsgleichung für die Relativbewegung eines Massen- 
punktes, wie sie z. B. bei Grammel?) zu finden ist. ; 

Durch Grenzübergang zu einem in I° festen Körper (®’—0) kommt man aus (2) zu 


der Motorgleichung: 


TF+HBTB=P, 
mid+ [ur] +fup])+[ufur] =P . : : : 2 20202 .(228), 


mtftsol)+ [fu] + Tutlu(TW)=-M .. . . 0. 0...(22b) 


®?) Grammel: Kinetik der Massenpunkte. Handb. d. Physik. Bd. V, Berlin 1927 Kap. 7, S. 330, 


deren Vektorzerleguug 
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ergibt. Das sind die ebenfalls bekannten Gleichungen für die allgemeine Bewegung eines 
Körpers in einem Inertialsystem, wie sie z. B. bei Winkelmann und Grammel'’) zu finden 
sind. Dieselben Gleichungen bekommt man natürlich auch durch den Grenzübergang zu 
einem in 2" festen Bezugssystem I° (5 =0). nur daß dann gestrichene Größen an Stelle der 
ungestrichenen in (21) treten. 

Aus (22) lassen sich weiterhin die bekannten Eulerschen Gleichungen für die Bewegung 
des starren Körpers ableiten, wenn man entweder den Körper in einem Punkte festhält und 
diesen Punkt zum Koordinatenursprung wählt (vd —0) oder wenn man den Ursprung in den 
Schwerpunkt legt (r,=0). In beiden Fällen bleibt von (22b) die Vektorform der Euler- 
Gleichungen: 

Tu+[u(Tw]=M. 


5. Anwendung auf das physische Pendel von einem Freiheitsgrad. 


Als einfaches Anwendungsbeispiel sollen nun die strengen Bewegungsgleichungen für 
ein Pendel aufgestellt werden, das um eine in I° feste Achse drelibar gelagert ist. Wir 
wählen den Aufhängepunkt O0 des Pendels zum Ursprung der 
Systeme 2* und F° und legen die z-Achse von 2° sowie die 
(-Achse von 2* so, daß sie in die Richtung der Pendelachse 
fallen. Die A-Achse von 2* legen wir senkrecht zur (-Achse in 
die durch die Pendeldrehachse und den Pendelschwerpunkt be- 
stimmte Ebene. Die B-Achse wird schließlich senkrecht zu den 
Achsen A und C gewählt. 

Als Ausgangseleichungen der weiteren Rechnung nehmen 
wir die Vektorgleiehungen (19), brauchen uns jedoch um die 
Gl. (19a) nicht zu kümmern, da sie nur die hier nicht weiter 
interessierenden Lagerkräfte ergibt. Auch die A- und B-Kompo- 
nenten der Vektorgleichung (19b) können unberücksichtigt bleiben. 
da sie ebenfalls nur Lagermomente bestimmen. 

Bei der ©-Komponente fallen einige Glieder heraus, nämlich 
diejenigen Vektorprodukte, deren einer Faktor uw’ ist. Da u’ stets 

3ild 3. in der C-Richtung liegt, haben die mit u” gebildeten Vektor- 
produkte keine Komponente in der C-Richtung. 

Auch die beiden letzten Glieder der linken Seite von (19b) fallen heraus, da sie sich 
gegenseitig aufheben. Man bekommt schließlich eine Gleichung von der Form: 





Kür + ma (vg + Wera uite]— PrlWd)ıa — Pi) + Iclu) 
+ uy! D, nr Op Up PD, ur) up\ wir Ur PD, Up PD; Ur) M, . 


Darin ist z. B. (v), die in Richtung der B-Achse genommene Komponente des in 2° abgeleiteten 


Vektors d. Zwischen den nach 2° und den nach I* genommenen Komponenten bestehen die 
folgenden Transformationsgleichungen, bei denen zur Abkürzung \u-di=g gesetzt wurde: 


/ 
u=u,cosp—+tu,sing, 

NYg=—u,sinp + u,cosq, 

un = U, 


entsprechend: 
(= u,cosp +tu,sing, 


(u): u.sing + u,cosp, 
tle=Wu,, 
(v)p v.sin g u BER Uy COS p. 


Auch in dem äußeren Moment Mc kommt die Unbekannte 9 bzw. u. vor. Sei das 
Moment z. B. verursacht durch ein Schweremoment und durch eine von u abhängige Dämpfung, 
so hat man, wenn die Komponenten der Erdbeschleunigung in der x- bzw. y-Richtung g; 
bzw. g, sind: 

M=— mg,asing +mg,acosy — flur). 


Durch Einsetzen dieses Wertes und der angegebenen Transformationsformeln bekommt man 
unter Berücksichtigung trigonometrischer Umrechnungsformeln die Pendelgleichung: 


10) Siehe [4], S. 380, 








a 
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che + fliup)+ Esing + Feosp + Heos2p + Ksin2g GR - 
mit den Abkürzungen: 


E= ma(v, uU,0,y+ Uyvd,) PD ,lu;u, Ux) Pplu. uU. ru) tt Mmgxd, 


F=ma(lv,+ u. — ur.) — Pylu.u. tu, + Prlu; u, — ur) — m gya, 


H= PDelu), u) + (Ip I ut, 


K=—-20deu,u,+ . (91 - MW Wi). 
Darin sind die Komponenten der Bewegung von 2° gegen I’ also u. u,u,,vxr,r, sowie 
deren Ableituneen in £° bekannte Funktionen der Zeit, also sind die EFHK als Funktionen 
von t gegeben. Ebenso ist auch f(«.) als Funktion von ur. gegeben. 

Die Bewegungsgleichung (23) zeigt augenfällig, wie kompliziert die allgemeinen 
Gleichungen selbst in so einfachen Fällen wie dem vorliegenden werden können, wenn die 
Bewegung von einem beliebig bewegten System aus beurteilt wird. 

Wir wollen noch eine naheliegende Spezialisierung der Gl. (23) vornehmen und sie auf 
ein auf der drehenden Erde mit horizontaler Achse schwingendes Pendel anwenden. Dann 
können u und v als konstant angenommen werden. Ferner wollen wir uns auf eine Be- 
trachtung kleiner Auslenkungen des Pendels aus der Vertikallage beschränken und folglich 
die Gleichungen für kleine Winkel 9 linearisieren. Man erhält so die Gleichnng: 


%#+(E+2K)P=—(F+H), 


deren Koeffizienten nunmehr konstant sind. Die Lösung ergibt harmonische Schwingungen 
von der Schwingungszeit: 


Or 
n__ 9) ( 
En EL>OK' 
um eine Gleichgewichtslage: 
F+H 


( == - >. 
ro E+2KR 
Um noch einige Vereinfachungen in den Konstanten zu bekommen, wollen wir ferner dahin- 
gehend spezialisieren, daß wir 2* als Hauptachsensystem voraussetzen, also P,= P, =D =0 
setzen. Außerdem wollen wir üblicherweise die von der Erdanziehung herrührende Kompo- 
nente der Beschleunigung mit der durch die Erddrehung verursachten Zentrifugalbeschleunigung 
zusammenfassen: 
7 Ix Uy U; + U; Uy 2 j S, ’ 
Es bleibt nur diese Komponente übrig, da die y- und 2-Komponente wegen der Wahl von 2° 
verschwinden (z-Achse vertikal, y- und z-Achse horizontal). 
Bezeichnen wir ferner den Winkel zwischen der Horizontalebene aın Beobachtungsort 
und dem Drehvektor u mit 9, so kann man für die Komponenten von u schreiben: 
UY=—ucosV, 
uw,=usin®sina, 
u. —=usind cosa. 
Dabei ist « der Winkel zwischen der Nordrichtung und der positiven z-Richtung. Mit Ein- 
setzen der angeführten Werte bekommt man nun für Gleichgewichtslage und Schwingungs- 
zeit die folgenden Werte: 
(H-9)Ww 


o,=+ —iudemnma .. :.» . 3 ; 2 

Fo T Im a q n K r - . ; ‚ : (24), 

Ton w/R | (Ip -— 9 ,) u (cos? — sin? 9 sin’a)] (25) 
2 OR mga, 2mayg' | . 


Danach zeigt ein physisches Pendel von einem Freiheitsgrad bei horizontaler Aufhängeachse 
nur dann die genaue Vertikale an, wenn entweder 9 ,—= 9, oder a = O (Pendeldrehachse in 
Nord-Süd-Richtung) oder 9—=0 bzw. W° (Beobachtungsort am Äquator oder am Pol) gilt. Die 
Abweichung ist jedoch sehr klein, da sie mit dem Quadrat der kleinen Größe u eingeht. 
Beispielsweise beträgt die Abweichung eines Stabpendels (9 ,=0) von einem halben Meter 


+ I 


i ® } u 2 y - x . 
reduzierter Pendellänge a 0,5 m} mit Ost-West gerichteter Pendelachse (a —= 9%) für eine 


geographische Breite von #=45," nur 


9% = —2,710% Bogensekunden. 
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Der Wert ist so klein, daß er sich bei dem heutigen Stand der Meßtechnik jeder Beobachtung 
entzieht. 

Gleichermaßen ist die Änderung der Schwingungszeit so gering, daß sie mit den heutigen 
Hilfsmitteln nicht beobachtet werden kann. Das Zusatzglied von (25) beträgt maximal für 
das vorhin angegebene Pendel 1,3 10°'". Am Pol wird die Schwingungszeit um diesen Betrag 
verlängert, am Äquator wird sie in demselben Maße verkürzt, sofern die Pendelachse in der 
Öst-West-Richtung liegt. Schwingt das Pendel dagegen um den Nord-Süd-Meridian, so fällt 
das Zusatzglied fort. 


6. Systeme starrer Körper. 


Die allgemeinen Motorgleichungen des Abschnittes 2 lassen sich ohne Schwierigkeiten 
auf Systeme starrer Körper übertragen. Werden die Bewegungen mehrerer starrer Körper 
vom Führungssystem aus betrachtet, so hat man für alle den gleichen Motor & der Führungs- 
geschwindigkeit, aber verschiedene Motoren ®; der Relativgeschwindigkeit, sowie verschiedene 
Trägheitsmotortensoren T; anzusetzen. Außerdem hat man zum System der äußeren Kräfte 
(Motoren ®;) die zwischen den einzelnen Körpern wirkenden Reaktionskräfte und -momente 
hinzuzufügen, deren Motoren R; ; seien. Durch Summation über sämtliche Körper erhält man 
dann aus (5) die Gleichungen: 


En o° 
TEHRTBHTEH BT) | 
(26). 
; (i De: 
Die Matrix der Reaktionsmotoren R;; ist dabei schiefsymmetrisch, d. h. es ist stets R;; = 0 


und R;= — R;; (actio—reactio). Zwangsläufige Bewegungen einzelner Körper werden sich 
dabei in Zusatzbedingungen zwischen den Bewegungsmotoren ®; bzw. deren Komponenten aus- 
drücken. Wenn z. B. zwei Körper in einem Punkte fest miteinander verbunden sind, so gilt, 
falls dieser Punkt zum Ursprung der körperfesten Systeme gewählt wird, stets v/;—=d;,;, die 
Vektoren der Drehgeschwindigkeiten sind dagegen im allgemeinen voneinander verschieden 
w;=Fu;. 

Mit Rücksicht auf spätere Anwendungen wollen wir noch die Motorgleichung in der 
Form (9) für Systeme angeben. Wir erhalten: 


je} je} v 

G+HBII HI +HEHI-B- I R;=V0 . 2.2.2... 
a Ne 
7. Anwendung auf ein Drei-Körper-System. 


Kardanisch gelagerte Kreisel: werden bei fast allen Kreiselgeräten sowie auch bei 
Kreiselmodellen für den Unterricht verwendet. Bei der Berechnung der Geräte werden die 
Massen der Kardanringe jedoch meist ganz vernachlässigt oder nur in Näherungen berück- 
sichtigt. Wir wollen uns daher die Aufgabe stellen, mit Hilfe der Gl. (26) die Bewegungs- 
gleichungen für einen aus einem Kreisel (1), einem inneren Kardanring (2) und einem äußeren 
Kardanring (3) bestehenden Verband von drei Körpern in völliger Strenge anzuschreiben, eine 
Aufgabe, die meines Wissens bisher auch für nicht be- 





wegtes Bezugssystem nicht behandelt worden ist. # 
Wir legen neben dem Bezugssystem I° drei ver- m 

schiedene körperfeste I (Kreisel), 25 (innerer Kardan- /f NN 

ring) und 2% (äußerer Kardanring) derart in den Ver- \ 


band, daß die A-Achse von 2 zusammen mit der -Achse 
von S° in die in I'° feste Drehachse des äußeren Kardan- 
ringes fällt. Die B-Achsen von 25 und 27 sollen in die 
Richtung der inneren Kardanachse fallen, schließlich 
sollen die C-Achsen von 2% und 2% in die Figuren- 
achse des Kreisels fallen. Nach Bild 4 sollen die Kardan- 
ringe durch Federn derart an das Bezugssystem gefesselt 
sein, daß im Gleichgewichtsfall die innere Kardanachse 
mit der „Achse von £°, die Figurenachse des Kreisels 
mit der z-Achse von I° zusammenfällt. Bei Gleichgewicht 
fallen also dann die A- und B-Achsen von 25 und 5 
mit der x- und y-Achse von 2° zusammen, ferner gehen Bild 4. 


dann alle drei körperfesten C-Achsen in die z- Achse. 
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Zur Vereinfachung der Rechnungen wollen wir voraussetzen, daß die vorstehend ge- 
wählten körperfesten Systeme 2%, 2% und 2% Hauptachsensysteme sind, und daß ferner der 
Schwerpunkt S aller drei Körper in den Schnittpunkt der Figurenachse mit Jen Kardan- 
achsen fällt. D: k: ler Schwerpunkt als in I°° fester Punkt : Ursprung der Systeme 
achsen fällt. Dann kann der Schwerpunkt als in I°° fester Punkt zum Ursprung der Systeme 
2°, 5*, $* und I# gewählt werden. Damit vereinfachen sich die Ausgangsgleichungen ganz 


erheblich. Wir brauchen nur die Momentkomponenten der einzelnen Glieder von (26) zu be- 
rücksichtieen und bekommen wegen t,=0 die Bewegungsgleichungen in der Vektorform: 


TE + (T;w)]+ Tu + [u(T;uw)] + [ui (T;u)] 
+ [u(T;w)] + (T;fuui)) -M;—- VRR; = 0 i=1,2,3) 
I 
Darin ist: 


| 94, 0 0 | 

4 0 Op, 0 (genommen nach IF), 
| 0.00 
| O4, 0 0 

A 0 Op, 0 (genommen nach 2%), 
| ) 0 Se, | 
| 9,0 0 

A 0 Ip, 0 (genommen nach 27) 
| 0 0 Ic, | 

M, — 0 0) 

ne ® , 0, 0) (nach 2%). 

Mn. —=(Mı, Mir, Mc;,) 


Dabei ist M.,, das von den zwischen äußerem Kardanring und Gehäuse ausgespannten Federn 
ausgeübte Moment, Mg, und Mc, sind die Lagermomente. 
Die Matrix der Reaktionskräfte hat das Schema: 
ie RR; 0 
sy x 





t;; N, 0 Rt, 
lo n,o) 
mt Re NR. Rs =—NR,,. Ferner R,=R,,=0, da der Kreisel und der äußere Kardan- 
ring nieht aufeinander einwirken. Für die Komponenten hat man: 
Rule Zn 8 7 (nach 2%), 
n.=(Mi Main Mca) (nach 33). 


Die Vektoren der Drehgeschwindigkeiten haben die Komponenten: 


“=(u4, wB, u) (nach 2*), 

h wuelns, Wu, We (nach F%) 

wies, ©; 0) (nach 2%), 

i u =(u4, uB, uc ) (nach I) 

(ü4,, WB;; C,) (nach 27), 

; (nä.,; MB, Uc,) (nach 2%). 
h Bevor wir nun diese Werte in die Gl. (28) einsetzen, müssen wir noch die kinematischen 


Beziehungen zwischen den in den verschiedenen Systemen genommenen Komponenten auf- 
stellen. Jede beliebige Lage des betrachteten Systems läßt sich durch drei nacheinander 
ausgeführte Drehungen um die Achsen 

x — A,-Achse (äußere Kardanachse), 
B,= B,-Achse (innere Kardanachse), 
,— (,-Achse (Kreisel-Figurenachse) 


2 


erreichen. Wir bezeichnen die Winkel der um diese Achsen ausgeführten Drehungen mit 
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und können folglich die Lage des Systems durch diese drei Winkel angeben. 
Die drei Winkel stimmen nicht mit den meist verwendeten Eulerschen Winkeln über- 
ein, sind Jedoch dem vorliegenden Fall besonders angepaßt, da sie am Modell unmittelbar ge- 


messen werden können. 

Jeder der genannten Drehungen entspricht eine orthogonale Transformation, Werden 
die Drehungen hintereinander ausgeführt, so entsteht wieder eine Transformation, deren 
Matrix durch Multiplikation der Matrizen der Einzeltransformationen entsteht. Wir werden 
daher zuerst diese Einzeltransformationen bestimmen. 

l. Drehung um die äußere Kardanachse (.- bzw. 

Das entspricht der Transformation der Komponenten zwischen 


mit einem beliebigen Vektor p: 


A,-Achse), siehe Bild 5. 
Yo und &#, Wir schreiben 


(29) 


ns 


(Po —T,(p)ss 


oder in Komponenten : 
Px — P4s 


Pu = PB; COS a PC3 sin «a, 


Pz = PB, Sina pc, Cosa. 


Die Transformationsmatrix ist also: 


| 1 0 0 
I 0 cos a sin a ; 
| 0 sin a cos a | 


2. Drehung um die innere Kardanachse (B,- bzw. B,-Achse) ergibt in analoger Weise 


(siehe Bild 6): 
13 


nf 
Be. - 
5%) 











2:4; A ——n 
ee 
8; ” E 
. & 7 / 
Bi kr, 
/ 
A404 5 45 Er 2 / y‘ 
Bild 7. 


Bild 


Bild 5. 
Ms;=?T,(ps; 


mit der Matrix 
| cosß 0  sinp | 


„= 0 1 V 
| sin? 0 cosp | 


3. Drehung um die Figurenachse (C,- bzw. C,- Achse) ergibt (siehe Bild 7): 
EEE EEE 


PMs=T(p)s; 


ti 


mit der Matrix: 
cos gq sin g 0 | 
%,=1sing cos gq 8. 
| 0 0 | | 


Bei der weiteren Rechnnng werden nun noch die Transformationen zwischen den Systemen 
2% und 3, 2° und 2% sowie zwischen 2° und 2F benötigt. Die Matrizen dieser Trans- 


formationen seien: 


T,=>TT, (27 nach 2%), 
22%, (2° nach 23), 
2.= 2, 2,%, (2° nach 2%). 








XUM 
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Nach den Regeln der Matrizenmultiplikation bekommt man: 


| cos P coSg cos $ sin y sin ß | 
T, sing cos y 1) a ee a 
| sin P cos q sin sing cos | 
| cos ß 0 sin ß | 
= sinasin  cosa sinacosß |. . (33) 
| cos a sin p sina COS 4 COS P 
| COS PCOSg cos sing sin p 
T, sinasinPcosp + cosasing sinasin?sing + cosa cosy sinacosßt. . . (34. 
| cosasinPßcosy + sinasing cosasindsing—+Sinacosg ae 


Nach diesen Vorbereitungen kann das System (28) in seine Komponentengleichungen 
zerlegt werden. Im Falle i=1, also für den Kreisel selbst sollen alle Komponenten nach 2 
genommen werden. Zu diesem Zwecke muß die im dritten Ausdruck von (28) vorkommende 
Ableitung von u nach 2° vermöge der Transformation (34) nach 27 transformiert werden. 
Wir schreiben für die Komponenten in 27 in unmittelbar verständlicher Weise: (u).4, (u): 
(w)e- In gleicher Weise müssen die nach 2% genommenen Komponenten von R,=—NR,, 
vermöge (31) nach 27 transformiert werden. Man bekommt auf diese Weise die drei Kom- 
ponentengleichungen: 

94, [a + (4 + (Ic, — IB.) (up uc+upüuc+ wceug + up ur) 
+ 94, (upue - tet) = M4,€0sy—+ MR, Sing 


Op, [tn + Wp] + (91, — 9e,) (Wauc+usuc+tusurtucn)) (35) 


+ Ip, (ve uw, —uıur) M..,sing + Mp,,cosgq 


Ar, ne + (We) + (OR, — 941) (upuatupustupu,+uuu,) 





+9, (us up up u) 
Im Falle i=2, also für den inneren Kardanring sollen alle Komponenten nach 27 ge- 
nommen werden. Man hat dann die Komponenten von u vermöge (33) und die Komponenten 
von R,;: R,, vermöge (30) nach 2# zu transformieren. Die Ausrechnung ergibt dann die 
gleichen Ausdrücke wie im System (35), nur überall mit dem Index 2 an den Trägheits- 
momenten und Komponenten der Drehung. Die Momente sind dagegen etwas anders. Man 
bekommt: 
94,4 + (W).1,] + (ec, — 98,) (uR, uc, + up, uc, + uc, Up; + BR, UC,) 


+ 94, (UB, WC, — UC, up.) = Ma, + Ma,, cos # — Mec,, sin p 
Op, [tr + (WB,] + (94, —- 9e,) (u, ue, + U, uc, + u, uc, + uc, u)4,) (36) 
U), 
+ Op; (uc, u, U, Uc,) = MR, + MB; 


9; [uc; + (w)e,] + (Op, — 94,) (WB: Wi + Up, UA, + Up; WA; + UA, UB,) 





+ Ic, (1.4, UB, — UB, u4,) = M a,, sin + Me,, cos 





t Schließlich werden im Falle i=3, also für den äußeren Kardanring alle Komponenten 

t nach 2# genommen. 

h Dann hat man vermöge (29) die Komponenten von u nach 2% zu transformieren. Man 

bekommt: 

i 91,[u4 + (u).1,]+ (9, — 98,) uc, ur = Ma, + Ma;: 

; Ip; w)B, Tr I 43 IKc;) U, + ( O4; + OB; 7 3) U4 UC, = M»L + MB;; | . 8. 
Ic; \ u)C, +(Op; - 94) UB;, UA, — (Os OB; + 9c;) UAs u, = Mc; + Me 32 


Dabei bestehen zwischen den Komponenten der Eigendrehung der drei Körper noch 
vier Beziehungen, die aus den Transformationen (31) und (32) folgen: 


un u 1cosy —- upsing 


a TEEN 


vB, =u4sing—+ uprcosg EB Midland Loc SE 7.2 Se 
uc, =uc 
"= u4cosßcospy— upcosßsinpg+tuesinßp . 2. 22.202... 











XUM 
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Man hat somit in den neun Bewegungsgleichungen (35), (36) und (37) und den vier 
kinematischen G]. (35) und (39) insgesamt dreizehn Gleichungen für die dreizehn Unbekannten 
4, UB, UC, Ur, UBz, UC,, WAz, MAsı; MBaı, MAss, Mc;,, MB;, Mc;. 

Durch schrittweise Elimination der Lagerkräfte sowie der Drehungskomponenten u4,, 
UB,, uc, und u4, könnte man das Problem auf insgesamt drei Gleichungen zurückführen, 
entsprechend den drei Freiheitsgraden des betrachteten Systems. 

Die Ausrechnung wird jedoch sehr bald so undurchsichtig, daß wir im allgemeinen Falle 
darauf verzichten wollen. Selbst in dem ganz einfachen Sonderfall eines z. B. nur um die 
z-Achse von X konstant drehenden Bezugssystems vereinfachen sich die Formeln nicht 
wesentlich, da die Komponenten von u in den drei körperfesten Systemen im allgemeinen 
von Null verschieden sind, also nur wenige Glieder fortfallen. 

Wir werden uns deshalb im folgenden darauf beschränken, die Bewegungsgleichungen 
für den in unseren Betrachtungen mit enthaltenen Sonderfall anzugeben, daß sich der 
kardanisch gelagerte Kreisel in einem Inertialsystem befindet, also u =0 wird. 


Dann folgt aus der zweiten Gleichung von (37) sofort das Lagermoment Ma; MB3s>. 

Eine Unbekannte ist damit ausgeschieden. Aus der ersten und dritten Gleichung folgt: 
Ma.» = 04, 44, — Ma,, 
Mc„=— Mcı;- 
Diese Werte hat man in (36) einzusetzen, wobei zu berücksichtigen ist, daß M;; M;; ist. 
Neben den bekannten von den Federn ausgeübten und also von den Winkeln « und 9 ab- 
hängigen Momenten M4, und Mp,, kommen jetzt noch die Momente M4,,, Mp,, und Me, 
vor. Die dritte Gleichung von (36) enthält nur noch die Unbekannte Mc,, die damit er- 
mittelt ist, und deren Wert in die erste Gleichung von (36) einzusetzen ist. Aus der ersten 
und zweiten Gleichung von (36) lassen sich schließlich die Momente M,,, und Mr,, aus- 
rechnen und in das System (35) einsetzen. Damit ist dann die Elimination der im allgemeinen 
nicht interessierenden Lagerkräfte durchgeführt. Wir wollen jedoch das Ergebnis in der bis- 
herigen Bezeichnungsweise nicht erst angeben, da sich die Gleichungen nach Einführung der 
Impulse wesentlich übersichtlicher fassen lassen. Wir führen ein: 
%, = Impuls des Kreisels, 

\$> = Impuls des inneren Kardanringes, 
\; = Impuls des äußeren Kardanringes. 
Unter Berücksichtigung der auch für die Komponenten dieser Impulse geltenden Transfor- 
mationsformeln (29) bis (34) bekommt man durch den angegebenen Eliminationsprozeß aus 
(35) die Gleichungen: 





Sı)at+&Heur (Sı)zuc={M a Solar [+ Se. + pr, u, (Ie)azun.|sin’ßy 
n S/ 
+ MB; sing (Se)ı (Ile (u)g + Vr)r ur, 2. 
Sı)a Fr (1)a u — (Sc Pr (Ma, (K)a, Kar + (2)B> Udo (N2).As uB,] sin? P} (40). 
+ Mp,, €0S 4 (Ion (2); (U;)e + (Nele Was 
(Fi) + (lau —Aı)aup=V 


Dies sind also die vollkommen strengen Bewegungsgleichungen für einen kardanisch 
aufgehängten Kreisel. Sie reduzieren sich im Falle verschwindender Masse der Kardanringe 
auf die bekannte Eulersche Form. Es mag betont werden, daß selbstverständlich auch die 
kinematischen Eigentümlichkeiten der kardanischen Hängung in diesen Gleichungen wieder- 
gegeben werden. So deutet das Glied mit cos ö# im Nenner an, daß nach einer Drehung des 
inneren Kardanringes um 9° die Gleichungen ihren Sinn verlieren, weil dann die Zwangs- 
kräfte (nämlich die Lagerkraft M +.) gegen unendlich geht. In der Tat verliert der Kreisel 
dann einen seiner Freiheitsgrade. 


Auch der sogenannte kinematische Fehler kardanisch aufgehängter Kurskreisel (kräfte- 
freier, als Richtungshalter verwendeter Kreisel), der darin besteht, daß der äußere Kardan- 
rahmen (Kursrahmen) bei bestimmten Drehungen des Bezugssystems mitdreht, ist in den all- 
gemeinen Gl. (36) bis (39) enthalten. 
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8. Die Näherungen. 

Die strengen, im körperfesten System angeschriebenen Gleichungen bereiten große 
Integrationsschwierigkeiten, so daß ihnen wohl keine besondere praktische Bedeutung zu- 
kommt. Ahnlich, wie schon die einfachsten Formen dieser Gleichungen, die Euler-Gleiehungen, 
praktisch wenig Bedeutung erlangt haben, da ihre Integration bisher nur für wenige Sonder- 
fälle, z. B. den kräftefreien Kreisel, den schweren symmetrischen Kreisel und den Kowa- 
lewskischen Kreisel gelungen ist. Das Beispiel des vorigen Abschnittes führte die Schwer- 
fälligkeit der körperfesten Rechnung besonders deutlich vor Augen. Die Gleichungen mußten 
dort in drei verschiedenen Systemen angeschrieben werden. Ihre strenge Integration bringt 
selbst in einfachsten Fällen solche mathematischen Schwierigkeiten mit sich, daß die Lösung 
hoffnungslos erscheint. Meines Wissens sind die strengen Bewegungsgleiehungen für einen 
starren Körper auf bewegtem Bezugssystem bisher in keinem einzigen Falle wirklich gelöst 
worden. Dagegen gibt es eine Anzahl von Fällen, in denen Näherungslösungen angegeben 
worden sind, die sich aus den strengen Gleichungen durch Anwendung der Methoden der 
„kleinen Schwingungen“ herleiten lassen, z. B. die Theorie des Kamerlingh-OÖnnesschen 
Pendels bei Winkelmann und Grammel'). Dabei ist es jedoch meist zweckmäßiger, 
von den in 2° angeschriebenen Gleichungen auszugehen und nicht erst die körperfesten 
Gleichungen als Grundlage zu nehmen. Zwei Beispiele dieser Art sollen im folgenden den 
dabei einzuschlagenden Weg aufzeigen. 

Die Schwierigkeit des Anschreibens der allgemeinen Motorgleichungen in 2° besteht 

wie bereits anfangs erwähnt wurde darin, daß der Trägheitstensor zeitlich veränder- 
liche Anteile besitzt. Man kann nun entweder so vorgehen, daß man den Trägheitstensor 
in I'° ausrechnet, oder daß man auf die explizite Form der Ausgangsgleichung verzichtet 
und statt dessen mit den Gl. (9) bzw. (27) rechnet, in denen nur die Impulse auftreten. 
Dieses letztere Verfahren erweist sich insbesondere bei der Berechnung schneller Kreisel als 
sehr fruchtbar und führt zu Bewegungsgleichungen von großer praktischer Bedeutung. 

Der Vollständigkeit halber sei erwähnt, daß sich die allgemeinen Gleichnngen voll- 
kommen streng für den Fall eines im Schwerpunkt aufgehängten Körpers mit kugelförmigen 
Trägheitsellipsoid sowohl in I’ als auch in F° sowie in 2* anschreiben lassen. Der Trägheits- 
tensor 7 reduziert sich dann einfach auf die skalare Zahl ©. 

a) Kleine Schwingungen pendelnder Körper. Zunächst sei der Fall betrachtet, 
daß der Körper kleine Schwingungen um eine in 2° feste Gleichgewichtslage ausführt, wie 
dies z. B. bei den Schwingungen eines elastisch gelagerten Motorblockes auf bewegtem Flug- 
zeug der Fall ist. Die die Schwingungen beschreibenden Variablen, d. h. die drei Verschie- 
bungen längs der «-, y- und z-Achsen sowie die drei Drehwinkel a, ß, y um diese Achsen 
werden als klein von erster Ordnung, angenommen. Die Systeme 2° und 2* seien so ge- 
wählt, daß sie im Gleichgewichtsfall zusammenfallen, als Ursprung wird zweckmäßigerweise 
der im Gleichgewichtsfall mit dem Schwerpunkt zusammenfallende Punkt gewählt. Die Auf- 
gabe besteht nun darin, den Trägheitstensor für den gegen 2° etwas verschobenen und ver- 
drehten Körper als Funktion dieser Ortskoordinaten in 2° aufzustellen. Dabei sollen alle 
Größen zweiter Ordnung (also Quadrate und gegenseitige Produkte der Verschiebungen und 
Verdrehungen) vernachlässigt werden. Wir können daher von vornherein die durch Ver- 
schiebungen verursachten Änderungen der Trägheitsmomente vernachlässigen, da sie auf Grund 
des bekannten Steinerschen Satzes mit dem Quadrat dieser Verschiebungen eingehen. Auch 
die Deviationsmomente ® sind bei Verschiebungen von zweiter Ordnung klein, da sie mit 
dem Produkt je zweier Verschiebungen gehen. Es genügt also Verdrehungen zu betrachten, 
die wir uns durch die kleinen Winkel «a, f, y gegeben denken. Wird F* gegen 2° verdreht, 


l 5 
so bestehen zwischen den nach 2° und den nach 2* genommenen Impulskomponenten die 
FF -Jpy+JcP | 


Beziehungen: m. 
- Du 
R 4h,= Juartda Kr) ::. . . 
Fe | 


C 


Jıß+Jpa+Jc 
| Diese lassen sich entweder unmittelbar aus Bild 8 ablesen oder 
| auch aus der Transformationsmatrix (34) herleiten, wenn man dort 
die Winkel a, $ und y=g als klein annimmt, also sina=a, 
cos a—=1 setzt usw. Für die Impulskomponenten nach 2* hat 


man zu setzen: 





I un, — Peup— Pr. 
BE 5, l < 
va A»=-Deu, + hup— Pıuctk: - - - (42). 
Zild 8. ‚ ‚ ; 
zo Je: DPrwWs- DPyupr+9ec te 


1!) Siehe [4], S. 449 
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Ferner gilt zwischen den Komponenten der Dreligeschwindigkeit die gleiche Beziehung (41), 
wie zwischen den Impulskomponenten. Wir schreiben sie diesmal in der reziproken Form: 


> 


Tr u. tTuy— aß 
2 Kam. 55 din a6 "U, „75 2 ME. ' 
, u.P —Uuya+tı;, 


Setzt man nun (43) in (42) und danach (42) in (41) ein, so bekommt man die Impulskompo- 
nenten in 2° als Funktionen der Drehgeschwindigkeitskomponenten in 2°. Wegen der 
Beziehung ‘$= Tu kann man daraus die Glieder des Trägheitstensors T in 2° bestimmen 
und bekommt bei Vernachlässigung der von zweiter Ordnung kleinen Glieder: 


| O,+2P, y 20, pP PD, +9 4 - Ip)Y - P,P+ D,a D5 + OP Dea+ Dar) 
- DPe+9a- Ip)y + Ppa- Daß Op+2P ,a-20ey Pu+(9p- I )a+PePp- Darf 
| PD, (I-9)P+ Dıy Dr a PD, (IM-Ke) a+ Dep Dr y Or +2 DaP 20, a | 


Würde man darin auch die von erster Ordnung kleinen Glieder vernachlässigen, so be- 
käme man natürlich genau die Trägheitsmatrix für das System 2*. Ist I°* Hauptachsen- 
system, so vereinfacht sich der Trägheitstensor zu: 


vr (94—- Myr I% O)P | 
T=1(0,—- 9,)y 9% (Hp-9)ar. 
(Ic—-9ı)P (IM—-I)a I. 





Mit Hilfe dieser Matrix lassen sich nun die Vektorgleichungen (15) oder (19) nach 2° zerlegen, 


und damit die gesuchten Näherungsgleichungen gewinnen. Dabei treten bei der Bildung des 
Impulses ‚= Tu Glieder von der Form w;a, «;ß, u;y usw. auf, also nichtlineare Ausdrücke. 
Man wird sich daher in jedem einzelnen Falle Rechenschaft über die Größenordnung der 
einzelnen Glieder ablegen müssen, um zu sehen, ob sie als Größen erster oder zweiter Ordnung 
angesehen werden können. 

In den allermeisten Fällen sind die Winkel harmonische, oder aus harmonischen An- 
teilen zusammengesetzte Funktionen der Zeit mit endlichen Frequenzen. In diesen Fällen 
sind mit den Winkeln auch deren Ableitungen, die Drehgeschwindigkeitskomponenten als 
klein zu betrachten, so daß die Produkte von Winkel und Drehgeschwindigkeitskomponente 
vernachlässigt werden können. Man kommt so zu Endgleichungen, die man auch durch Ver- 
nachlässigung der von erster Ordnung kleinen Glieder im Trägheitstensor erhalten kann. Es 
ist deshalb zweckmäßig, in allen Fällen, in denen keine abnorm hohen Frequenzen zu er- 
warten sind, von vornherein den Trägheitstensor in 2° gleich dem in 2* anzusetzen, also 
die von erster Ordnung kleinen Glieder zu vernachlässigen. 

Auf diese Weise erhält man zwar Gleichungen, die sich von den raumfest angeschrie- 
benen Gleichungen formal nieht unterscheiden, die jedoch gänzlich andere Bedeutung haben 
und nur für „kleine Schwingungen“ gültig sind. Es mag auch an dieser Stelle betont werden, 
daß bei der Anwendung derartiger Näherungsmethoden stets eine besondere Vorsicht angebracht 
ist, da die Lösungen der angenäherten Gleichungen nur für „hinreichend kleine Zeiten“ eine 
Näherung für die wirklichen Verhältnisse darstellen. 


b) KleineSchwingungen von Kreiseln und Kreiselverbänden. Bei Kreiseln, 
also Körpern, die auch im Gleichgewichtsfall noch eine Eigendrehung besitzen, führt eine 
ganz ähnliche Rechnung zum Ziel. Wir nehmen an, daß der Körper im Gleichgewichtsfall um 
eine sowohl in 2* als auch in F° feste Achse rotiere und wählen diese Achse z. B. zur 
C-Achse von 2* und zur z-Achse von 2°. Bei Schwingungen um diesen stationären 
Bewegungszustand kommen zur Eigendrehung noch zusätzliche 


Drehungen um die zur Drehachse der Eigendrehung senkrechten 


i Paz ; u : : ı|2 
Achsen. Die Winkel dieser Zusatzdrehungen wollen wir als klein er 
betrachten und werden folglich ihren Cosinus gleich eins und || 
ihren Sinus gleich dem Winkel selbst setzen. Dann fällt die AB- Br 


Ebene von 2* näherungsweise mit der xy-Ebene von 2° zu- 
sammen (Bild 9). Die C-Achse weicht von der z-Achse um einen 
kleinen Winkel ab, dessen Projektionen auf die yz- bzw. x2-Ebene « 
bzw. sind. 

Wie man aus Bild 9 oder aus dem Transformationsschema (34) h 
für kleine Winkel a und ableiten kann, gelten zwischen den Im- r 
pulskomponenten die Beziehungen: Bild 9. 











ee * 





Pe . 8 , P nu Z. angew. Math. Mech. 
354 Magnus, Über die allgemeinen Bewegungsgleichungen starrer Körper Bd. 4 Nr. 4 Dez. 1945 
J.—=J jcosp — Jzsng+JcePp | 





AR Jı sin g + Jp cosq Jc a (44). 


J,=J ,lasınd - P COS) + Jplacosg + sing + J | 


Darin hat man sich für J,, Jp. Je die Werte (42) eingesetzt zu denken. Meist läßt sich das 
körperfeste System so legen, daß 2°* Hauptachsensystem wird. Dann bleibt von (42) nur 
noch übrig: 

e Jı= Iıwu,, 


j Ip un (45) 


JI=%r u, | 


Das läßt sich stets erreichen, wenn die Drehachse der Eigendrehung Hauptträgheitsachse ist. 
Wenn es jedoch darauf ankommt, bei der Rechnung auch den Einfluß von Unwuchten zu er- 
fassen, die durch eine Abweichung der Drehachse von der Hauptträgheitsachse hervorgerufen 
werden, so läßt sich die Verwendung der allgemeinen Impulsausdrücke (42) nicht umgehen. 


Der Kunstgriff der weiteren Rechnung besteht nun darin, daß die im allgemeinen be- 
kannte und meist sogar konstante Eigendrehung des Körpers von den Komponenten der Zu- 
satzdrehungen abgesondert und für sich behandelt wird, d. h. die Komponenten der Zusatz- 
drehungen werden so angesetzt wie bei einem nicht um die C- Achse drehenden Körper: 


ws: u.cosy + u,sing + u; (asing —Pcosgp), 
"Bp=— u, sing +u,cosyg tu. (acosp+Psing), 
0= wP-watn,. 


Man sieht daraus, daß das letzte Glied der ersten und zweiten Gleichung vernachlässigt werden 
kann, da es von zweiter Ordnung klein ist. Das folgt sofort durch Einsetzen des Wertes 


für «, aus der dritten Gleichung. Wir können also weiterhin: 
’ 4 f ’ > 
u, u,cosp + uy,sing, 
u" = — using + u), cos q 


setzen und erhalten dann für die Impulskomponenten (44) mit den Werten (45), also für den 
Fall, da& 2* Hauptachsensystem ist: 


ar : er 
„+ > (9 N O5) sın 29 Uy ._ , Ue pP 


mr —_ Q, cos? q E= e77 sin? q ) uU, 


J,= 5 (94 — %)sin2g u, + (9 ,sin’p + 9% cos’y)u,— Ic tea 
(46). 





a . r . 9 | ’ 
J.= 15 ( O, 9,) sın 29 - P (9, cos” q .- Ip sın“ q )| u, 





> 1 
+ 5 (d— 9,)sin29+a(9,sin’g + Ogeos’p)u, + Icu; 


Damit hat man die in die allgemeinen Gl. (9) bzw. (20) einzusetzenden Komponenten 
des Relativimpulses nach 2° gefunden. Diese Werte stimmen um so besser, je größer die 
Eigendrehung «,. gegenüber den anderen Drehgeschwindigkeitskomponenten ist, je mehr also 
der Gesamtimpuls des Körpers in die C-Achse fällt. Das ist sicher dann der Fall, wenn 
u,=uya-— u,ß gegenüber «,. vernachlässigt werden kann, was bei hinreichend schnellen 
Kreiseln immer möglich ist. Da u‘ eliminiert ist, kann man jetzt aus (46) nicht wie im Falle 
der kleinen Pendelschwingungen den Trägheitstensor nach 2° bestimmen. Man kann also 
auch nicht die expliziten Gl. (18) oder (19) verwenden. 


Bei der Berechnung des in (9) bzw. (20) vorkommenden Führungsimpulses \j; kann man 
sich meist damit begnügen, die Glieder ohne die kleinen Faktoren a bzw. 5 mitzunehmen, 
wenn man nicht sogar den Führungsimpuls ‘7; vollkommen gegen den Relativimpuls \% ver- 
nachlässigen kann, was bei nicht zu schnell drehendem Bezugssystem sicher erlaubt ist. 
hat bei Vernachlässigung der von erster Ordnung kleinen Glieder die Komponenten: 
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; a 1 E 
Jx=(9 ,cos®’p+ 9, Sin’ y)u,. +5 (9, — 9,)sin2pu, | 


R Br ) “ 17). 
Jy=5(9, - )Ssin2g u, +(9 sing + 9,cosgy)u, 1x) 


Jr - w?7 U> 


Mit den Werten (46) und (47) bekommt man nun aus (20) die sehr allgemeinen Nähe- 
rungsgleichungen für die Bewegungen eines Kreisels auf bewegtem Bezugssystem. Sie bieten 
im allgemeinen sehr viel weniger mathematische Schwierigkeiten als die strengen Gleichungen 
und sind in Sonderfällen (z. B. symmetrischer Kreisel auf konstant drehendem System) sogar 
elementar lösbar. In der allgemeinen Form (20) umfassen die Gleichungen auch noch den 
Fall, daß der Körper keinen festen Punkt im Bezugssystem besitzt, also selbst noch irgend- 
welche Bewegungen gegen I°° ausführt. Ist dagegen ein Punkt des Körpers in Y'° fest, so 
reduziert sich (20) auf: 

© ° 
sugar tngeR . :. 252.4 240 


Aus dieser Form lassen sich gleich einige bekannte Sonderfälle herausnehmen. Z. B. bekommt 
man mit \\y=0 und symmetrischem Kreisel die von Schuler [5] für langsam drehendes 
Bezugssystem aufgestellten Kreiselgleichungen. Nimmt man außerdem u=0 (Kreisel im 
Inertialsystem), so kommt man zu den seit langem bekannten linearisierten Kreiselgleichungen, 
die häufig die Föpplschen Gleichungen genannt werden, aber schon vor Föpp!l z. B. bei 
Routh'’) zu finden sind. 


Weiterhin bekommt man aus (48) ebenfalls für u—=0 (Kreisel im Inertialsystem) aber 
unsymmetrischen Körper die Formeln von Fuchs [5] und Weigand [9]. Diese letzteren 
Formeln sollen kurz hergeleitet werden, um die Form der bei unsymmetrischen Körpern auf- 
tretenden Glieder an einem Beispiel zu zeigen. Mit u=0 reduziert sich (48) auf den 
Impulssatz: 

je} . 
a 
Je 
Setzt man darin die Werte (46) ein, so erhält man aus den beiden ersten Komponenten- 
gleichungen nach einigen trigonometrischen Umformungen: 


1 u. u 
5 l9 4; +9) +19 1 — Ip)eos2yplu, +5 (9 4 - %)Ssin2pü,+ eur u), 


| | . ’ ’ 
u uc| > (9 - 9 )sin2pu, +19 ,— 9) c0s2g u) s= MM. 
1 | (49). 
39, 9,)sin2pu,+ 59149) —(9 ,— 9p)ec0os2Yp] u, — 9 ur u) 


Y 





| 2 ‚ 
+ue[i®, O,)eos2pu,.+519,- 9)sin2yu,|=M, 


Dabei ist jedoch der Fall angenommen, daß u —=g konstant ist, also «.—=0. Dieser 
Fall kommt praktisch sehr häufig vor, weil der Kreisel meist durch einen Antrieb auf kon- 
stanter Drehzahl gehalten wird. Die Koeffizienten haben beim unsymmetrischen Kreisel die 
Periode 2 u, weil das Trägheitsmoment während eines Umlaufes zweimal seinen Höchst- und 
Tiefstwert annimmt. Die Gl. (49) entsprechen genau den Bewegungsgleichungen, die Fuchs [8] 
und Weigand [9] für den unsymmetrischen Kreisel verwendet haben. 


Im Falle des symmetrischen Kreisels (9 ,= ©, = ©) geht (49) in die bekannten Föppl- 
Gleichungen über: 

u, +%ucuy=M; | E 

22 er (50). 

O9u,--Ocuru,=M, 

Seine ganze Anpassungsfähigkeit erweist das angegebene Näherungsverfahren aber erst, 

wenn man Veıbände starrer Körper zu untersuchen hat. Man braucht dann nur die Einzel- 

impulse der Körper zu summieren und hat damit in (9) bzw. (27) einzugehen. Wir wollen 

das an dem im Abschnitt 7 streng berechneten Beispiel des kardanisch gehängten Kreisels 

zeigen, beschränken uns dabei jedoch der Einfachheit halber auf den Fall eines symmetrischen 
Kreisels. Für die Summen der Impulskomponenten hat man: 


12) Siehe [1], Kap. V. 


23° 





EiE? Dar or er er ereune 
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N”K=J, tIo + Ju =Our + 9cucp+9ur + sur. 
Ju, +J.+J, = 9 u, — % wca+ Ou,u, +0, 


”I.=J,+0+0. 


Die ersten Glieder folgen aus (46) mit 9,= 9%,. Setzt man das Trägheitsmoment des 
gesamten Systems um die .«- bzw. y- Achse 


«=9+ + Os. 


Iy + u, 
an, so hat man die Impulskomponenten: 


ee. u + IucPß. 


Ju, = 9,uy — cur a. 
woraus nach dem Impulssatz sofort die gegenüber (50) allgemeineren Föpplschen Gleichungen 


O,u,.+ Ocuruy= M,;, 


O,uy — Icuru, = My 
folgen, bei denen die zum Kreisel hinzugekoppelten Massen mit berücksichtigt sind. 


Die Einflüsse der mit dem Kreisel verbundenen und zwangsläufig mit ihm gekoppelten 
Zusatzmassen lassen sich auf diese Weise auch für den unsymmetrischen Kreisel stets dann 
erfassen, wenn die Zahl der Freiheitsgrade des Verbandes durch die hinzugekoppelten Massen 
nicht erhöht wird. Bringen jedoch die Zusatzmassen neue Freiheitsgrade, so hat man dem- 
entsprechend mehr Gleichungen mitzunehmen, in denen dann auch die von den einzelnen 
Teilkörpern aufeinander ausgeübten Reaktionsmomente zu berücksichtigen sind. Man muß 
dann wieder auf die allgemeine Form (27) zurückgreifen, bzw. auf deren Vektorzerlegung, die 
sich durch Summierung und Hinzufügung der Reaktionsmomente aus (20) ergibt. Diese 
Gleichungen sind insbesondere zur Berechnung von gekoppelten Kreisel- und Pendelsystemen 
geeignet, wie sie in der Technik häufig anzutreffen sind, z. B. beim Zwei- und Dreikreisel- 
kompaß, bei dem -mit Pendeln gekoppelten Ein- und Zweikreiselhorizonten und vielen anderen 
Geräten. Bei der Berechnung des Verhaltens derartiger Geräte auf bewegtem Fahr- oder 
Flugzeug sind die Näherungsgleichungen unentbehrlich geworden, weil sie selbst in solchen 
Fällen noch gewisse Aussagen zu geben imstande sind, in denen eine strenge Lösung an 
mathematischen Schwierigkeiten scheitert. 
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Fehlerabschätzung für das Iterationsverfahren zur Auflösung 
linearer Gleichungssysteme. 
Von 1. Collatz in. Karlsruhe. 
Vorgetragen auf der GAMM-Tagung 1942 in München. 


Für das Iterationsrerfahren zur numerischen Auflösung linearer Gleichungs 
systeme wird eine Fehlerabschätzung durchgeführt unter der Voraussetzung, 
daß die Hauptdiagonalelemente genügend stark überwiegen. Die Gaufßsche 
Transformation ist nicht immer geeignet, um aus einem beliebigen Gleichungs- 
system ein solches mit überwiegenden Hauptdiagonalgliedern zu erzeugen. 


1. Fragestellung. 


Vorgelegt sei ein lineares Gleiehungssystem 


ı 


ee 20 20) 5-5 20 SEE © 


N 


Ir 


N 


oder in Matrizenschreibweise 
12% u EEE a a ;  ' 
die Determinante a;;| sei ==0. 
Bei dem „Iterationsverfahren in Gesamtschritten“ ') verbessert man Näherungswerte für 
die .«; schrittweise und berechnet aus der v-ten Näherung x;"” eine (+ 1)-te Näherung x; +” 
nach der Vorschrift 


n 
i l y Ba 
Era r; \ d;y u. Ar ee were 
d;; a 
k 1 
k ! 


! l n 
b | v “ Y x 
a,” zii (v; \ d;y 2, De \ d;y Aal ta a a 
d;; — — | 
k=] k=i-+l 


Diese beiden Iterationsverfahren konvergieren nicht immer, und wenn sie konvergieren, nur 
dann für die numerische Rechnung rasch genug, wenn die Hauptdiagonalelemente a;; die 
übrigen Elemente «a;,. genügend stark überwiegen. 
Im folgenden wird auf die beiden Fragen eingegangen: 
I. Wie stellt man zu einem beliebigen vorgegebenen Gleichungssystem (1), bei dem die 
a;; nicht genügend überwiegen, ein neues Gleichungssystem mit genügend über- 
wiegenden Hauptdiagonalelementen her? 
2, Wie kann man aus den nach den Iterationsverfahren berechneten Näherungswerten 
ihren Fehler abschätzen ” 
2. Hilfsbetrachtung: 2 Unbekannte. 
Bei 2 Gleichungen mit 2 Unbekannten 
da, rAda,—t,, | 


4,X%,+4,32,=r, 


(5) 


lassen sich die Verhältnisse vollständig überblicken. Für die Fehler &;”” der »-ten Näherung 
g;’ r 2,” — %; G= 3 


erhält man beim Iterationsverfahren in Gesamtschritten 


e;" D gg, PD 
wobei zur Abkürzung 
A,. dl;, 
a 
@,, don 


I, R. v. Mises-H. Pollaezek-Geiringer: Z. angew. Math. Mech. Bd. 9 (1929), S. 58 bis 77. 
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gesetzt ist, während beim Verfahren in Einzelschritten 


v+» 
; = 


) 


eo 
€ aAE; 


wird: Man sieht also: Beide Iterationsverfahren konvergieren genau dann, wenn la!<1 ist; 
man kann a als Gütezahl für die Konvergenz auffassen, der Fehler gelit auf den Bruchteil «a 
herunter, und zwar beim Gesamtschrittverfahren in 2 Schritten, beim Einzelschrittverfahren 
in 1 Schritt. Das Gesamtschrittverfahren benötigt zur gleichen Genauigkeit im vorliegenden 
Fall genau doppelt so viele Schritte als das Einzelschrittverfahren. 
Bei drei Gleichungen mit drei Unbekannten liegen die Verhältnisse nicht melr so ein- 
2 3 


fach. Bei einem Gleiechungssystem (2) mit der Matrix A \ 2 2 >| konvergiert (bei be- 
| 1. 2/ 
liebigem Vektor r) das Einzelschrittverfahren, aber nicht das Gesamtschrittverfahren und bei 
4 u: 9 5) 
einem Gleichungssystem (2) mit der Matrix W =|1 1 | konvergiert das Gesamtschritt- 
.) ») 1 


verfahren, aber nicht das Einzelschrittverfahren. 
3. Die Gaußsche Transformation. 
Bei der Gaußschen Transformation multipliziert man Gl. (2) von links mit der aus A 
durch Spiegelung an der Hauptdiagonalen hervorgehenden Matrix W und erhält als neues 
Gleichungssystem 











DES mE BERN: BETEN 


Die Matrix ® ist symmetrisch. Die zu ihr gehörende quadratische Form ist positiv definit 
und daher!) konvergiert das Einzelschrittverfahren bei dem neuen Gleichungssystem (6). Die 
Theorie hat dadurch eine große Abrundung erfahren, für die praktische Rechnung hat jedoch 
die Gaußsche Transformation einige unangenehme Eigenschaften. Das neue Gleichungs- 
system (6) gehört zwar zu einer positiv definiten quadratischen Form, aber die Hauptdiagonal- 
elemente überwiegen oft nur in so schwachem Maße, daß das Iterationsverfahren sehr langsam 
konvergiert und für praktische Zwecke völlig unbrauchbar sein kann. Ferner kann die 
Gaußsche Transformation die Konvergenzgüte verschlechtern. 

Zur Erläuterung wird ein Gleichungssystem mit 2 Unbekannten betrachtet, Gl. (5). 
Hier hat das Glejehungssystem nach Nr. 2 als Gütezahl 


Ad,ad;, 
a= : 
#0 


während die Gütezahl des nach Gl. (6) berechneten transformierten Gleichungssystemes 
(a,a,+ a 


(a,,’ nu A,,') (G,;° + dy5') 


21 A,5)" 


B= 


lautet. Wie es sein muß, ist | <1,. Fülrt man die Quotienten 


(da 


HK = : . y= 


. d,, £ d, 


ein, so kann man in einer .«-„y-Ebene die Bereiche 
aufsuchen, für deren Punkte die Gaußsche Trans- 
formation die Konvergenz des Iterationsverfahrens her- 
beiführt (\a! >1), verbessert (1>!a!>|ß|) oder ver- 
schlechtert (1 > >|a|). 
Bei positivem ., y sieht man, daß nur die beiden 
Fälle |a|>1 und 1> | > al auftreten. Beschränkt 
7. man sich auf positive Werte der Quotienten u, —. 
T=z, 11 22 
"so sind also nur die folgenden 2 Fälle möglich (ab- 
gesehen vom Grenzfall a=1): Entweder ist beim Aus- 
gangssystem a >1, dann führt die Gaußsche Trans- 
formation die Konvergenz herbei, oder es ist O<a<1l, 
dlann verschlechtert die Transformation die Konvergenz. 
. Wendet man also die Gaußsche Transformation mehr- 
de. ie alte ee Weiler mals hintereinander an, so kann sie möglicherweise 
wird die Konvergenz des Ilterationsverfahrens das erste Mal Konvergenz herbeiführen, hei jedem 
weiteren Male jedoch verschlechtert sie die Konvergenz 
des Iterationsverfahrens. 








rn 


U] herbeigeführt [ITD verbessert EZ verschlechtert 
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Als einfaches Beispiel betrachten wir das Gleichungssystem 


0, —22,=0 
) - >} 
x X, . 
mit der Lösung »,=2, .„,=1 und der Gütezahl a=4. Bei der Gaußschen Transformation 
entsteht das neue Gleichungssystem 
> HK - { ss == 6. 
{ x, 7 > Jg > 3 


z 16 5 Pr 
mit der Gütezahl 5==0,64 und bei nochmaliger Gaußscher Transformation das Gleichungs- 
. +) 


system 
Hr 2, =2, 


400, +412,= — 39 


i 16 ® : Bar i i 
mit der Gütezahl (681 > 0,95181. Bei dem letzten Gleichungssystem geht beim Iterations- 
x 


verfahren der Fehler bei jedem Schritt also nur um 5°/o herunter, d. h. das Verfahren ist 
praktisch unbrauchbar. 


4. Das Umformungsverfahren von C. Runge. 


Nach dem Vorangehenden erscheint die Gaußsche Transformation trotz der theoretischen 
Abgerundetheit für die numerische Rechnung noch nicht als das ideale Mittel zur Erzeugung 
überwiegender Hauptdiagonalglieder. Es sei hier vielmehr auf ein bereits von C. Runge 
angegebenes Verfahren?) hingewiesen; man stellt durch Linearkombination der gegebenen 
Gleichungen neue her, bei denen außerhalb der Hauptdiagonale nur noch kleine Zahlen 
stehen. Daß das stets möglich ist, beweist das Gaußsche Eliminationsverfahren, bei dem 
man ja durch Linearkombination erreichen kann, daß außerhalb der Hauptdiagonale überall 
Nullen stehen (damit ist dann das Gleichungssystem aufgelöst). Beim Rungeschen Ver- 
fahren verzichtet man auf die Forderung, außerhalb der Hauptdiagonale streng Nullen zu er- 
zeugen, sondern begnügt sich damit, daß außerhalb kleine Zahlen stehen, so daß das Iterations- 
verfahren gut konvergiert. Gegenüber dem Gaußschen Eliminationsverfahren hat das Runge- 
sche Verfahren die Vorteile: 

l. Die Erzeugung kleiner Zahlen außerhalb der Hauptdiagonale ist sehr viel leichter als 

die Erzeugung von Nullen, oft kann man einem Gleichungssystem unmittelbar an- 
sehen, welche Zeilenkombinationen schnell zum Ziel führen. 
Dadurch, daß man beim Rungeschen Verfahren die Zeilenkombinationen so vor- 
nehmen kann, daß man nur mit ganzen oder zumindest genauen Zahlen multipliziert, 
addiert und subtrahiert, hat man keinen Genauigkeitsverlust und kann z. B. ein 
Gleichungssystem mit dem Rechenschieber auflösen, während man beim Gauß- 
schen Verfahren zur Erzeugung von Nullen mit im allgemeinen ungenauen Zahlen 
multiplizieren muß und bei größerer Anzahl von Unbekannten nicht mehr mit dem 
Rechenschieber reelınen kann. 

Das Rungesche Verfahren wird so lange fortgesetzt, bis die Hauptdiagonalelemente 
„zenügend* überwiegen, d. h. bis das „Zeilensummenkriterium“ 


n 


diil > b; mit b; = Uik a E 
e—=ı 
ki 


für alle © erfüllt ist. Es werde noch 


ı 
u A;;| b; | 8 
ER N a ne 
er, 


gesetzt. Dann ist o ein Maß für das Überwiegen der Hauptdiagonalelemente. Je kleiner o, 
desto besser konvergiert das Iterationsverfahren. Als grober Anhaltspunkt für numerische 
Rechnungen sei die Forderung o < 2 genannt. Wie beim Eliminationsverfahren hat man auch 
beim Rungeschen Verfahren die bequeme Kontrolle für die Summen aller Koeffizienten 
und rechten Seiten einer Zeile. 





®,C, Runge H. König: Numerisches Rechnen, Berlin 1924, S. 187. 
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5. Fehlerabschätzung’). 
Es gilt der 
Satz: Ist bei einem linearen Gleichungssystem 


” 
Ad; Ar f; (4 l 3 ı) 
h l 
die Bedingung 
h; ee d;y x Idi; 
KB 
k ! 


für alle ö erfüllt, so gilt für die nach der Vorschrift (4) des Iterations- 
verfahrens in Einzelschritten®!) berechneten Näherungen .r;”” die Fehler- 
abschätzung 


at ER a ee EL ARE eye 
.: b; er) > . x 
Dabei ist o—= Max - n und 5° der Betrag der maximalen Anderung 
ä i A;; ; 
beim letzten Iterationsschritt ö”’ = Max .r; "+", 
Beweis: Es sei Z” der größte Betrag der n Fehler 
= 2;® rn; 
also 
3” = Maxle," 
i 
Für den Fehler gilt nach (4) 
i l n 
l e ” N 
a + I \ d;, BE - 1. \ dj, ers F { ? ‚ £ Ä E (10). 
UN —— 
k 1 7 i l 
on . iu 0 . ® = 
Führt man nun die Abkürzung «=, ein, so ist nach Voraussetzung 
9; l+o 
b;Sula;; 
mit «<1 und es gilt weiterhin 
er 
denn für i=1 liest man |e,”+P)< uZ®’ unmittelbar aus (10) ab, da dann die erste Summe 
. . 1 ’ . ’ N ’ 
leer ist, und für die weiteren i kann man rechts bei der ersten Summe |. "+P’|ZuZ” <zZ” 


einsetzen und erkennt so (11) für alle ö als gültig. Nun sollte Z”*+»” das Maximum der n 
Zahlen |2;”*"| sein, muß also unter diesen Zahlen vorkommen, d. h. es gilt 


bi + 1) u Ag 
und damit 
ar u F Au 


Das ist die bekannte Konvergenzaussage Z” >0 für »>x. Nun sei 0°” der größte Betraz 
der beim Übergang von der v-ten zur (v—+1)-ten Näherung aufgetretenen Än lerungen 


Pa Max x,” 4-1) x ıy) — Max 6” #) m 


! D 
? ? 


und es sei x, diejenige Unbekannte, für die |e;"”, den Maximalwert Z” annimmt: \*, "= Z", 
Für dieses o gilt: 


y y 4 ‚ 1 r r 3 
PL > & Be FR 1) » ge.” er ı) Al 2# ı) 


1 | I ze+n 
0 
Damit ist 
2” + 1) o PL 
bewiesen. g 


3) Bei R. v. Mises-H. Pollaezek-Geiringer: Z. angew. Math. Mech. Bd. 9 (1929), S. 65, findet sieh eine 
Fehlerabschätzung, allerdings nur für das Iterationsverfahren in Gesamtschritten und zudem bei Erfülltsein des Spalten 
summenkriteriums. Dieses hat gegenüber dem hier zugrundegelegten Zeilensummenkriterium die Nachteile, daß man 
zur numerischen Anwendung erst die Zeilen dureh die Hauptdiagonalelemente durehdividieren muß und daß man eine 
Aussage über die Summe der Fehlerbeträge und nicht direkt über den Maximalfehler erhält (höchstens dureh die grobe 
Abschätzung: Maximalfehler <Summe der Fehlerbeträge). 

4) Wie man leicht beweist, bleibt der Satz richtig, wenn man „Vorschrift (4)“* und „Einzelschritte“ dureh ‚‚Vor- 
schrift (3)* und „Gesamtsehritte‘ ersetzt, 
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Zur Untersuchung der Güte der Fehlerabschätzung werde folgendes einfache Beispiel 
betrachtet, bei dem sich zeigt, daß die Fehlerschranke (9) und der tatsächliche Maximalfehler 


sich beliebig nahe kommen können. 


(l+o)x2,+o%2, 


BIS TOR. 2, =, 0, 
u |. 
Die exakte Lösung lautet 
0 
I, = | +0 
der Maximalfehler beträgt mithin N = 
u 6 


Geht man bei dem Gleichungssystem 


so erhält man bei einem Iterationsehritt 


Der Parameter o stimmt mit dem in (8) eingeführten o überein, man erhält daher mit 


ö=1 die Fehlerschranke o-ö 
Maximalfehler unterscheidet. 


o, die sich für o>0 prozentual um beliebig wenig von dem 


117 


KLEINE MITTEILUNGEN 


Zur Konvergenz des Newtonschen Ver- 
fahrens für Gleichungssysteme. (\us dem 
Institut für aneewandte Mathematik und Darstel 
lende Geometrie der Technischen Hochschule 
Braunschweig.) In einer Dissertation ®). die wegen 
des Krieges nicht zedruckt werden konnte, hat 
Herr Karl Bussmann eine Reihe von Sätzen 
und Fehlerabschätzungen gezeigt. durch die einige 
Sätze von A. Ostrowski?) und Fr. A. Wil 
lers°®) wesentlich verallgemeinert werden. Von 
ihnen seien die wichtigsten in etwas anderer For 
mulierune ohne Beweis hier mitgeteilt. 

Zu lösen sei das Gleichungssystem 

U E TORR- RAR u 9% An v8 


Ln)=UV. %= 


Die y, seien Funktionen des Vektors x = (2,,'*'.2,,) 
und sollen stetige partielle Ableitungen bis zur 
2. Ordnung besitzen. Es sei g der Vektor mit den 
Koordinaten 9.9%s....91. Zur Abkürzung werde 
gesetzt: 

zZ 0° gu 


Je» Ir) 


Matrix {9.3=®, 
Or, dxu . In 


O7 


Determinante von ® = @. 
SG, CWW)=E@. 


Adjunkten von ©: Gx;, 
IK)EH, Su 
Das Newtonsche Verfähren besteht in der 
folgenden Iteration: Man geht von einer Näherung x, 
v=0(,1,?,...) der gesuchten Stelle x aus, rechnet 
sich den Vektor q, mit den Koordinaten 9, (x,) aus 
und gewinnt durch lineare Transformation dieses 
Vektors einen neuen Vektor x,+,. nämlich 


b4, =b +5 9%. 


Dabei ist &, die an der Stelle x, genommene Matrix 
6) - - 7", also 


4 =t |. 


1) Eingereicht bei der Technischen Hochschule Braun 
schweig 1940. 








> 


2, A. Ostrowski: Konvergenzdiskussion und Fehler 
abschätzung für die Newtonsche Methode bei Gleichungs- 
systemen. Commentarii Mathematici, Helvetici, Bd. 9 
(1937/38), S. 79 bis 103. 

3) Fr. A. Willers: Zur Konvergenz des Newtonschen 
Näherungsverfahrens. Z. angew. Math. Mech., Bd. 18 
(1938), S. 197 bis 200, 


sezeichnungen: 1. a sei eine Lösungsstelle 
des Gleichungssystems g(r) =. 

2, Die Koordinatenbeträge des Vektors ı ad, 
also die Zahlen #2 — 4x heißen die n Ab 
weichungen der Näherung y.. 

3. Die größte dieser n Abweichungen nennen 
wir den Maximalfehler ö, der Näherung y.. 

4. Die Koordinatenbeträge des Vektors u +, Xı. 
also die Zahlen |, + ,. Tyz.. heißen die n Koı 
rekturen der Näherung d.. 

5. Die größte dieser n Korrekturen nennen wiı 
dieMaximalkorrektur d» der Näherungx 

6. Einen abgeschlossenen achsenparallelen Würfel 
mit dem Mittelpunkt x, »=0,1,.2,...| und deı 
Spannweite (d.h. Kantenlänge) s bezeichnen wir 
durch das Symbol W [x.. s]. 

7. In einem vorgegebenen abgeschlossenen Ge 
biet bezeichnen wir eine obere Schranke für eine 
Funktion |g(x) | mit g', eine untere mit ‚g; eine obere 
Schranke für die Werte der rn Funktionen g, 

., 9n bezeichnen wir mit gx usw. 

Satz 1 (Ein Konvergenzkriterium). 

Es sei schon bekannt. daß im Würfel W |xo.s 
wenigstens eine Lösung aliege. Im Würfel W[r,. 2 
schätzen wir G und sämtliche 9... gmau. 
ab und berechnen aus den auftretenden Schranken 
eine dem Würfel WIx,. s| zuzuordnende Kenn 
zahl o, 

Zeigt sich nun, daß die Spannweite 

s< ae, 
so gilt zweierlei: 
I. a ist die einzige Lösung in W ix. |. 
II. Die zu Io eehörende New to nsche Folee kon 
vergiert gegen a. 
Zur Ermittlung von a werden in W|1.2s] be 


stimmt ,G. gu: 2, 0x: 3 Dann wird 
4-,G 


2 n—ı 7 
nen! ga" Guam 














Satz 2 
mation). 


Konvergenz bei genügender Approxi 


Wir betrachten eine Würfelfolge W, = W Ix.s,) 
v—=0,12,...,für del. > >5,>...,2.5>0, 
3. W,+, ganz in W,. 4. a in jedem W,. Dann 
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eibt es in dieser Würfelfolge einen ersten Würfel, 
dessen Spannweite s,< als die zugeordnete Kenn- 
zahl @, ist, wobei @, nach Satz 1 zu berechnen ist. 
Die zu seinem Mittelpunkt x, gehörende Newton 
sche Folge konvergiert also nach Satz 1 gegen a. 
D. h.: Das Newtonsche Verfahren konvergiert 
stets. wenn die Ausgangsnäherung x, die Lösung a 
schon „genügend” approximiert. 

Satz3 (Fehlerabschätzung für den Fall, daß die 
Existenz einer Lösung bekannt ist). 

In W [xo. s] liege eine’ Lösung a und es sei s« 
als die Kennzahl «a |Satz 1). Man gehe in der 
Newtonschen Folge x, so weit. bis 


s\2 n 


| ; 
a In — | 
Dann ist der 
! R Bo B88n—1d,) 
Maximalfehler d, +, d?'1+ s 
a N 1 a 


Satz 4 (lteration mit konstanter Matrix). 

Es sei schon bekannt, daß in W |x.. s) eineLösung a 
liege. Aus Abschätzungen von G(Xo) » Gx2(ko 
und ferner von G,„zu(x) in W 1x. 2s| errechnen 
wir eine Kennzahl 3. Zeigt sich, daß die Spann- 
weite 


sam 


so gilt zweierlei: 
l. a ist die einzige Lösung in W |x.. S'. 
Il. Die Vektorfolge 
bu, =ub— &°'%; 
bei der also die Matrix ®5! bei jedem Schritt die 
gleiche bleibt. konvergiert gegen a. 

Die Berechnung von 3 erfolgt so: Es sei G (1,) 
=/G,|, |Q@xi (0) v, und in W[x,2s] sei 
Iniulf)| S Alain. 

Dann ist 
G,| 


In? dv. Gasiı 





Satz 5 (Erster Existenzsatz). 

Es sei aus x, bereits die Näherung x, berechnet. 
In W[x.4do] rechnen wir uns eine unten ange- 
eebene Kennzahl y aus: ferner berechnen wir 
G (x) = G,„. Zeigt sich, daß die Maximalkorrek 
tur d, von x, (also die Viertel-Spannweite jenes 
Würfels) 

du,<Y % 
so gilt dreierlei: 
I. In W [x..4 d,| liegt genau eine Lösung. 
II. Die zu x, gehörende Newtonsche Folge kon 
vergiert gegen a. 


Ill. Für den Maximalfehler gilt 





ER - 5 Be 
speziell also 
Be 
Die Berechnung der, Kennzahl y erfolgt so: In 





Wfx.,4d,) werden bestimmt 9,.; und .; 
Dann ist 
1 


nen! a" "Qosau 





Satz 6 (Zweiter Existenzsatz). 


Unter denselben Voraussetzungen wie im Satz 5 
rechnen wir uns im Würfel W [x,. 4 d,| eine Kenn- 


zahl o aus. Zeigt sich, daß die Maximalkorrek- 
tur d, von Xu 


so eilt dreierlei: 
I. Es gibt in W[x,.4d,| wenigstens eine 
Lösung. 
Il. Die zu x, gehörende Newtonsche Folge kon- 
vergiert gegen eine solche. 
Ill. Für den Maximalfehler gilt 


20? 


a 
: A 7 d? 





Dabei wird o in W |xr,.4 d,| berechnet dureh 
G 


n® Gy; Qua. 


Satz 7 «(Dritter Existenzsatz; Verschmelzung von 
Satz 2 und Satz 6). 

Unter denselben Voraussetzungen wie im Satz 6 
rechnen wir uns im Würfel W |x,,8 do] eine Kenn- 
zahl r aus. Zeigt sich, daß die Maximalkorrek- 
tur d, von X, 


so gilt zweierlei: 
l. In W [x 4d,) liegt genau eine Lösung a. 
Il. Die zu x, gehörende Newtonsche Folge kon- 
vergiert gegen a. 
t wird bestimmt durch Absehätzungen in W [yo ,8 dy]. 
nämlich 
(r 
T=—_ —=— . 
n® Gy 3 Juin 
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Praktische Interpolation. | chrstelligen 
lafeln, wie sie zur Durchführung von Rechnungen 
ınit größerer Genauigkeit verwendet werden, ist 
die lineare Interpolation nicht mehr ausreichend. 
Die Wiederholung von Interpolationen höheren 
‘irades für mehrere Werte eines Argumentes ist 
umständlich und zeitraubend. Im folgenden soll 
ein Verfahren zur direkten Berechnung der häufigst 
vorkommenden elementaren Funktionen mit 
Hilfe eines größeren Tafelwerkes angegeben wer 
len, las in der Rechenpraxis große Vorteile bietet. 

Es sei since für ein beliebiges mehrstelliges € zu 
bestimmen. Wir zerteilen das gegebene x in 

E Fr 


ERFYT HT 


n 
oder 
zZ I, no € 
€ & 96 s 
== 


und setzen voraus, daß die Sinus- und Cosinus 
funktionen für die Argumente z,, &, &.... tabu- 
liert vorliegen. e, wird dann eine so kleine 
Zahl sein, daß 


sın #, = E&,. COS E, >= ® 


Der Fehler 4 sine—e für kleine Werte des 


S 


Argumentes s ergibt sich zu 


x 
der Fehler 
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erößer; so er- 


{=} 


ist für dieselben Werte & wesentlich 


hält man für e= 00-001: 


Wir entwickeln nun fortlaufend 


sinr=sin (tr, ZEo)=sinz, cos&E + cos 7, sine. 


COS 2 = 008 (7, + &) = 008 .L, CO8E sin .r, sin €: 


sine= sin &, 7 ö sin &, COS ö + c08 £, sin Ö, 


cos E=cos(ie, + COS &, 6080 -— Sin &, Sind; 


und setzen schließlich 
sin &; & 
voSsE$; ; 


Die erforderlichen Multiplikationen und Additio 
nen sind auf einer größeren Rechenmaschine ohne 
\ufschreibung von Zwischenergebnissen durchzu 
führen. Analoge Entwicklungsformeln, nur mit 
durchwegs positiven Vorzeichen, gelten für die 
Hyperbelfunktionen Sin. und Coj«. 


Die Zweckmäßigkeit und Schnelligkeit des Ver- 
fährens soll das folgende Beispiel erläutern. Die 
Zahlenwerte wurden dem bekannten Werk von K. 
Hayashi: „Sieben- und mehrstellige Tafeln der 
Kreis- und Hyperbelfunktionen“ (Berlin 1926) ent- 
nommen. Es seien die Funktionen sin U. cos .«. 
Sin und Coj für das Argument «© — 3.5472398 
auf 8 Ziffern genau zu berechnen. Für die gefor 
lerte Genauigkeit ist es ausreichend, cos 0,0398 
hzw. Koi 0.032398 = 1 zu setzen. Es empfiehlt sich 
folgende Aufschreibung: 








Arg. sin cos Zin Goj 
x =%54 0,33755092 — 0,92168003 17,218953 17,247966 
#11, 0272 0,027199038  0,987408 0,027200062 1,082592 
d — N,04398 0,03398 1 0,4398 1 


0,027239737 0,087379 0,027.39863 


1,042621 


— 0,39461348—0,91884721| 17,344277 | 17,373081 





Die Werte der Funktionen in den ersten beiden 
/Jeilen werden dem Tafelwerk.entnommen, dann be- 
ıechnet man die Funktionen für das Argument 
‘4. Zeile) und dann schließlich die Werte der letz 
ten Zeile. 


Zum Vergleich sei sin © für denselben Wert von 
" auf die übliche Weise unter Benutzung der 
Newtonschen Interpolationsformel für äquidistante 
\rgumentwerte berechnet. Diese lautet mit der 


. . LT To 
Abkürzung = 
Ir 
i Wr 
Y YrAdiyıTr > Pyt 
Z(A—-VD(i—2 d. n—1 
k & j s Yo 
n! 


Mit 1 erhält man die Formel für lineare Inter 
polation, n = 2 ergibt quadratische Interpolation 
usw. Man entnimmt aus dem angeführten Tafel- 
werk die Werte der zweiten Spalte und berechnet 
damit die Differenzen 1 y;: 
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I; 7 Iy; ? 17%; 
3.54 0,387 95092 

919725 
00 307 14817 + 3972 

315753 3 
3.56 106 30570 + 4062 

Y1isbi + 93 
3.57 115 42261 t- 4155 


907536 
3.58 124 44797 


Nun berechnet man 


0.0? 72398 
0,72398 . 


001 
AY% 0.02669863 . 
A (h BR 2 207 
= [2 Yo 007397. 
(4 1) (/ 2 
| 13 4, {- 0,074 
6 i 


und erhält schließlich. wenn man 


linear interpoliert n 1): sin «a 0,394 60955 
quadratisch .. n 2): sin 0,394 61352 
kubisch e n 3):sina 0.304 61348 


Krst die kubische Interpolation ergibt den genaue 
Wert. In derselben Weise wären die Funktionen 
eos .r, Zin z und Cvj x zu berechnen. 

Die Funktionen tz x und Tg .c werden am ein 
fachsten durch Berechnung der zugehörigen Sinuss« 
und Cosinusse interpoliert. 

Für die Berechnung der weniger häufie vorkom 
inenden Umkehrfunktionen empfiehlt sich folgende:ı 
Voreane: Es sei Y are sin.2 für ein mehrstelliees 
v zu berechnen. Wir wählen ein dem Argument .ı 
‚zunächst liegendes ,, für welches ein Tabellenwert 
y, vorlieet und setzen Y Yı te. Bei der Ent 


wicklung von sin Y ct wird für sin und für 
cos e = | „ gesetzt und daraus die Verbesserung 
rl, 
€ €, 7 €” 
) 2 
vyi—ı 


gereelmet. Dieser Ausdruck wird durch Iteration 
ausgewertet. Als erster Näherungswert wird am 
„weckmäßigsten 


genommen. 

Beispiel: Aus den Tafeln entnehmen wir arcsin 
0.751 — 0,8495752. Für das um ein Tafelintervall 
tiefere «, = 0,750 ist 








Y, = are sin 2, = 0,8480621.. 


Nun berechnet man 


vos Yı = ] 1 — £,* V.H6614378 
und damit 

&, = 0.0?151358. 
Nach der ersten Iteration verbessert sich diese 
Wert in & = 0.02151 316, womit sich YV = Yı + & 
0,8459 5752. übereinstimmend mit dem obigen Wert. 
ergibt. 
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Entsprechend errechnen sich die Korrekturen e für 
lie Umkehrfunktionen 


aärc 608 2; 


Br E N, 1 
& €, ie &, 
2yı1 yi-a, 
Ar Sin: 
r c 1 
Bee, &? 6, : 
2 yVre?+1 Ver-+1 
Ar oje: 
Bin ! N, 
€, i t € : 
2 yı’ 1 vr 1 
aretga 
I L 
er. 
&e,I1— € 
AN 2 I+rz, 
Ar Igu 
€” r 7 
1 + 
‚| 5 ig er 
wenn wieder Y Yı re gesetzt wird, worin 


ein zum Argument 7, gehöriger Tafelwert ist. 
Sehr einfach gestaltet sich die Interpolation bei 
der e-Funktion. 


wobei für ein 


1: e’n j 


E)e 
Mit einer Rechenmaschine mit Rückübertragung ist 
lie Berechnung dieser mehrfachen Produkte sehı 
schnell durchführbar. 
Im Tafelwerk von Hayashi ist im Anhang (S 
274) ein Beispiel der Berechnung von e füı 
a EEE 
r 0.628/31/85307 


.) 


lureheeführt. Wir teilen die Zahl x so ab, wie es 
(lie Striche anzeigen und setzen 


eN,05853 7 1.085307. 


Mit den aus der Tafel entnommenen Werten für 
e,62S und eN,P®3l ist dann 


e» 1°873859 111 - 1033 10048 - TOPS5B1 = 1.87445E0SS 
Mit kubischer Interpolation ergab sich 
e> 1.8744560875 . 


Der wahre Wert von e> ist 1,87445 60875 85. 
Schließlich sei yv=In& zu interpolieren. Es ist 


In(a, +: nz,11+-—)=Iinz, +in(1+5 


£ . 
wenn ds = eesetrt wird. 
r 


Nach Reihenentwieklung des In (1 +6) bis zur 
Iritten Potenz von Ö und Anordnung der Glieder 
nach Art des Hornerschen Schemas erhält man die 
für die praktische Berechnung gut auswertbarv 
Formel 


Inr=Inz, t 1 (05 03 Rh) 


Graz. Hans Ereer. 3% 


Zur Interpolation von Kurvenscharen. 


Ausgehend von bereits vorliegenden Arbeiten, 
welche affine Transformationen zwischen deı 
Funktions- und einer Interpolationsebene b« 
nutzen, wird die Leistungsfähigkeit solchen 
Transformationen besprochen. Ein bereits eı 
probtes, nicht notwendig affines Verfahren lie 
tert oft sünstige re Umre chnunge N. 


I. Einführung 

In der mathematischen Praxis werden häutig 
Kurvenscharen verlangt. deren Einzelkurven nur 
‚dureh umständliche und zeitraubende Rechengänge 
zu ermitteln sind. Zur Herabsetzung des Arbeits 
aufwandes ist es dabei zweckmäßig. nur einige 
wenize Grundkurven nach ihrem Bildungsgzesetz 
„u errechnen, die übrigen Kurven der Schar da 
xeren durch Interpolation zwischen diesen Grund 
kurven aufzusuchen. Zu dieser Aufgabe. insbeson 
dere in ihrer Anwendung auf die Fluzbahnscharen 
von Geschossen und Bomben, sind in den letzten 
„wei Jahren mehrere Beiträge !) ?) ?) *) °) geliefert 
worden. 

Im folgenden wird zunächst auf die von 
R. Sauer?) gemachten Ausführungen zurück 
veeriffen und schließlich ein Verfahren beschrie 
ben. das bereits vor einigen Jahren in der Bomben 
ballistik zur Anwendung webracht wurde ®) 7) ° 
und dem Einzelfall angepaßte, nieht mehr lineare 
Transformationen zwischen der Originalebene und 
ler Interpolationsebene benutzt. 


2. Aufgabe 
Vorzelegt seien für die Funktionen 
r=er(t.‘) Y y(t.i De 2 1 


mit «er Hilfsvariablen ? und dem Parameter / 
einiee diskrete. noch näher zu kennzeichnende 
Einzelwerte, aus denen durch geeignete Interpola 
tion weitere Funktionswerte im Bereich f A 
Ja—r=/, abgeleitet werden sollen. Für diesen 
Interpolationsbereich möge die Darstellung (1) in 
einem reehtwinkligen Koordinatensvstem .r, y ein 
Kurvennetz mit den Kurven / konst. und f 
konst. festlegen. in welehem durch jeden Punkt 
venau eine 4-Kurve und eine #-Kurve geht. 

Der von Sauer benutzte Leitgedanke ist der, die 
Interpolationen in einem Koordinatensystem 


vorzunehmen. dessen Kurvennetz / konst.. / 
konst. für diesen Zweek besonders zeeigenet ist 
und das mit dem .r-. y-System durch affine Trans 
[ormationen verknüpft ist. 


1) H. Athen, Interpolationsverfahren zur Berechnung von 
Fiugbahnseharen und ihrer Veränderung durch Variation 
der ballistischen Grundwerte. Z. angew. Math. Mech. Bd. 1% 
1939), S. 361 bis 371. 

2) H. Athen, Der Bombenwnrf aus dem Sturzflug und 
dem Gleitflug. Luftf. Forschg. Bd. 17 (1940), Heft 7, 
S. 216 bis 220. 

3) R. Sauer, Über Interpolation von Kurvenscharen mit 
Anwendung auf die Bereehnung von Geschoßflugbahnen. 
Z. angew. Math. Mech. Bd. 20 (1940), S. 280 bis 284, 

4) K. Stange, Zur Berechnung einer Flugbahnschar nach 
dem Athenschen Verfahren. Z. angew. Math. Mech. 
Bd. 20 (1949), S. 350 bis 357. 

5) H. Athen. W. Correll, H. Athen mit „Zuschriften an 
den Herausgeber“ in Z. angew. Math. Mech. Bd 21 (191]), 
S. 63 bis 64. 

6) H. Knobloch, Grundzüge der Rallistik des Bomben- 
wurfs. In „Beiträge zur Ballistik und technischen Physik‘, 
hrsg. v. H. Schardin. Leipzig 1998. 

?) H. Knobloch, Näherungsformeln für die Ballistik 
pendelnder Bomben. Z. ges. Schieß- und Sprengstoffw. 
3d. 33 (1938), S. 346 bis 348. 

8) H. Knobloch, Die Ballistik in der Luftwaffe (Bd. 12 
der Luftfahrt-Lehrbücherei). Berlin 1942. 
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3. Die affinen Transformationen 


Eine (für festgehaltenen Parameter 4) aftine 
Transformation zwischen den Systemen r. y und 
=. y läßt sich schreiben in der Form 
(A 74, 


4 E.V b,(4 bh 


e(t,A a 4)» &(t,% a,(4)-n(t,2 


1) (8.7) + b,(2)- nit. 


Bei der Festlegung der zunächst noch frei verfüg 
haren sechs Koeffizienten in (3) lassen sich sechs 
Forderungen erfüllen. welche man an das &» 
System zur Herbeiführung günstiger Interpolations 
bedinzungen stellen kann. 


Im einfachsten Falle leistet ein Ansatz etwa 
der Art 
tt,‘ t 7(t,.% v0 
EN n(t,. 4) 1 
E(t,,d t, n(t,, 7) 0 


zute Dienste, wobei man /,, f, in die 


Nähe der 
Bereichsgrenzen #,. f, und 7, in die Mitte von 
’3 —t, legt. In den Fällen jedoch. wo die durch 
wenige diskrete Parameterwerte 4; (i .®& 
ausgewählten ” Grundfunktionen 


l 


weh: vi; [7] 


Y t.%; 


aus denen durch Interpolationen die Funktionen 
für beliebiges # innerhalb des Bereichs hergeleitet 
werden sollen. aus Differentialgeleichun 
zen gewonnen wurden oder wo zumindest die 
Differentialquotienten zahlenmäßig leicht bestimm 
bar sind. lassen sich mit Vorteil Festlegunzen der 


NE N”: 
. OS en o’S ( 
\bleitungen und _, und auch —,, und 
ot of Of? ( 


ı 
‚ ver- 
t? 


N = 
wenden (z.B. alle bis auf \ gleich Null). 
of 


Die Koeffizienten a und b in Gl. (3) sind nun 


von 4 über die Zwischenfunktionen &(t,. A). 
Ö5 r " - 2 ’ 

(t,.A).... abhäneie. müssen also. da diese Größen 
DE; i 


nur für einzelne Parameterwerte 4: (i 1.2 ,..9) 
errechnet wurden. selbst zunächst über 4 inter 
poliert und tabelliert werden. Diese Notwendigkeit 
stellt den Wert des Verfahrens sehr in Fraee. 
Sauer?) tindet in dem Spezialfall. daß x und y 
rechtwinkliee Koordinaten und die laufende Zeit 
darstellt. die Ableitungen von r und y nach f also 
physikalische Bedeutung haben. einen brauchbaren 


Ausweg dadurch. daß er die zünstiesten Inter 
polationsbedingungen nicht in die Mitte des / 
Intervalls legt, sondern an den Anfang. Diese 


Möglichkeit, die Abhängigkeit der Koeffizienten « 
und 5 vom Parameter analytisch anzugeben. 
liegt aber für den Allgemeinfall nicht vor. Im 
foleenden wird daher ein anderer Weg beschritten. 


{. Erweiterung 

Der Grundgedanke des beschriebenen Verfahrens 
ıst der, den Haupteinfluß der beiden Parameter / 
und 7 getrennt zu erfassen und Funktionen von 
wei Veränderlichen nur zur Korrektur einzufüh 
ren. um diese dann leicht und großmaßstäblich 
interpolieren zu können. 

Häufie kennt man nun bereits. z. B. aus physi 
kalischen Überlegungen, die Abhängigkeit einer 
Funktion von ihren Parametern in mehr oder weni 
ver euter Annäherung. In solchen Fällen wird man 
meist viel einfacher gebaute Transformationen 
zwischen der .r, y-(Original-)Ebene und der £, >» 
(Interpolations-)Ebene verwenden können, die dann 
natürlich auch nicht mehr linear in x und y bzw. 
= und 7 zu sein brauchen. Solche Verfahren wurden 
bereits in der Bombenballistik zur Anwendung ge 


bracht wobei die Art der Transformationen 
aus den ewonnen wurde, die man 
beim Ansatz einer linearen Abhängiekeit der Luft 
widerstandsverzögerung von der Bahngeschwindig 
Näherungsformeln erhält. Beim Bomben 
wurf aus dem Horizontalfluz findet man z. B. 
für die Abhängigkeit der horizontalen Wurfweite X 
und der Fallzeit 7 von der Flugzeug - Eigen 
seschwindiekeit ”, und der Abwurfhöhe AH die 
Transformation 


Ergebnissen & 


keit als 


\uf dem gleichen Wege erhält man beim Wurf aus 
dem Sturzflug für die Abhängiekeit der Wurf 
weite X und der Fallzeit 7 von Flugzeugr-Eiren 
xeschwindigekeit ®,. Fallhöhe H und Bahnneigungs 
winkel # die Transformation 





T 1 
4 
V2-g9H (1,01 +Hm)+ (v, sin A)? — r,sin A > 
X: eos/T(1 TE 
mit der Umkehrung 
f Y 
& — +11 A —— 
T | U, COS %; 2 
15) 
u er v,sin2 T — 1.01-H | 


wobei die &. 7 nach (6) in Abhängigkeit von v,. H 
und 4 zu tabellieren sind. Die Kurvennetze deı 
Interpolationsebene 5, 7 liefern dann zusammen 
mit der Rechenvorschrift (5) die Wurftafel der be 
treffenden Bombe für beliebige Flugzeug - Eigen 
reschwindigkeit. beliebige Abwurfhöhe und be 
liebigen Bahnwinkel im Abwurfpunkt. 

In Bild 1 und 2 sind als Beispiele Interpolations 
netze für den Bombenwurf aus dem Horizontalflug 





Pa] 
m’S) 
70 *34, 
29 
7 
| 
3 26 un 46.7 
£ 
& 
Bild 1. Interpolationsebene zum Formelsystem (4). Wurf 


tafel für den Bombenwurf aus dem Horizontalfing. 
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Bild 2. Interpolationsebene zum Formelsystem (5). Wurf 
tafel für den Bombenwurf mit der Flugzeug-Eigenge 
schwindigkeit v, 100 m/s. 


und dem Sturztlug wiedergegeben. Durch Differen 
tiation der zugehörigen Formeln (4) und (5) ge 
winnt man einen Überblick über die Interpolations- 
senauiekeit. Danach liert die Unsicherheit in 
der Fallzeit unter 0.05 s und in der Wurfweite 
unter 1°/o. wenn S und n in Bild 1 mit einer Ge 
nauigkeit von 0.6- 10° s/m bzw. 0,8.10°° s/m, in 


Bild 2 mit 0.3.10 s”1 bzw. 0.9.1076 m 1 Ge 
nauigekeit interpoliert werden. Die Unterlagen 
zu Bild 1 und 2 wurden für = 0° dem Aufsatz ® 


entnommen. 


»), H. Knobloch, Beitrag zur Lösung des außenballisti- 
schen Hauptproblems. Z. ges. Schieß- und Sprengstoffw. 
Bd. 35 (1949), S. 245 bis 248 und S. 270 bis 279. 


5. Erzebnis 


Es ergibt sich somit für die Praxis folgendes 


Bild: 


Die Funktion (a, P....) der unabhängigen Ver 
änderlichen «. f.... liegt für verschiedene Einzel 


werte «a; Pj.... errechnet vor: sie soll in «. ß.... 
interpoliert werden. 

Durch physikalische Überlegungen oder durch 
probeweises Aufzeichnen ermittelt man eine analy 
tisch einfache oder tabellarisch genau wegebene 


z. B. durch sinus. Quadratwurzeln .. .) Näherungs 
funktion «(a, ß....) und ergänzt sie dureh ein 


Korrekturglied 


) j} Pr 
FE er LreB.:; Re N re ı ı 


Da die Zahlenwerte dieses Korrekturgliedes 
sind gegen die Originalfunktion (a, P....). so 
wird man 5 meist ohne weiteres gut in «, ß.... 
interpolieren können, besonders dann. wenn die 
\äherung «ala, ß,...) so gewählt wurde, daß = 
etwa in der Mitte des (resamtbereichs verschwin 
det. Die Funktion :r(a, ß....) liegt also in der 
Form 


s 
Fr 


1er, +) al RE le 
vor, wobei jetzt nur die Korrekturfunktion = 
aufzutragen und zu tabellieren ist. während 
ala,3,...) analytisch oder durch engmaschige Ta 
bellen gezeben ist. Die Interpolationsgenauigkeit 
für die Originalfunktion r ist dann, absolut ze 
nommen, die gleiche wie für das Korrekturglied S. 
welches jedoch gegenüber der Originalfunktion nur 
einen ganz wesentlich kleineren Funktionsbereich 
ausfüllt und sich daher graphisch gut auftragen 
läßt. In manchen Fällen wird an Stelle von (8 
die Form 


2882...) l&B,..)- TE nla, B,. ..) 


oder ein ähnlieh sinnvoll aufzebauter Ansatz vor 
teilhafter sein. 


Berlin. H. Knobloch. 39 


BUCHBESPRECHUNGEN 


\bhandlungen der Internationalen 
Vereinieung für Brückenbau und 
Hochbau, Bd. 6. 1940/41. Zürich 1942, heraus- 
eereben vom Generalsekretariat in Zürich. Preis 
brosch. 15 M. . 

Der 6. Band enthält 10 Abhandlungen in deut 
scher, 2 in französischer und 4 in englischer 
Sprache, ein Teil hiervon gehört zu den Arbeits 
eebieten des für 1940 in Aussicht genommenen 
Kongresses. Der Inhalt des Bandes ist in mehr- 
facher Hinsicht von Bedeutung. Zunächst durch 
(tie Einbeziehung von Ansätzen und Anschauungen 
der Plastizitätstherorie in die Baustatik. Zu 
nennen sind: P. P. Bijlaard über die Theorie 
(ler örtlichen plastischen Formänderungen und der 
plastischen Stabilität: J. Fritsche über die 
Balkenbiegung im plastischen Bereich. — Für die 
Statik ist hervorzuheben die Arbeit von E. Am- 
stutz, Erweiterte Theorie des Seilpolygons zur 
Berechnung von Hängebrücken: darin wird gezeigt. 
daß sich die Differentialgeleichung zweiter Ordnung 
für das Biegemoment im Versteifungsbalken in 
zwei simultane Gleichungen erster Ordnung auf- 
spalten läßt, aus denen die erweiterte 
„ideelle* Querkraft und das Biegemoment durch 
Seilpolygone ermittelt werden kann. Die Ver- 
steifungsträger von Hängebrücken haben als Gleich- 


nis das gespannte, elastisch unterstützte Seil. die 
gewöhnlichen Biegeträger das freie gespannte Seil. 
Mit dieser steht die Arbeit von F. Stüssi, 
Statik der Seile in engem Zusammenhang. — Der 
Elastizitätstheorie gehören die Arbeiten von 
7. Bazant über diekwandige zylindrische Rohre 
und von Ch. Massonet über die Theorie des 
Kniekens durch Verdrehung an. Die übrigen 
Arbeiten befassen sich mit anderen Knickproble 
men und mit Aufgaben aus der Baustatik. 
Karlsruhe i. B. Th. Pöschl. 429. 


Dr. habil. ALBERT KOCHENDÖRFER, Dozent 
für Physik a. d. Techn. Hoehschule Stuttgart, Pla 
stische Eigenschaften von Kristal 
ien und metallischen Werkstoffen. 
Reine und angewandte Metallkunde in Einzeldar 
stellungen, herauszegeben von W. Köster. Nr. 7.) 
XII + 312 S. m. 91 Abb. Berlin 1941. Verlag Julius 
Springer. Preis geb. 2850 M. 

Der physikalische Vorgang für die Entstehung 
bleibender Formänderungen in kristallischen Werk 
stoffen ist die plastische Gleitung. Ihre Unter 
suchung rückt immer mehr in den Mittelpunkt der 
praktischen Metallkunde, da sie zugleich auch die 
physikalische Grundlage für die zahlreichen tech 
nischen Verformungsmöglichkeiten bei der Ver 
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arbeitung der Werkstoffe darstellt. Die einfachen 
umkehrbar-ein-eindeutigen jeziehungen zwischen 
den Spannungs- und Verzerrungsgrößen, die im 
elastischen Gebiet Geltung haben. werden bei Über 
schreitune einer oft sehr unscharf vorhande 
nen Grenze von einer Vielheit von Ersceheinun 
zen abgelöst. die heute sowohl qualitativ, als auch 
quantitativ nur sehr unvollkommen bekannt sind. 
Ihre Erforschung hat in der Technik zur Aus 
bildung eines uneeheuren 
\usmaben geführt, aber immer stärker wird das 
Bedürfnis nach wesicherten theoretischen Grund 
lagen für das ganze Gebäude empfunden. Auch 
der Verfasser des vorliegenden Buches ist sich 
Bedeutung dieser Aufrabe und auch der Schwieri: 
keiten bewußt, die seiner Lösung entgegenstehen. 
\ls einfachstes Gebilde dient der Einkristall bei 
homogener und niehthomogener Verzerrung: die 
Bezeichnungsweise erfolet im wesentlichen nach 
dem Vorbilde von G. Tammann. sie wird aber 
nach Ansicht des Ref. den charakteristischen Un 
stetigkeiten «des plastischen Gleitvorganges nicht 
erecht. in («denen sich «(lie Mosaikstruktur der 
Stoffe äußert. Aus dem Einkristall wird das Veı 
halten des vielkristallinen Werkstoffes gewonnen. 
Die Einteilune des Buches ist demnach (die fol 
eende: 1. Einsinnige Verformung. 
II. homogene Wechselverformung. III. inhomogene 
Verformune von Kristallen. IV. einsinnige Ver 
formung, V. Wechselbeanspruchung metallischer 
Werkstoffe. In VI. werden die mathematischen 
Zusätze eeeeben. 

In dem Buche sind vielfach eigene, experimen 
telle Untersuchungen des eereben. 
außerdem ist ein umfangreiches Schrifttum. vor 
allem im Hinblick auf die erhaltenen Versuchs 
ergebnisse verarbeitet. An Hilfsmittel ler Theorie 
vird eigentlich nur der einachsige Spannungszu- 
stand und auch bei der Biegunze und Torsion 
werden nur die elementarsten Ansätze verwendet. 
Dagegen sind die physikalischen Theorien. die auf 
diesem Gebiete zur Ausbildune welaneten. die 
Mosaikstruktur, Verfestigung, die  Versetzungen 
u. del. in großer Ausführlichkeit dargestellt: die 
Bedeutun«e der Thermodynamik für die Theorie 
tritt jedoch nach Ansicht des Ref. zu sehr in den 
Hintergrund. Wer die Sachlage kennt, wird ab 
schließende Ergebnisse nicht erwarten. und größer 
als das bisher Erreichte ist das. was in der Zu 
kunft zu tun übrig bleibt. 


Versuchswesens von 





homogen® 


Verfassers 


Karlsruhe i. B. Th. Pöschl. 4830. 


Dr. E. BOLLE, Einführungindieinnere 
Ballistik. Nach dem Vorlesungsmanuskript und 
den Arbeiten von E. Boll@ bearbeitet von Dr. 6. 
Seitz. VI+139 S. m. 48 Abb. Braunschweig 
1942, Verlage Fr. Viewer & Sohn. Preis eeb. 10 M. 

In dem vorliegenden Werk sind die Veröffent 
liehungen, Vorträge und Vorlesungen des verstor 
benen Oberregierungsrats an der Chemisch-Techni- 
schen Reichsanstalt Dr. Erwin Bolle zusammen 
sefaßt und durch neuere Erkenntnisse in seinem 
Sinne ergänzt worden. 

In sechs Kapiteln werden u. a. behandelt: Meß 
methoden der inneren Ballistik. empirische Formeln 
zur Berechnung der innerballistischen Werte, Ver 
brennungsweise des Pulvers, Lösung des innerbal 
listischen Hauptproblems durch geschlossene For- 
melsysteme, Einpreßwiderstand und Geschoßrei 
bune, Bereehnune der Geschützrohre (Festigekeits 
lehre). Explosion und Detonation. In allen Kapiteln 
finden sich zahlreiche Zahlenbeispiele zur Erläute- 
rung und zur kritischen Beurteilung der theoreti 
schen Methoden. In dieser Beziehung kann (die 
Darstellung der empirischen Formeln und der 
mathematischen Lösungsmethoden zur Berechnung 
der innerballistischen Werte als besonders gelungen 


bezeichnet werden. Außer den bereits von Boll« 
behandelten Formelsystemen von Hevdenreich, 
Leduc und Ledue-Bolle, sowie den Lösun 
(sossot-Liouville, Charbonnier und 
wird die neue Lösune von Langeweiler einerehend 
behandelt. 

Das kleine Werk wendet sich an alle, die das 
Studium der inneren Ballistik beruflich betreiben. 
aber nach den Beeleitworten von Prof. Dr. Ri 
marski auch an diejenigen, die über die inner: 
Ballistik und die Explosionsvorgänge unterrichtet 
sein wollen. um ihre eigenen Kenntnisse und ih 
Wissen zu vertiefen, 
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Dr. habil. WILHELM BECKER, Dozent an deı 
Universität Wien. Sterne ınld Stern 
systeme. Wissenschaftliche Forschungsberichte. 
Naturwissenschaftliche Reihe. herausgereben von 
Dr. Raphael Ed. Liesegane, Bd. 55.) XII + 394 S 
m. 94 Abb. Dresden und Leipzig 1942, Verlag Theo 
lor Steinkopff. Preis geb. 30 M 


Das vorliegende Buch soil. wie (er Verfasser im 
Vorwort schreibt. eine Zwischenstellung zwischen 
einem rein wissenschaftlichen und einem populären 
Buch einnehmen. Gleichzeitie soll dem Leser dure! 
ausführliche Angaben von Zahlenmaterial und Lite 
raturquellen ein übersichtliches Handbuch für das 
dargestellte Gebiet der Astronomie gewreben werden. 
Dieses weitzeesteckte Ziel hat der Verfasser in An 
betracht der Fülle des Stoffes auf verhältnismäßig 
eneem Raum erreicht. Zwar ist das Buch als Gan 
zes weit mehr als wissenschaftlich wie als populär 
zu bezeichnen: trotzdem wird auch der allgemein 
naturwissenschaftlich interessierte und vorzebildete 
Leser mit Nutzen und Freude sich daraus über den 
eerenwärtigen Stand der Erforschung des Weltalls 
orientieren können. 

Der Hauptteil des Buches behandelt las Milel 
straßensystem. Dem Milchstraßensystem wird dann 
in einem zweiten Teil die Behandlung der außeı 
galaktischen Nebel zegenübergestellt. die ja. soviel 
wir heute wissen, dem Milchstraßensystem koortli 
niert sind. 

Das Buch stellt. auch was das Zahlenmaterial be 
trifft, einen nach dem heutigen Stand der Foı 
schung zuverlässigen Führer durch das weite Ge 
biet der Stellarastronomie dar. Derjenige, der Ge 
legenheit hatte. die Entwicklung dieses Zweigs deı 
\stronomie während der letzten 20 Jahre zu veı 
folgen. wird stark beeindruckt sein von der Fülle 
des Beobachtungsmaterials. das in diesen Jahren 
eesammelt wurde und in dieser kurzen Zeit zu einer 
so erundlegenden Umeestaltung unserer Auffassung 
vom Bau des Sternsvstems und der Welt der außer 
oalaktischen Nebel geführt hat. Man vergleiche 
etwa das hier entwickelte Weltbild mit dem Welt 
bild. das die 1920 führenden Astronomen in (der 
6. Auflage (des Neweomb-Engelmann dargestellt 
haben! 

Für den Leser dieser Zeitschrift scheint es mir 
wesentlich festzustellen. daß in dem vorliegenden 
Buch theoretische Entwicklungen nur zestreift weı 
den. Sie sind ja auch, jedenfalls soweit es den 
\ufbau des Sternsystems betrifft. seit Seeliger und 
Schwarzschild ganz verlassen worden, und haben 
mehr oder weniger rein empirischen Methoden Platz 
vemacht. \ndererseits sind die Beobachtunges 
ergebnisse für die Aufstellung einer dynamischen 
Theorie des Milehstraßensystems noch viel zu dürf 
tie, als daß es sieh lohnen würde, eine ins einzelne 
eehende Theorie zu entwickeln. Kosmologische 
Theorien werden in dem Kapitel über die räumliche 
Verteilune der außergealaktischen Nebel nur kurz 
berührt. 4 

Da sich das Buch an einen weiteren naturwissen 
schaftlieh interessierten Leserkreis wendet, hätte 
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in dem Kapitel über die Bewegung und Entfernun 
een der Fixsterne etwas näher auf die großen 
Schwieriekeiten eineezeaneen werden können. (die 
für den Astronomen in der Festlegung eines ze 
eieneten Koordinationssystems bestehen. Mißt er 
doch die Positionen und Bewegungen der Sterne 
von der beweeten Erde aus. (leren Rotationsachse 
und Bahnebene infolge der Anziehung von Sonne, 
Mond und Planeten eine komplizierte Bewerung 
m Raum ausführt. Die Frage. wie gelangt der 
Astronom überhaupt in aufeinanderfolgenden Nä 
heruneen zur Festlegung eines Koordinaten 
systems. das als Inertialsystem im Sinne der klassi 
sehen Mechanik eelten kann. würde sicher auch 
manchen Leser dieser Zeitsehrift interessiert 
haben. 

Dem Buch. (das mit zahlreichen zuten und ein 
dlrucksvollen Abbildungen ausgestattet ist, kann 
eine weite Verbreitung gewünscht werden. 

Potsdam. P.tenBruggencate. 422 

Ferner sind bei der Schriftleitung folgende Bücher 
eingegangen (ausführliche Besprechung bleibt vor 
behalten): 

Dr.-Ing. GUSTAV NIEMANN, Prof. a. d. Techn. 
Hochschule Braunschweig, Schneckentriebe 
mit flüssiger Reibung, Abhängirkeit 
der übertragebaren Leistung und des 
Reibwertes von Zahnform. Abmes 
sune, Prehzahl und Sehmierzähig 
keit. mit einem Anhane mathematischer Ab 
leitungen von Prof. Dr.Ing. Constantin 
Weber. Dresden. (VDI-Forsehunesheft 412. Bei 


lage zu „Forschung auf dem Gebiete des Ingenieur 
wesens“, Ausgabe B, Bd. 13, Januar/Februar 1942.) 
29 Ss. m. 37 Bildern und 3 Tafeln. Berlin 1942, 
VDI-Verlage G.m.b.H. Preis broseh. 5 M. für VDI 
Mitelieder 4,50 M. 


Dr.-Ing. E. MELAN, o. Prof. a. d. Techn. Hoch 
schule Wien, Dipl.-Ing. R. SCHINDLER, Obering. 
der Waagner-Biro A.-G.. Wien. Die genaue Be 
rechnung von Trägerrosten, mit 39 Ta 
feln für die praktische Anwendung ausgearbeitet 
von Dipl.-Ing. R. Rothmayer. V + 134 SS, m. 
101 Abb. u. 9 Zahlentafeln. Wien 1942, Springen 
Verlag. Preis brosch. 21.10 M. 

Prof. Dr. A. WITTING, Repetitorium und 
\ufgabensammlung zur Differential 
rechnung. 2. neubearb. Aufl. (Sammlung Gö 
schen Bd. 146). 144 S. m. 57 Fig. u. 415 Beispielen 
und Aufgaben. Berlin 1942, Verlag Walter de 
(Gruyter & Co. Preis geb. 1.62 M. 


Dr. phil. BRUNO TIEDEMANN, Über Boden 
untersuchuneen bei Entwurf und 
\usführune von Inzenieurbauten. 
2. verb. Aufl. 40 S. m. 22 Abb. und 5 Tafeln. Ber 
In 1942. Verlae Wilhelm Ernst & Sohn. Preis 
brosch. 2 M. 


Dr. H. GREINACHER, o. Prof. a. d. Universität 
3ern, Ergänzungen zur Experimental 
physik, einführende exakte Behand 
lung physikalischer Aufgaben, Fra 
gen und Probleme. X+ 181 S. m. 79 Abb. 
Wien 1942. Springer-Verlag. Preis brosch. 6.60 M. 


NACHRICHTEN 


Ludwig Schiller. 


\m 11. Dezember 1942 vollendete Prof. Ludwig 
Schiller. Vorstand der Abteilung für angewanidte 
Mechanik und Thermodynamik beim Physikalischen 
Institut der Universität Leipzig sein 60. Lebensjahr. 
Das in München beeonnene Studium schloß er 1911 
mit der Doktorarbeit über die Änderung der Dielek- 
trizitätskonstanten des Kautschuks bei Zur senk 
reeht zu den Kraftlinien ab. Im Weltkrieg erwarb 
er sich hohe Auszeiehnuneen, u. a. EK. 1 und 2. 

Für «das wissenschaftliche Lebenswerk Schillers 
bestimmend war ein Aufenthalt im Prandtlschen 
Institut in Göttingen 1919 bis 1920. Hier entstanden 
die in seiner Habilitationsschrift ') niedergelegten 
Untersuchungen über das Anlaufproblem der Rohr 
strömune, hierin findet sich erstmalig die Angabe 
ler unteren Schranke (der kritischen Revnoldsschen 
Zahl bei beliebige eroßen Störungen (Re — 1160). 
Von hier ab bestimmte die Frage nach den ver 
wiekelten Voreäneen beim Übergang der Laminar 
strömung zur Turbulenz die eine Riehtung der Schil 
lerschen Arbeiten, und man darf wohl sagen, dab 
er diesen Weeen mit eanz besonderer Liebe nach 
eine. Die zahlenmäßige Erfassung des Betrags der 
zum Strömunesumschlag erforderlichen Störungs 
enereie ebnete den Were zu vielen sehönen Erkennt 
I issen. 

Von den sehr zahlreichen übrigen Arbeiten Schil 
ers?) seien hier nur die Bestimmung der Wider 
stands- und Wärmeübergangszesetze bei hohen 
RRevnolisschen Zahlen in Rohren unter verschie 
denen Betriebsbedineungen. sowie Wärmeüber 
traeunesversuche bei freier Konvektion und Lei 


1, Forsch.-Arb. Ing.-Wes. Heft 248 (1922); auch Z. angew. 


Math. Mech. Bd. 2 (1922), S. 96. 
?\ vgl. z.B. Phys. Z. Bd. 26 (19251, 8. 566: Forseh. Ing. 
Wes. Bd. 10 (1939), S. 246 


tung. ferner grundlegende Überlegungen zu Zähig 
keitsmessungen und aerodynamische Arbeiten er 
wähnt. Als Herausgeber des vierteilizen 4. Bandes 
les Handbuches der Experimentalphysik hat er uns 
unter Mitwirkune namhaftester Fachkollegen ein 
Fundamentalwerk für «die Forschungsarbeit ve 
schenkt. 

Neben seiner Tätiekeit als Hochschullehrer und 
Forscher steht Schiller seit 1928 der Versuchsstelle 
für Hauswirtschaft des Deutschen Frauenwerks als 
wissenschaftlicher Oberleiter vor. 

Das Kennzeichnende der Schillerschen Forschungs 
art ist die ins Kleinste gehende Untersuchung der 
letzten Einzelheit. und keine Erscheinung läßt ihm 
Ruhe. deren Deutung noch nicht restlos geklärt er 
scheint. Selbst als Experimentator mit den feinsten 
Tücken der Versuchstechnik vertraut. erzieht eı 
seine Schüler zu eiserner Selbstkritik und was 
vielleicht das Wertvollste ist. das seine Schüler von 
ihm in ihren Lebensweg mitnehmen, zu jenen 
Unermüdlichkeit, die befähigt, im Kampf mit den 
widerstrebenden Geistern der Versuchsanordnung 
Sieeer zu bleiben. Und wie im Versuch. so eilt 
auch der Deutung der Ergebnisse oft stundenlange 
Diskussion mit immer neuen intuitiv gegebenen An 
reeungen. 

Wir wünschen dem Geburtstagskind. daß noch 
recht viele Schüler aus seinem reichen wissen 
schaftlichen und versuchstechnischen Erfahrungs 
schatz schöpfen mögen. 


\achen. \.Naumann. 431 


Persönliches. 

Prof. Dr. M. La _qA | ly m Prof. der anee 
wandten Mathematik an der Techn. Hochschule 
Dresden, wurde von seinen amtlichen Verpflich 
tungen entbunden (emeritiert 
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In Weuaujlage liegt vee: 


Der Ultraschall 1'wissensaat ung Teanix 


Von Dr. Ludwig Bergmann, o. Professor für Physik 
und Direktor des Physikalischen Instituts der technischen Hochschule Breslau. 


Dritte, völlig überarbeitete und erweiterte Auflage. 1942. 
DIN A 5, 457 Seiten mit 269 Bildern und 47 Tafeln. 


Gebunden 25.— RM, für VDI-Mitglieder 22.50 RM. 


Die außerordentlich günstige Aufnahme dieses Ultraschallwerkes durch die Fachwelt hat bereits 
drei Jahre nach dem Erscheinen der zweiten Auflage eine Neuauflage notwendig gemacht. 
Die neuen Erkenntnisse und Beobachtungen auf dem Gebiete des Ultraschalles haben zu einer 
völligen Überarbeitung des ganzen Stoffes und zu zahlreichen Ergänzungen und Verbesserungen 
geführt; alle Abschnitte, insbesondere die über die praktische Anwendung des Ultraschalles, wurden 
wesentlich erweitert. So stieg der Umfang des Werkes um rund 90 Seiten auf fast 460 Seiten, 
während die Zahl der Abbildungen um über 40, die Schrifttumsstellen sogar um 370 vermehrt wurden. 


Der „Bergmann“ stellt in seiner dritten Auflage zur Zeit die umfassendste Darstellung auf dem 
Gebiete des Ultraschalles dar und wird dem Wissenschaftler und dem praktisch arbeitenden Ingenieur 
ein willkommenes Hilfsmittel und eine Quelle weiterer Anregungen sein. Infolge der Fortschritte 
der letzten Jahre ist die Neuauflage auch für den von Nutzen, der eine frühere Auflage besitzt. 
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Neue VDI-Forschungshefte 


Heft 412 Schneckentriebe mit flüssiger Reibung. Abhängigkeit der übertrag- 
baren Leistung und des Reibwertes von Zahnform, Abmessung, Drehzahl und 
Schmierzähigkeit. 

Von @. Niemann VDI. 
Mit einem Anhang mathematischer Ableitungen. Von C. Weber. 
DIN A 4, 31 Seiten mit 37 Bildern und 3 Zahlentafeln. 1942. 


Heft 413 Nichtstationäre Wärmeübertragung in Gasen, insbesondere in 
Kolbenmaschinen. 


Von Dr.-Ing. habil. P. Pfriem VDI, Berlin. 

Aus der deutschen Versuchsanstalt für Luftfahrt (DVL), Berlin- Adlershof, In- 
stitut für motorische Arbeitsverfahren und Thermodynamik. 

DIN A 4, 31 Seiten mit 9 Bildern. 1942. 


Heft 414 Gesetzmäßigkeiten der freien Turbulenz. 
Von Dr. phil. H. Reichardt VDI, Göttingen. 


Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung, Göttingen. 
DIN A 4, 24 Seiten mit 15 Bildern. 1942. 


Heft 415 Der Wärme- und Stoffaustausch im Trocknungsgut. Die analytische 
und graphische Behandlung der Trocknung poriger hygroskopischer Güter. 
Von Prof. Dr.-Ing. O0. Krischer, Darmstadt. 
DIN A 4, 24 Seiten mit 18 Bildern und 5 Zahlentafeln. 1942. 


Heft 416 Die Berechnung des Wärmeübergangs in der laminaren Grenz- 
schicht umströmter Körper. 
Von Dozent Dr.-Ing. habil. E. Eckert VDI, Braunschweig. 
Aus dem Institut für Motorenforschung der Luftfahrtforschungsanstalt Hermann 
Göring, Braunschweig. 
DIN A 4, 26 Seiten mit 29 Bildern, 3 Arbeitsblättern und 12 Zahlentafeln. 1942. 


Heft 417 Die Notlaufeigenschaften der Gleitlagermetalle in Maschinen der 
Feinmechanik. 
Von Dr.-Ing. H. Lüpfert VDI, Stuttgart. 
DIN A 4, 32 Seiten mit 25 Bildern und 3 Tafeln. 1942, 


Heftpreis 5.— RM., für VDI-Mitglieder 4.50 RM. 
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